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ABSTRACT 

’Iliis  report  presents  a computer  package  designed  to  facilitate  the 
assessment  and  use  of  a decision  maker's  utility  function  for  multiple 
objectives.  Tlie  package  provides  routines  for  (1)  specifying  the  decision 
maker's  pi'eferenccs  over  multiple  criteria,  (2)  treating  uncertainty  in 
the  consequences  resulting  from  a decision,  (3)  ranking  alternative  courses 
of  action  in  order  of  preference,  and  (4)  studying  the  effects  changes  in 
preferences  or  uncertainty  estimates  may  have  upon  the  ranking  of  alterna- 
tives. Tlie  routines  are  designed  to  be  applicable  in  a variety  of  prob- 
lem contexts. 

ITie  paper  is  organized  as  follows.  The  decision  analysis  approach 
which  provides  the  theoretical  basis  for  the  program  is  summarized.  This 
is  followed  by  a description  of  existing  methods  for  multiattribute  utility 
function  assessment  and  use.  Tlien  the  computer  package  is  presented  and 
compared  with  the  aforementioned  methods.  Applications  of  the  package  to 
several  problems  are  illustrated  and  areas  for  future  improvement  and  re- 
search are  suggested. 
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1.  INTRODUCTION 


Many  decision-making  problems  are  characterized  by  two 
sources  of  complexity.  First,  there  are  multiple  objectives 
on  the  basis  of  which  the  de«jision  should  be  made.  In  weigh- 
ing alternative  actions,  the  decision  maker  must  consider  the 
tradeoffs  between  the  degree  of  achievement  in  one  objective 
and  the  degree  of  achie>^ement  in  others.  Second,  there  is 
often  uncertainty  about  the  consequences  which  will  result 
from  any  particular  action. 

Because  of  these  complexities,  there  is  a need  for  a 
formal  approach  to  help  in  evaJuating  alternatives.  Decision 
analysis  is  an  appiccch  whx.:h  explicitly  addresses  the  multi- 
ple objective  and  uncertainty  issues.  The  theoretical  basis 
for  this  is  well  established.  However,  many  practical  prob- 
lems arise  when  one  tries  to  apply  decision  analysis  in  parti- 
cular situations.  This  thesis  describes  a computer  package 
for  overcoming  some  of  these  difficulties. 

1.1 The  Decision  Analysis  Approach 

Raiffa  (14)  discusses  ne  philosophy  and  techniques  of 
decision  analysis  in  detail.  We  can  think  of  the  decision 
analysis  approach  as  consisting  of  four  steps: 

1.  structuring  the  problem, 

2.  quantifying  the  uncertainties  involved. 


10- 


3.  quantifying  the  decision  maker's  preferences,  and 

4.  combining  the  first  three  steps  to  evaluate  the 
alternatives . 

Structuring  includes  identifying  the  decision  maker  and 
the  problem  objectives.  Measures  of  effectiveness  (attributes) 
indicati'  j the  degree  to  which  each  objective  is  achieved  are 
also  fo:  lulated.  Let  us  designate  our  set  of  attributes  as 
X^/***#^  and  use  to  indicate  a specific  amount  of  attribute 
X^.  For  example,  X^  may  be  profit  in  1975  measured  in  thou- 
sands of  dollars  and  x^  may  be  188.  A consequence  will  be 
denoted  by  x = (x^,X2r . . • ,x^^)  and  indicates  the  level  x^  of 
each  attribute  which  results  given  that  consequence. 

Quantifying  uncertainties  involves  describing  the 
uncertainty  in  the  possible  consequences  of  any  alternrtive. 

For  each  alternative  A^ , a probability  distribution  p^  (x) 
indicating  which  consequences  might  occur  and  their  likelihood 
of  occurrence  is  required.  The  p^  may  be  derived  by  means  of 
some  analytical  or  simulation  model  or  by  subjective  assess- 
ments. 

Quantifying  preferences  means  ^“^ssessing  the  decision 
maker's  utility  function  u(x)  = u (x^/X^/ . . . which  assigns 
a number  to  each  of  the  possible  consequences.  This  function 
is  called  a multiattribute  utility  function  and  will  be 
referred  to  by  the  mnemonic  MUF.  A MUF  has  two  properties 
which  make  it  useful  in  addressing  the  issues  of  uncertainty 
and  tradeoffs  between  objectives.  These  properties  are: 
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1.  u(x*)  > u(x”)  if  and  only  if  x*  is  preferred  to  x" 
and 

2.  in  situations  with  uncertainty,  the  expected  value 
of  u is  the  appropriate  guide  for  making  decisions, 
i.e.,  the  alternative  with  the  highest  expected 
value  is  the  most  preferred. 

This  second  property  follows  from  the  axioms  of  ration- 
al behavior  postulated  first  in  von  Neumann  and  Morgenstern [18 ]. 

Evaluating  alternatives  involves  calculating  the 
expected  utility  for  each  of  the  alternatives  using  the  and 
u from  the  previous  steps.  Various  parameters  of  the  probabi- 
lity distributions  and  the  utility  function  can  be  varied  to 
see  how  these  affect  the  expected  utility  of  the  alternatives, 
i.e.,  how  "sensitive"  the  results  are  to  changes  in  the  para- 
meters . 


1.2 Statement  of  the  Problem 

A major  practical  problem  arises  when  one  tries  to 
obtain  a MUF  that  is. "tractable"  yet  appropriate  for  a parti- 
cular situation.  The  general  approach  has  been  to  postulate 
assumptions  about  the  decision  maker's  preferences  and  derive 
the  restrictions  they  place  on  the  functional  form  for  u.  Then, 
for  any  specific  problem,  the  adequacy  of  the  assumptions  must 
be  verified  and  the  parameters  for  the  utility  function 
assessed  and  checked  for  internal  consistency.  Ideally,  the 
functional  form  of  the  MUF  would  have  the  following  properties : 
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1.  be  general  enough  to  apply  to  many  real  problems, 


2»  require  a minimal  number  of  assessment  questions 


to  be  asked  of  the  decision  maker, 


3,  require  assessments  which  are  reasonable  for  a 


decision  maker  to  consider,  and 


4.  be  easy  to  use  in  evaluating  alternatives  and 


conducting  sensitivity  analysis  with  respect  to 


various  parameters. 


Even  with  a convenient  functional  form  for  the  MUF, 


the  nature  and  magnitude  of  a problem  can  make  the  assessment. 


bookkeeping,  and  use  of  quantitative  preference  information  a 


formidable  task.  The  computer  package  described  in  this 


thesis  is  designed  to  handle  this  task  for  a variety  of  prob- 


lem contexts. 


Organization  of  the  Thesis 


Chapter  2 summarizes  the  theoretical  development  of 


the  functional  forms  for  MUF’s  upon  which  the  computer  package 


is  based.  Chapter  3 discusses  exijting  methods  for  assessing 


and  using  MUF's  and  their  difficulties.  Chapter  4 describes 


the  computer  package  and  the  manner  in  which  it  alleviates  the 


difficulties  mentioned  in  Chapter  3.  Chapter  5 presents 


several  applications  of  the  package  to  different  problems 


illustrating  the  use  of  the  various  package  routines.  Chapter 


6 discusses  suggestions  for  improving  the  package  and  for 


■ 
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future  research.  Chapter  7 contains  a summary  and  conclusions 
of  the  thesis. 

Five  appendices  contain  detailed  information  concerning 
understanding  and  use  of  the  computer  package.  Appendix  A is 
a concise  summary  of  the  package  commands.  Appendix  B is  a 
listing  of  the  program.  Appendix  C describes  some  of  the  algo- 
rithms used  in  several  of  the  package  routines.  Appendix  D 
contains  a discussion  of  the  overall  program  design.  Appendix 
E explores  the  tradeoff  properties  among  the  attributes  implied 
by  the  functional  forms  used  for  the  multiattribute  utility 
function.  It  serves  to  explain  the  design  and  use  of  some  of 
the  package  routines. 
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2*  THE  ADDITIVE  AND  MULTIPLICATIVE  UTILITY  FUNCTIONS 

This  chapter  states  the  conditions  which  imply  that  a 
MUF  is  either  additive  or  multiplicative.  None  of  the  condi- 
tions require  the  decision  maker  to  consider  preference  trade- 
offs between  more  than  two  attributes  simultaneously  or  to 
consider  lotteries  (specifying  various  x and  the  probabilities 
of  receiving  them)  with  the  level  of  more  than  one  attribute 
being  varied.  Furthermore,  the  assessments  needed  to  specify 
an  n-attribute  utility  function  are  n single-attribute  utility 
functions  and  n scaling  constants.  Some  properties  of  these 
forms  are  discussed  as  well  as  their  applicability  to  different 
classes  of  problems. 

2.1 The  Basic  Assumptions 

The  two  basic  assumptions  which  we  use  for  both  addi- 
tive and  multiplicative  utility  functions  are  referred  to  as 
preferential  independence  and  utility  independence.  These  are 
defined  as  follows: 

Preferential  Independence:  The  pair  of  attributes 

(Xj^,X2>  is  preferentially  independent  of  the  other  attributes 
(X^,...,X^)  if  preferences  among  (Xj^,X2>  pairs  given  that 

are  held  fixed,  do  not  depend  on  the  level  where 
(X2#...#X^)  are  fixed. 

Preferential  independence  implies  that  the  tradeoffs 
between  attributes  and  X2  do  not  depend  on  X^^.-.^X^. 
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Utility  Independence;  The  attribute  is  utility 
independent  of  the  other  attributes  if  preferences 

among  lotteries  over  X^  (i.e.,  lotteries  with  uncertainty 
about  the  level  of  X^^  only)  given  X^r.-./X^  are  fixed,  do  not 
depend  on  the  level  where  those  attributes  are  fixed. 

The  main  result  can  now  be  stated. 

Theorem  1.  For  n ^ 3,  if  for  some  X.,  (X.,X.)  is  pre- 
ferentially  independent  of  the  other  attri- 
butes for  ail  j 5^  i and  X^  is  utility 
independent  of  all  the  other  attributes, 
then  either 


or 


n 

u(x)  = r k.u. (x. ) , 

- i=i  i i 1 


n 

1 + ku(x)  = n [1  + kk.u. (x.)]  , 

“ i=l  ^ ^ 


(1) 


(2) 


where 

(i)  u and  u^  are  utility  functions  scaled  from  zero 
to  one, 

(ii)  the  k^*s  are  scaling  constants  with  0 < k^  < 1, 
and 

(iii)  k > -1  is  a non-zero  scaling  constant  satisfying 
the  equation 

n 

1 + k = n (1  + kk.) 
i=l  ^ 


(3) 
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The  proof  of  this  result  is  found  in  Keeney  [ 9 ) . Alternative 
sets  of  assumptions  leading  to  either  form  (1)  or  (2)  are 
found  in  Fishburn  [ 3 ],  Poliak  [121,  and  Meyer  [11].  The  func- 
tional form  (1)  is  referred  to  as  the  additive  utility  function 
and  (2)  is  the  multiplicative  utility  function.  For  the  case 
of  two  attributes,  the  following  is  proved  in  Keeney  [7  ]: 

Theorem  2.  For  n = 2,  if  is  utility  independent  of 
X2  and  ^2  utility  independent  of  X^,  then  the 
utility  function  u(x^,x^)  is  either  additive  or  multi- 
plicative. 

n 

Using  either  (1)  or  (2) , if  Z k.  = 1,  the  utility  function 

is  additive,  and  if  Z k.  1,  it  is  multiplicative.  When 

n " 

Z k.  > 1,  then  -1  < k < 0,  and  when  Z k.  < 1,  then 
i=l  ^ i=l  ^ 

0 < k < ®.  To  use  either  the  additive  or  multiplicative  form, 
we  need  to  obtain  exactly  the  same  information.  We  have  to 
assess  the  n single-attribute  utility  functions  u^(x^)  and  the 
n scaling  constants  k^. 


2.2  Nesting  Utility  Functions 


The  results  concerning  the  functional  forms  above  are 
valid  regardless  of  whether  the  X^’s  are  scalar  attributes  or 
vector  attributes.  This  means  that  the  x^*s  can  be  either 
scalars  or  vectors.  In  the  former  case,  the  component  utility 
function  u^  is  a uniattribute  utility  function,  whereas  in  the 
latter  case,  u^  is  itself  a multiattribute  utility  function. 
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If  is  a vector  attribute,  it  is  possible,  subject  to  satis- 
fying the  requisite  assumptions,  to  use  Theorems  1 and  2 con- 
cerning u^.  In  such  a case,  we  will  say  u^  is  a nested  MUF. 
That  is,  u^  is  a MUF  nested  within  the  MUF  u.  Our  interest  in 
nesting  utility  functions  will  become  more  apparent  in  the  dis- 
cussion concerning  the  applicability  of  the  functional  forms. 

2.3 Applicability  of  the  Functional  Forms 

In  terms  of  the  required  assessments,  the  additive  and 
multiplicative  utility  functions  appear  to  be  the  practical 
ones  for  say  n ^ 4.  Discussions  on  this  and  the  reasonableness 
of  the  assumptions  can  be  found  in  Keeney  [ 9 ] . Even  when  the 
requisite  assumptions  do  not  precisely  hold,  it  may  be  a good 
approximation  to  assume  they  do.  Furthermore,  by  nesting  one 
MUF  inside  another,  additional  flexibility  in  the  preference 
structure  can  be  achieved. 

The  effect  of  nesting  multiplicative  forms  is  to 
create  an  extra  degree  of  freedom  in  the  problem  by  having  an 
extra  independent  constant.  Without  nesting,  the  number  of 
independent  scaling  constants  is  equal  to  the  number  of  single 
attributes.  However,  suppose  u^  is  a MUF  nested  within  u and 
that  u^  has  three  single  attributes.  Then  one  would  need  n 
scaling  constants  for  the  "outer  MUF"  and  three  for  the  "inner 
MUF"  for  a total  of  n + 3,  even  though  there  are  only  n + 2 
single  attributes,  X^,...,X^^^  and  the  three  single  attributes 
in  u^.  The  degree  of  freedom  afforded  by  the  extra  parameter 
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fjermits  tradeoffs  between  two  attributes  to  be  dependent  on  a 

third.  Specifically,  tradeoffs  between  any  of  the  single 

attributes  in  u^  and  those  not  in  depend  upon  the  levels  of 

the  other  single  attributes  in  u . This  is  discussed  in  detail 

n 

in  Appendix  E. 

Jsing  various  nesting  schemes,  enough  extra  constants 
could  be  provided  to  model  situations  in  which  tradeoffs  bet- 
ween many  pairs  of  attributes  depend  on  the  level  of  other 
attributes.  That  is  to  say,  situations  in  which  the  preferen- 
tial independence  assumption  does  not  hold  for  all  the  single 
attributes  can  still  be  modeled  using  nesting. 

In  case  of  utility  independence  violations,  the  parti- 
cular problem  may  be  far  mere  sensitive  to  the  scaling  con- 
stants or  tradeoffs  among  the  attributes  than  to  the  condi- 
tional single-attribute  utility  function  variations.  Thus, 
even  in  these  cases,  the  additive  or  multiplicative  form  may 
provide  an  adequate  model  for  the  problem. 

In  summary,  the  additive  and  multiplicative  utility 
functions  are  simple  enough  to  be  tractable  and  yet,  especially 
with  nesting,  robust  enough  to  adequately  quantify  preferences 
for  many  problems.  In  practice,  however,  assessing  and  using 
such  MUF*s  is  "easier  said  than  done." 
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3.  DIFFICULTIES  WITH  EXISTING  METHODS  FOR  ASSESSMENT  AND  USE 

In  this  chapter,  existing  methods  for  assessing  and 
using  MJF*s  are  discussed.  Difficulties  encountered  with  these 
methods  include: 

(1)  the  necessity  to  ask  "extreme  value”  questions  to 

the  computational  requirements  for  specifying 
a utility  function  to  a manageaole  level, 

(2)  the  tedium  of  calculating  component  utility  func- 
tions and  scaling  constants  even  in  this  case, 

(3)  the  lack  of  immediate  feedback  to  the  decision 
maker  of  the  implication  of  his  preferences, 

(4)  the  absence  of  convenient  procedures  for  "updat- 
ing” the  decision  maker's  preferer;  s and 
conducting  sensitivity  analysis. 

In  all  that  follows,  we  will  assume  that  the  assump- 
tions implying  that  the  MUF  is  either  additive  or  multiplica- 
tive hold.  The  discussion  is  developed  in  terms  of  the  steps 
customarily  followed  in  assessing  and  using  a MUF. 

3.1 Specifying  the  Preference  Functions  over  the  Single 

Attributes 

Techniques  for  assessing  single-attribute  utility 
functions  have  become  fairly  standard  (Raiffa  [14]),  and 
sophisticated  computer  programs  have  been  developed  for 
fitting  single-attribute  utility  functions  (Schlaifer  [16]). 
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Such  programs  provide  quick  feedback  which  assists  the  decision 
maker  in  checking  if  his  assessments  and  their  implications 
appear  reasonable.  There  is  difficulty  in  using  these  programs 
for  multiattribute  utility  applications,  since  at  present,  they 
do  not  exist  in  conjunction  with  a multiattribute  utility 
assessment  package. 

3 . 2 Assessing  the  Tradeoffs  among  Attributes 

The  issue  of  tradeoffs  among  the  attributes  is  add- 
ressed by  assessing  the  k^*s  in  the  utility  functions  (1)  or 
(2) . In  theory,  the  general  method  for  doing  this  is  very 
simple.  If  there  are  n attributes,  we  want  to  assess  the  n 
unknown  k^^'s  by  creating  n independent  equations  with  the  n 
unknowns  and  solving.  An  equation  is  created  by  (i)  having 
the  decision  maker  indicate  two  options,  where  an  option  is 
either  a consequence  or  a lottery,  between  which  he  is  indif- 
ferent, and  (ii)  equating  the  expected  utility  of  these  options 
using  either  (1)  or  (2) . For  instance,  if  the  decision  maker 
finds  x'  and  x"  indifferent,  then  u(x*)  = u(x")  provides  one 
equation  with  at  most  n unknowns. 

Because  of  the  difficulty  and  tedium  in  manually  solv- 
ing n equations  (which  are  not  necessarily  linear)  with  n 
unknowns,  current  practice  in  assessing  the  k^*s  usually 
requires  sets  of  equations  which  are  simple  to  solve.  This 
basically  limits  the  assessment  questions  to  two  types.  To 
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indicate  these,  let  us  define  x*  = (x, * . ,x  *)  and 

*s  (Xj^®  ,X2®  / . . . ,x^® ) as  the  most  desirable  and  least  desir- 
able consequences.  Then,  because  of  the  scaling  conventions 
given  in  Theorems  1 and  2, 

u(x*)  = 1 , u(x®)  = 0 , (4) 

and 

Ui^Xi®)  = 1 , u^(x^®)  = 0 , i - 1,2, ...,n  . (5) 

Question  I.  For  what  probability  p are  you  indifferent  bet- 
ween 

(i)  the  lottery  giving  a p chance  at  x®  and  1-p  chance 
at  X® , and 

(ii)  the  consequence  (x®^, . . . ,x®^_^,x®^ ,x®^^^, . . . ,x® ^) . 

If  we  define  the  decision  maker’s  answer  as  p^,  then  using 
(4),  the  expected  utility  of  the  lottery  is  p^,  and  using 
either  (1)  or  (2) , the  utility  of  the  consequence  is  k^. 
Equating  the  expected  utilities,  we  find 

= Pi  (6) 

The  second  type  of  question  is  illustrated  by 

Question  II.  Select  a level  of  X^,  call  it  x^*,  and  a level 
of  Xj,  call  it  x^ * , such  that,  for  any  fixed  levels  of  all 
the  other  attributes,  you  are  indifferent  between 

(i)  a consequence  yielding  x^^*  and  Xj®  together,  and 

(ii)  a consequence  yielding  x ^ * and  x^®  together. 
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Using  (5)  and  either  the  multiplicative  or  additive  utility 
function,  the  utilities  of  these  two  indifferent  consequences 
can  be  equated  to  yield 

k.Ui(Xi’)  = k.u. (X.*)  (7) 

Once  the  single-attribute  utility  functions  u^  and  u^  are 
assessed,  both  u^(x^)  and  u^Cx^)  are  easily  found,  so  (7)  is 
a simple  linear  equation  expressing  the  relationship  between 
and  k ^ . 

A major  shortcoming  of  questions  of  both  types  is  the 
use.  of  the  extreme  levels  of  the  attributes,  that  is  the  x^*'s 
and  x^®*s.  ”.:.nce  the  range  from  to  x^*  must  cover  all  the 

possjjjie  Xj^'s,  the  implications  of,  and  h 'nee  preferences  for, 
the  extreme  levels  are  usually  very  difficult  for  a decision 
maker  to  consider.  A further  difficulty  with  Question  I is 
the  fact  that  the  effect  due  to  varying  all  n attributes  simul- 
taneously in  a lottery  must  be  considered.  Hence,  for  computa- 
tional ease,  we  must  force  the  decision  maker  to  .respond  to 
questions  much  more  difficult  to  evaluate  than  would  be 
theoretically  necessary. 

A common  practice  in  assessing  the  k^*s  would  be  to  use 
a question  I to  evaluate  the  largest  k^,  and  then  use  type  II 
questions  to  evaluate  the  magnitude  of  the  other  k^'s  relative 
to  the  largest  k^.  Once  we  have  the  k^'s,  the  additive  form 
holds  if  they  sum  to  one.  Otherwise,  the  k^'s  are  suostitutod 
into  (3)  to  evaluate  the  k for  the  multiplicative  form.  This 
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last  task  in  itself  can  be  difficult  using  only  a calculator. 

3.3  Evaluating  Alternatives  and  Sensitivity  Analysis 

Manual  calculations  are  clearly  impractical  for  evalu- 
ating alternatives.  With  uncertainty,  we  need  to  evaluate  the 
expected  value  of  u using  the  probability  distribution  describ- 
ing the  possible  consequences.  Even  with  probabilistic  inde- 
pendence among  the  X^*s.  the  computational  task  is  large.  It 
is  also  clear  that  sophisticated  sensitivity  analyses  are  out 
of  the  question  without  major  computational  help. 

On  the  other  hand,  it  is  a large  requirement  to  develop 
a special  computer  program  to  accommodate  a particular  problem. 
Such  programming  is  often  inflexible  because  of  the  special 
problem  nature  for  which  it  is  done.  For  instance,  it  would 
be  difficult  to  add  more  attributes,  to  try  different  "nesting" 
schemes,  or  explore  the  preference  structure  for  "hints"  of 
creative  new  alternatives  to  generate. 

3.4  Summary  of  Existing  Methods  and  Their  Difficulties 

Current  methods  for  assessing  and  using  MUF*s  require 
asking  very  difficult  assessment  questions,  yield  little  feed- 
back once  given  the  responses  requested  and  are  tedious  to 
implement  computationally.  These  drawbacks  can  often  result 
in  abandoning  the  decision  theoretic  approach  in  favor  of  less 
explicit  and  theoretically  well-established  but  more  expedient 
methods  for  dealing  with  specific  problems.  The  computer 
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package  to  be  described  in  the  next  chapter  is  designed  to 
remedy  some  of  r.hese  drawbacks. 
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4.  THE  COMPUTER  PACKAGE 

This  chapter  describes  the  major  features  of  a computer 
package  designed  to  alleviate  some  of  the  shortcomings  with 
existing  methods  for  assessing  and  using  multiattribute  utility 
functions.  The  package  is  referred  to  by  the  mnemonic  MUFCAP 
standing  for  "multiattribuce  utility  function  calculation  and 
assessment  package.'’  Steps  customarily  followed  in  obtaining 
and  using  a MUF  are  presented  with  a description  of  the  MUFCAP 
commands  appropriate  in  performing  the  particular  step.  Com- 
mand usage  is  illustrated  in  Chapter  5.  A concise  summary  of 
these  commands  is  in  Appendix  A and  the  program  listing  is  in 
Appendix  B. 

4.1 Commands  to  Structure  the  Utility  Function 

Structuring  a utility  function  consists  of  specifying 
a functional  form,  its  attributes,  and  the  ranges  for  each  of 
the  attributes.  MUFCAP  has  several  commands  for • structuring 
a preference  function.  The  INPUT  command  requests  a name  for 
the  utility  function  and  asks  for  the  number  of  attributes 
which  are  arguments  of  this  function.  The  package  then 
requests  a name,  and  a range  for  scalar  attributes.  The  range 
consists  of  two  numbers  which  bound  the  amounts  to  be  consi- 
dered for  each  attribute.  To  specify  a vector  attribute,  one 
inputs  a range  with  one  bound  equal  to  the  other  bound  such 
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as  0,0,  MUFCAP  recognizes  this  as  a signal  for  a vector  attri- 
bute and  notes  that  the  associated  with  that  attribute  is  a 
nested  MUF.  The  package  then  requests  the  number  of  attributes 
which  are  arguments  of  this  nested  MUF.  For  each  of  these,  a 
name  and  range  is  solicited.  Further  levels  of  nesting  could 
be  specified  if  desired  and  the  information  requested  would  be 
analogous  to  the  material  above.  After  a nested  MUF  is  com- 
pletely specified,  the  program  returns  to  ask  for  the  names 
and  ranges  for  whatever  attributes  have  not  yet  been  covered 
in  the  outer  MUF.  When  all  the  attributes  have  been -input, 
the  structure  is  complete  and  MUFCAP  requests  a new  command 
from  the  user. 

The  INPUT  command  provides  for  all  the  bookkeeping 
which  will  be  necessary  for  information  to  follow.  Each  k^ 
and  u^  (including  those  in  a nested  MUF) , can  be  assessed 
using  the  name  of  the  attribute  with  which  it  is  associated. 

The  INPUT  command  is  quite  flexible  in  having  no  logical  limit 
to  the  degree  of  nesting  allowed. 

In  addition  to  INPUT,  the  package  has  commands  for 
adding  or  deleting  attributes  to  or  from  the  utility  function. 
It  also  has  a command  to  facilitate  "regrouping”  of  the  attri- 
butes into  various  "sub-MUF*s."  In  this  way,  a model  for  a 
problem  can  be  conveniently  altered  in  terms  of  different 
nesting  schemes. 
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4 . 2 Commands  to  Specify  the  Single  Attribute  Utility 

/unctions 

The  next  step  in  assessing  a MUF  involves  specifying 
the  u^*s  for  the  single  attributes.  As  noted  in  Section  3.1, 
sophisticated  computer  programs  do  exist  for  assessing  single 
(scalar)  attribute  utility  functions.  One  could  incorporate 
these  into  MUFCAP.  However,  for  simplicity,  several  less 
sophisticated  routines  for  assessing  unidimensional  utility 
functions  (referred  to  as  UNIF's)  were  developed. 

MUFCAP  has  a command  UNISET  for  ecifying  any  of 
three  UNIF  types;  linear,  exponential,  and  piecewise  linear. 
Pratt  [13]  considers  the  implications  of  these  forms.  The 
linear  utility  function  implies  risk  neutrality.  This  form 
requires  no  more  information  than  the  range  of  the  attribute. 
The  exponential  form  implies  constant  risk  aversion  or  con- 
stant risk  proneness.  It  requires  the  specification  of  a 
certainty  equivalent  for  a single  lottery.  Given  this,  the 
exponential  form  is  fitted  and  scaled  automatically  by  the 
program.  The  piecewise  linear  utility  function  is  specified 
by  providing  the  abscissa  and  ordinate  values  for  n points 
(3  £ n £ 15)  of  the  utility  function.  This  form  can  be  used 
for  non-monotonic  or  S-shaped  utility  functions.  These  three 
types  provide  the  user  with  the  means  of  specifying  a UNIF 
appropriate  for  many  situations.  More  forms  can  easily  be 
added  to  the  package  in  the  future. 
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MUPCAP  also  has  commano.  which  enable  a user  to  quickly 
display  the  assessed  UNIF  for  purposes  of  checking  its  appro- 
priateness. The  conunand  UNICAL  calculates  the  utility  for  one 
or  a series  of  attribute  levels.  INVERSE  calculates  the  attri- 
bute level  corresponding  to  a given  utility  value.  LOTTERY 
evaluates  the  certainty  equivalent  for  any  lottery  with  n con- 
sequences and  their  associated  probabilities  over  that  attri- 
bute,  where  2 £ n < 15. 

To  summarize,  MUFCAP  has  commands  to  conveniently  set 
those  u. *s  which  are  UNIF*s  and  to  display  them  for  feedback 
purposes  to  check  on  their  reasonableness. 


i-1-2 to  Specify  the  Scalincr  Constants 

Using  the  attribute  names  as  identifiers,  MUFCAP 
allows  the  user  to  set  the  scaling  constants  in  the  MUF  corres 
ponding  to  each  attribute.  If  x.  is  a vector  attribute,  the 
associated  with  it  is  a MUF  with  its  own  internal  scaling 
constants.  By  referring  to  the  name  of  this  vector  attribute, 
the  user  can  specify  the  internal  scaling  constants  for  the 
associated  nested  MUF.  when  all  the  k.’s  for  a particular 

MUF  have  been  set,  the  program  automatically  calculates  the 
corresponding  k using  (3). 

Once  the  u^’s  have  been  evaluated,  the  package  has 
several  commands  useful  for  assessing  the  k. -s  in  any  parti- 
cular MUF.  The  command  INDIF2  takes  as  input  two  pairs  of 
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two  indifference  consequences  each.  These  consequences  can 
vary  only  in  terms  of  the  two  attributes  whose  k^'s  are  the 
object  of  assessment.  Then,  using  (2) , the  program  computes 
the  relc*.ive  k^*s  (i.e.,  the  ratio  for  attributes  i and 

j)  implied  by  the  indifference  pairs.  With  INDIF2,  the  user 
is  not  limited  to  choosing  consequences  which  have  one  attri- 
bute at  a least  desirable  level  in  order  to  determine  the 
relative  k^*s. 

Once  v;e  know  the  relative  k^'s,  we  can  assign  k^*s  in 
(2)  by  arbitrarily  setting  one  k^  to  a fixed  value  and  the 
others  in  terms  of  the  fixed  k^.  The  command  INDIFl  can  then 
be  used.  It  takes  as  input  a single  pair  of  indifference  con- 
sequences and  computes  the  k,  and  the  magnitude  of  the  k^^'s 
implied  by  that  pair  and  the  currently  assigned  k^*s.  It  does 
this  by  computing  the  factor  by  which  the  currently  assigned 
kj^'s  need  to  be  multiplied  to  be  consistent  with  the  indif- 
ference pair  just  given.  MUFCAP  provides  a routine  which 
allows  the  user  to  multiply  the  currently  assigned  k^'s  for 
any  MUF  by  any  factor  thus  resetting  them.  In  this  way, 
iwDIFl  enables  the  calculation  of  the  magnitude  of  the  k^'s 
using  an  indifference  relation  instead  of  a lottery  over  all 
the  attributes  at  once.  For  consistency  checks,  a new  indif- 
ference pair  of  consequences  can  be  input  using  INDIFl,  which 
then  computes  the  factor  described  above.  If  this  factor  is 
close  to  1,  the  indifference  pair  is  consistent  with  the  cur- 
rently assigned  scaling  factors. 
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Once  the  k^*s  for  a MUF  have  been  assigned,  an  indif- 
ference curve  (see  Appendix  E)  over  any  two  attributes  in  that 
MUF  can  be  calculated  with  the  command  IMAP.  IMAP  permits  a 
user  to  get  immediate  feedback  on  the  tradeoff  implications  of 
the  k^*s  or  indifference  pairs  which  he  has  specified.  He  can 
quickly  see  if  the  points  "claimed"  to  be  indifferent  really 
appear  so  to  him.  If  not,  the  k^*s  can  be  changed  or  other 
indifference  pairs  solicited  until  they  represent  more  accu- 
rately the  user's  preferences  for  tradeoffs  between  those 
attributes.  If  desired,  IMAP  can  be  used  in  conjunction  with 
INDIF2  and  other  commands  to  produce  indifference  curves  over 
two  attributes  before  all  the  other  k^^'s  have  been  assessed. 
This  is  discussed  in  Chapter  6 and  Appendix  E. 


4 . 4 Commands  for  Evaluating  Alternatives  and  Sensitivity 

Analysis 

Once  the  u^'s  and  k^*s  have  been  set,  the  utility 
function  is  completely  specified  and  can  be  used  to  evaluate 
alternatives-  MUFCAP  has  commands  for  specifying  two  kinds 
of  alternatives;  certain  and  uncertain.  For  certain  alter- 
natives, which  are  simply  consequences,  uniattribute  amounts 
are  solicited  until  the  alternative  is  completely  described. 
For  uncertain  alternatives,  at  present,  MUFCAP  assumes  proba- 
bilistic independence  and  requests  a probability  distribution 
function  for  each  scalar  attribute.  The  probability  distri- 
bution function  currently  used  is  a piecewise  linear 
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approximation  to  the  cumulative  probability  distribution  for 
X^.  The  user  supplies  n abscissa-ordinate  pairs,  where 
2 £ n £ 9 to  specify  the  cumulative  distribution.  The  cumula- 
tive distribution  was  chosen  rather  than  the  probability  den- 
sity function  because  the  fractile  method  of  assessing  probabi- 
lities (see  Schlaifer  (151)  yields  points  of  the  cumulative 
distribution  Other  forms  of  probability  distributions  such 
as  the  Gaussian  as  well  as  probabilistic  dependencies  could 
be  added  to  the  package  in  the  future. 

The  specified  alternatives  are  given  names  by  the 
user.  With  these  names,  the  user  may  add,  change  or  delete 
alternatives.  He  may  also  choose  the  ones  which  are  to  be 
evaluated  by  listing  their  names  with  the  appropriate  commands 
about  to  be  described. 

The  command  EVAL  is  used  to  evaluate  (i.e.,  compute  the 
expected  utility)  for  any  alternative  or  group  of  alternatives. 
EVAL  can  compute  the  expected  utility  for  the  overall  utility 
function  or  for  the  utility  function  associated  with  any  par- 
ticular attribute.  In  the  latter  case,  attribute  levels  in  an 
alternative  which  are  not  arguments  of  the  particular  utility 
function  are  ignored.  Typically,  EVAL  can  be  used  to  evaluate 
alternatives  for  the  current  multiattribute  model.  Parameters 
such  as  the  scaling  constants  or  probability  distributions  can 
then  be  changed  and  the  alternativer.  evaluated  again.  In  this 
way,  we  can  see  how  sensitive  the  rankings  are  to  changes  in 
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certain  parameters.  In  a group  decision-making  context,  dif- 
ferent utility  functions  and  probability  estimates  of  group 
members  can  be  used  to  evaluate  and  rank  the  alternatives. 


This  might  help  clarify  differences  of  opinion  and  suggest 
certain  creative  compromises  or  areas  where  more  precise  pro- 
bability estimates  may  be  needed. 

The  command  GRAD  evaluates  the  gradient  of  a utility 
function  at  any  number  of  specified  alternatives.  The  gradient 
is  defined  as  the  vector  r ^ ‘ ‘ ^ indicates 

the  direction  of  steepest  increase  in  the  utility  function  at 
a specified  point.  The  gradient  component  tells  us  which 
attribute  level  changes  would  yield  large  increases  in  utility. 
This  could  be  useful  in  generating  improved  alternatives  to 
the. current  one.  Of  course,  one  must  keep  in  mind  the  scales 
of  the  attributes  in  interpreting  the  gradient. 


In  addition  to  the  gradient,  GRAD  also  ccmoutes  the 

vector  ("1^  , * ’ * ' — I ’ Each  component  represents  the 

\ ^^2  ® n/ 

rate  of  change  of  u with  respect  to  a change  in  the  utility  u^. 
These  components  reveal  the  attributes  for  which  an  increase 


in  its  utility  will  yield  the  largest  increase  in  u.  The 
advantage  of  calculating  these  quantities  in  addition  to  the 


gradient  components  are  (a)  components  can  be  calculated  for 
MUF*s  as  well  as  UNIF’s,  and  (b)  the  unit  of  measurement  for  a 


uniattribute  does  not  distort  the  magnitude  of  the  component. 

Thus  in  some  cases,  • - might  give  a better  picture  of 

cu . 
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possible  improved  alternatives  than  . MUFCAP  makes  both 

V " 

available. 

Summarizing,  EVAL  permits  the  evaluation  of  alterna- 
tives, and  along  with  routines  which  alter  parameters,  provides 
for  sensitivity  analysis.  GRAD  makes  use  of  the  analytical 
formulation  of  the  problem  to  calculate  quantities  useful  in 
suggesting  improved  alternatives  to  the  currently  specified 
ones. 


4.5 General  Command  Format  and  Commands  for  Facilitating 

Use  of  the  Package 

MUFCAP  commands  are  designed  to  be  concise  and  are  for 
the  most  part  no  longer  than  three  words.  These  words  may 
initiate  a dialogue  when  more  information  is  necessary.  The 
input  format  is  free,  i.e.,  words  need  not  begin  in  a particu- 
lar position  on  the  page.  For  many  commands,  the  user  will  be 
prompted  if  he  has  left  out  a necessary  word. 

Mistyping  causing  invalid  numbers  on  input  is  handled 
automatically  by  the  program  and  a correct  number  is  requested. 
Provision  is  made  for  the  user  to  terminate  a lengthy  dialogue 
by  specifying  the  word  QUIT  for  the  next  number  to  be  input. 

A new  command  can  then  be  entered.  In  the  future,  a help 
command  could  be  easily  implemented  whir  i would  explain  the 
syntax  of  any  other  command,  give  definitions  of  terms  used 
in  the  program  and  make  suggestions  concerning  what  kinds  of 
steps  to  perform  in  assessing  and  using  the  MUF. 
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In  addition  to  these  features,  MUFCAP  has  the  facility 
for  saving  the  current  status  of  the  multiattribute  utility 
structure  and  the  current  alternatives  in  a file  of  the  user's 
choosing  to  be  read  in  at  a later  time.  This  gives  MUFCAP  the 
capability  for  filing  away  several  different  MUF  models  as  well 
as  a large  number  of  alternatives  for  the  same  problem.  It 
also  allows  the  user  to  build  up  his  model  over  many  different 
sessions  rt  the  terminal  and  restore  any  status  he  has  saved 
away  with  which  he  wishes  to  calculate  at  any  particular  time. 

Another  feature  of  MUFCAP  is  the  supplying  of  default 
settings  when  the  INPUT  command  is  used  to  structure  the  MUF 
for  the  problem.  After  INPUT,  the  default  for  all  MUF's  is 
the  additive  form,  with  all  the  k^'s  equal  to  each  other,  and 
for  all  UNIF's,  it  is  the  linear  utility  function.  With  these 
defaults,  the  user  is  set  to  calculate  immediately  after  input. 
Thus  feedback  can  begin  right  away  without  requiring  the  user 
to  completely  specify  everything  first.  Scaling  constants  and 
utility  functions  can  then  be  altered  after  observing  some 
feedback  to  refine  the  model  for  the  problem. 

Finally,  MUFCAP  provides  commands  to  print  out  the 
current  status  of  the  assessments.  There  are  routines  to  dis- 
play the  k^'s  and  k for  any  MUF,  the  range  and  type  for  any 
scalar  attribute  utility  function,  the  probability  distribu- 
tion of  any  attribute  for  any  alternative,  the  multiattribute 
utility  function  structure  (i.e.,  nesting)  and  the  currently 
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defined  alternatives.  Commands  are  also  provided  for  easily 

changing  parameters  such  as  individual  k.'s  or  the  components 
of  any  alternative. 
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5.  APPLICATION  OF  THE  PROGRAM  TO  DIFFERENT  PROBLEMS 

This  chapter  presents  several  applications  designed  to 
show  how  MUFCAP  can  be  used  in  practice.  Certain  application 
descriptions  contain  computer  printout  illustrating  the  use  of 
various  MUFCAP  commands.  Each  set  of  computer  printout  is 
followed  by  a comments  section  which  summarizes  the  pertinent 
leatures  illustrated  by  the  printout.  Reference  to  Appendix  A 
when  reading  the  printout  and  comments  is  recommended. 

Application  of  MUFCAP:  The  Mexico  City 

The  Mexico  City  Airport  problem  concerned  the  decision 
for  developing  the  city's  airport  facilities  in  the  most 
"effective"  manner  in  a multiobjective  sense.  The  analysis 
which  was  done  is  described  in  more  detail  in  Keeney  [ 8 ] . 

This  problem  was  approached  using  the  existing  methods  for  MUF 
assessment  and  utilized  special  computer  programming  to  aid  in 
the  calculations.  This  section  presents  what  might  have  been 
done  if  MUFCAP  had  been  available  then. 

5.1.1  Attributes  for  the  Problem 

The  Mexico  City  Airport  problem  was  defined  in  terms 
of  the  following  attributes: 

Xi  = total  cost  in  millions  of  pesos 

Xj  = the  capacity  in  terms  of  the  number  of  aircraft 
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operations  per  hour 

= access  time  to  and  from  the  airport  in  minutes 
= number  of  people  seriously  injured  or  killed  per 
aircraft  accident 

Xg  = number  of  people  displaced  by  airport  develop- 
ment 

Xg  = number  of  people  subject  to  a high  noise  level; 
(i.e.,  90  CNR  or  more) 

To  incorporate  time  effects  of  building  the  airport,  attri- 
butes were  defined  using  present  values  or  averages  where 
appropriate.  The  capacity  attribute  X2  had  to  be  made  a func- 
tion of  capacity  for  1975,  capacity  for  1985,  and  capacity  for 
1995,  and  thus  it  was  a vector  attribute. 

5.1.2  Summary  of  the  Method  Used  in  the  Problem 

After  verifying  assumptions  concerning  preferential 
and  utility  independence  and  ascertaining  the  appropriateness 
of  the  multiplicative  model,  assessments  were  begun.  First, 
the  fractile  method  was  used  to  obtain  probability  distribu- 
tions for  all  of  the  alternatives  under  consideration.  Pro- 
babilistic independence  was  assumed  to  simplify  calculations. 
Then  uniattribute  utility  functions  were  assessed  for  ail 
eight  scalar  attributes.  The  k^*s  were  assessed  using  the 
lottery  over  all  the  attributes  illustrated  by  Question  I in 
Section  3.2  for  bot^'  th?.  overall  MUF  and  nested  capacity  MUF. 
Consistency  checks  on  the  relative  involving  tradeoffs 
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of  two  attributes  at  a time  (see  Question  II,  Section  3.2) 
were  also  employed.  Special  computer  programs  and  graphic 
displays  were  developed  for  evaluating  alternatives  and  sen- 
sitivity analysis.  For  sensitivity  analysis,  the  program 
allowed  changes  in  (a)  the  endpoints  for  the  fractile  cumula- 
tive pro!  ability  distributions  and  (b)  in  the  scaling  factors 
k^.  The  shapes  of  the  utility  functions  or  the  cumulative 
probability  distributions  could  not  be  changed  without  pro- 
gramming adjustments. 

5.1.3  A MUFCAP  Approach  to  the  Mexico  City  Problem 

The  MUFCAP  approach  would  follow  the  existing  methods 
scheme  in  making  and  verifying  the  preferential  independence 
and  utility  independence  assumptions.  The  INPUT  command  would 
structure  the  multiplicative  function  giving  names  such  as 
"cost”  and  "access"  to  the  various  attributes  along  with 
ranges  for  the  attribute  amounts.  Capacity  would  be  put  in 
as  a nested  MUF. 

Alternatives  would  be  specified  by  inputting  the  nine- 
point  assessed  fractile  distribution  for  each  uniattribute  of 
a particular  alternative.  Utility  functions  for  single  attri- 
butes would  be  specified  using  any  of  the  three  forms  avail- 
able in  MUFCAP. 

Assessment  of  the  k^*s  could  be  accomplished  without 
depending  upon  the  supplying  of  the  probability  for  a lottery 
over  all  the  attributes  as  was  done.  Pairs  of  indifference 
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points  for  two  attributes  would  be  fed  into  MUFCAP  to  imme- 
diately produce  indifference  curves  for  examination  and  veri- 
fication by  the  decision  maker.  In  |:his  way,  the  relative 
k^*s  would  be  established  with  the  aid  of  feedback.  The  mag- 
nitude of  the  k^*s  could  be  established  using  INDIFl  (see 
Section  4.3),  so  a lottery  over  all  the  attributes  could  be 
avoided  for  this  purpose.  A good  consistency  check  would  be 
provided  by  comparing  the  magnitude  of  the  k^*s  implied  by 
each  method.  Using  MUFCAP,  all  of  the  initial  assessments 
could  be  made  and  stored  for  later  use.  The  assessments  would 
have  been  made  with  the  aid  of  immediate  feedback  and  with  no 
need  for  very  difficult  lottery  questions. 

After  the  initial  assessments,  alternative  evaluations 
and.  sensitivity  analysis  could  be  performed  immediately  with 
no  need  for  special  programming.  Fractile  distributions  and 
utility  function  shapes  could  also  be  altered  without  program- 
ming adjustments.  The  different  assessments  of  various 
individuals  and  groups  could  have  been  filed  away  for  later 
reference  using  MUFCAP 's  filing  capability. 

In  addition,  other  possibilities  could  have  been 
explored  with  a minimum  of  extra  effort.  New  attributes  such 
as  air  pollution  and  political  effects  could  be  added  into  the 
analysis  with  no  special  programming.  The  gradient  calcula- 
tion capability  may  have  been  used  to  support  other  alterna- 
tives for  exploration  and  development.  If  the  preferential 
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independence  of  some  attributes  are  questioned,  different 
nesting  schemes  could  be  tried  to  see  if  the  ranking  of  the 
alternatives  would  be  affected.  Thus  MUFCAP  could  have  pro- 
vided the  analysis  that  was  performed  with  no  special  program 
ming  and  might  have  been  used  to  explore  variations  of  more 
parameters,  other  multiattribute  nesting  schemes,  and  addi- 
tions  of  new  attributes. 
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5.1.4  Mexico  City  Airport  Illustrations 

lo^on  a Ian  size (300)  nono 
entf:p  p\5swoRn  for  mmi- 

M20225  .119/*0  ^Cf:0U^IT  \RE  LOW.  SEE  USER  \-rouMTS. 

ALAM  LOnO’l  l!l  PROGRESS  AT  10:33:«*0  ON  APRIL  29^  1975 
NO  BROADOAST  MESSAGES 
READY 

allocate  file (mexico)  dataset (mex Ico) 

READY 

call  nufcap 

TEMPNAME  ASSUMED  AS  A MEMBER  .‘lAME 

COMMAiID  WORD  AND  FILE  NAMES  MUST  BE  IN  CAPS  . 

COMMAND?  : 


Illustration  1 


READ  MEXICO 


COMMAND?  : PEBUn 

STRUCTURE  FOR  m«xlco 


cost 

0.4  80 

4.onnooE+'^3 

5.00000E+02 

capacity 

O.GOO 

cap75 

0.300 

S.OOOOOE^Ol 

1.30000E+02 

cap85 

0.500 

s.onoooE^oi 

2.00000E^02 

cap95 

0.400 

1.00300E+02 

2.50000E^02 

access 

0.100 

9.00000E-»-01 

1.20000E^01 

safety 

0.550 

1. 00000 E^ 03 

l.OOOOOE^OO 

dl spl acement 

0.130 

2 .500aOF>95 

2.50000E^03 

noise 

0.180 

1.50000E^03 

2.00000E+00 

COMMAND?  : 


Illustration  2 
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n I SPLAY  rnexlco 


L!STmn  OF  K 
cost 

capac  It y 
access 
safety 
HI spl aconent 
no  I se 


•ACTORS 

O.fiOn 
0.100 
0.?I50 
0.180 
0.180 


HMiK=  -0..T77  stin  K'r;  = i 890 
COMfrOMP?  . display  capncity 

LISTI.’IG  OF  K F\CTORS 
cap7S  0.300 

capS5  0.500 

cap 55  0.40n 

bigk=  -0.1*53  sun  k's  = 

COWnriD?  : DISPLAY  access 


1.200 


RAIIGE:  90.000  12.000 

''(X)=n(l-EXP(-CX)) 
RIGK  AVErL  VARIABLE  fJORHALIZEn 

COMMANP?  : 


Illustration  3 


UN ISFT  access  CR 

input  ANY  50-50  LOTTERY  |f|  THE  FORM  OF  C.E.,Q1  A 02. 
02  12  00 


PLEASE 


COMMAND?  : 

Ull  IC  AL 

access 

U( 

90. 

000)  = 

0 

.000 

U( 

74. 

400)  = 

0 

.304 

U( 

58. 

300)  = 

0 

.54  4 

U( 

43. 

200)  = 

0 

.733 

U( 

27. 

r>oo)  = 

0 

U( 

12. 

000)  = 

1 

.000 

COMMAND?  : 

INVERSE 

access  2 

.75 

77.463=IiIV( 

0* 

*250) 

41.61 

7=i;iv( 

0, 

.750) 

COMMAND?  : 

Illustration  4 


-43 


LOTTERY  access  3 
LOTTERY  EIIDPTS.  PLEASE? 

20  40  60 

CORRE3P.  PROB^BIL ITIES  PLEASE? 

.3  .4  .3 

CE  FOR  LOTVERY=  41.816 
COMMAND?  ; 

Illustration  5 


cost 

cap75 

cap85 

cap95 

access 

safety 

di splaceme 

noise 


ALTLIST 
al 1 one 

500.000 

130.000 

200.000 

250.000 

12.000 

1 .000 

2500.000 

2.000 

TABLE  FOR 


ftllhalf 

2250.000 

90.000 

140.000 

175.000 

51.000 
500.500 

126250.000 

751.000 

PROB  ALTErci 


500.000 

130.000 

200.000 

250.000 

12.000 

1000. 000 

250000.000 

1500.000 


CERT  EQUIV 
NO  PROB.  ALTERM 
COMMAND?  : 


Illustration  6 


al 1 one 
a 1 1 ha  1 f 
a3 

COMMAND?  : 
a 11  ha  1 f 

COMMAND?  : 
all  one 
a 1 1 ha  1 f 
a3 

COMMAND?  : 
al lone 
allhal f 
a3 

COMMMJD?  : 


EVAL  mcxlco 

1.000 
0.841 
0.855 

EVAL  mex Ico  al  1 hal  f 
0.841 

EVAL  capacity 
0.993 
0.805 
0.099 

EVAL  access 

1.000 

0.044 

1.000 


Illustration  7 
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KSei  mex \co  ADD 

PIGK«  0,000 
COMMANO?  : Pi  SPLAY  mexico 


listing  0*^  K FACTORS 


cost  0.254 

capacity  0.^17 

access  0.053 

s^^oty  0.1f55 

^displacement  0.095 
noise  O.095 

BIGK=  0.000  Sl.»f1  K*S 
COMflANP?  : EVAL  nexlco 
all  one  1.000 

allhalf  0.679 

0.624 

COMMAND?  : 


1.000 


Illustration  8 


READ  MEXICO 


COMMAND?  : APPALT  all -fourth  .25 

ALTERNATIVE  all-fourth  SPECIF. 
COMMAfJP?  : EVAL  mex  Ico  al  1 -fourth 
all-fourth  0.C16 

C0MM\NP?  : hropalt  all -fourth 

COMMAiNP?  : 


Illustration  9 


INPIFI  safety  cost 


INPUT  AN  INPIFFEREflCE  PAIR  PLEASE 
800  1000  300  2500 


IMPLIEP  tIE'./  K'S  FACTOR(S) 
iMPLiFn  riEi;  .oic;i(=  -n.<?59 
COMMAND?  ; 


0.970  { 


4.750) 


Illustration  10 
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INHIF2  sa^r-ty  cost 

INPUT  2 INDIFFFRFNCE  PAIRS  PLEASE 

• 

800  1000  300  2500 
♦:200  3500  750  2500 

BIGK-  -0. 2 07/ K( safety  ) 

INniF  PAIR  YIELDS  IIIFO  AROUT  REL  K'S 

REL  K CHECK.  CURREfIT  RATIO  cost  TO  safety-  1.571 

IMPLIED  RATIO  = 1.397 

COMMAND?  : 


Illustration  11 


IMAP  safety  cost 

INPUT  IIIDIF  PT.  THROUGH  WHICH  CURVE  WILL  PASS:  500  2500 


INPUT  HUMBER  OF  PTS . FOR  HAP:  5 


INPUT  safety  VALUES  FOR  HAP 

300  400  500  600  700 


INHIFFEREUCE  PTS 

( 

300.000, 

2922.530) 

( 

400.000, 

2715.855) 

( 

500.000, 

2500.002) 

( 

600.000, 

2272.636) 

( 

700.000, 

2050.779) 

UTIL 

FOR  CURVE 

V/ITH  OTHER  ATTR 

COMM.ANP?  : 


0.444 


Illustration  12 


INTER3K  mexico 

BIGK=  -0.453  INTERBK=  -0.526 
Illustration  13 


capacity 

coniiANn? 
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5.1.5  Comments  on  Mexico  City  Airport  Illustrations 
Illustration  1 

The  user  logs  in,  sets  up  a data  file  which  will  be 
used  and  invokes  MUFCAP. 

Illustration  2 

The  status  of  preferences  and  alternative  specifica- 
tions in  the  file  MEXICO  is  read  in.  The  multiattribute 
utility  function  structure  is  displayed. 

Illustration  3 

Characteristics  of  MUF's  and  UNIF's  associated  with 
various  attribute  names  are  displayed.  Mexico  and  capacity 
have  associdced  MUF’s  while  access  has  an  associated  UNIF. 

Illustration  4 

An  example. of  setting  a UNIF  is  shown.  The  UNIF  for 
access  is  assumed  to  be  of  the  constant  risk  type.  The  UNIF 
is  fitted  in  response  to  the  50-50  lottery  certainty  equiva- 
lent request.  UNICAL  tabulates  the  UNIF  for  various  amounts 
of  access.  INVEPSE  tabulates  the  amounts  of  access  having 
certain  utility  values.  The  amount  of  access  having  utility 
=.25  should  correspond  to  the  certainty  equivalent  for  the 
50-50  lottery  between  the  amount  of  access  having  utility  =.5 
and  that  having  utility  =0.  A check  with  Keeney  [ 8]  shows 
that  the  fit  for  access  appears  to  be  very  good. 
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Illustration  5 

An  example  using  the  LOTTERY  command  is  shown.  A cer- 
tainty equivalent  for  the  3-conse(iuence  lottery  is  output. 

Illustration  6 

Several  "certain”  alternatives  are  displayed. 

**allone"  has  all  the  attributes  at  their  best  levels.  "^3" 
has  cost,  capacity  and  access  at  their  best,  and  safety,  dis- 
placement and  noise  at  their  worst,  "allhalf"  has  all  the 
attributes  halfway  between  their  range  limits.  There  are  no 
uncertain  alternatives  in  this  current  status. 

Illustration  7 

This  illustrates  the  use  of  the  EVAL  command.  The 
overall  utility  function  mexico  is  evaluated  for  all  the 
alternatives;  and  then  only  for  allhalf.  The  MUF  associated 
with  capacity  is  evaluated  for  all  the  alternatives.  The  UNIF 
associated  with  access  is  similarly  evaluated. 

Illustration  8 

These  lines  illustrate  a little  sensitivity  analysis. 
The  K5ET  command  makes  the  overall  utility  function  "mexico" 
additive  but  maintains  the  same  relative  k^^*s.  The  alterna- 
tives are  then  evaluated.  Notice  the  change  in  rank  between 
"allhalf"  and  "a^"  with  the  additive  model  as  opposed  to  the 


original  model. 
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Illustration  9 

The  original  model  is  restored.  An  alternative  all- 
fourth is  added,  evaluated  and  dropped. 

Illustration  10 

A check  on  the  magnitude  of  the  k^'s  is  performed 
using  INDIFl  and  a single  ndif ference  pair.  The  check  shows 
that  the  current  k ' s agree  well  with  the  indifference-pair 
check. 

Illustration  11 

An  independent  check  is  made  on  the  relative  k^*s  con- 
cerning "cost"  and  "safety."  The  implied  ratio  agrees  well 
with  the  current  ratio. 

Illustration  12 

An  indifference  curve  is  tabulated  between  "cost"  and 
"safety. " 

Illustration  13 

A check  is  made  on  the  necessity  for  nesting  capacity 
as  opposed  to  using  the  attributes  cap75,  cap85  and  cap95 
along  with  the  others  in  a single  d-attribute  multiplicative 
form.  The  check  shows  that  without  nesting  the  apnroximation 


to  the  tradeoffs  among  the  attributes  would  be  pretty  good. 
(See  Appendix  E for  a more  detailed  explanation.) 
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5*2 Evaluation  of  a Computer  Time-Sharing  System 

This  section  concerns  an  example  relevant  to  a manager 
of  a time-sharing  system  in  formulating  a MUF  to  evaluate  dif- 
ferent courses  of  action.  The  data  and  formulation  is  based 
on  Grochow  [ 4 ] . This  problem  was  also  approached  using 
existing  methods  and  special  computer  programming.  A possible 
MUFCAP  approach  is  presented  here. 

5.2.1  Attributes  for  the  Problem 

The  following  attributes  were  used  in  the  time-sharing 

problem: 

A = Availability  measured  in  percentage  of  successful 
logins 

RT  = Average  response  time  to  majority  of  trivial 
requests  in  seconds 

RC  = Average  response  time  to  majority  of  compute- 
bound  requests 

5.2.2  Summary  of  the  Method  Used  in  the’  Problem 
The  first  stage  of  analysis  was  to  determine  what 

utility  independence  relationships  existed  among  the  attri- 
butes. It  was  found  that  RC  was  utility  independent  of  A,  and 
RT  was  utility  independent  of  A and  RC.  But  A was  not  utility 
independent  of  RT  or  RC,  and  RC  was  not  utility  independent  of 
RT.  Examination  of  the  attributes  showed  that  certain  forms 
of  independence  were  not  to  be  expected.  For  example, 
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tradeoffs  between  RC  and  A may  depend  on  RT  since  it  hardly 
pays  to  be  able  to  log  in  more  often  if  RT  is  very  bad. 

Grochow*s  approach  was  to  formulate  an  overall  utility 
function  involving  seven  conditional  one-attribute  utility 
functions  and  effectively  assessing  six  scaling  constants 
using  existing  methods. 

5.2.3  A MUFCAP  Approach 

A possible  MUFCAP  approach  to  this  problem  would  be 
to  try,  as  an  approximation,  the  following  nesting  scheme: 

u(a,rt,rc)  - u(u^,u^) 

o r 

where  u = u^(rt,rc)  and  u = u (a) 

This  is  the  multiplicative  form  with  u^  as  a nested 
MUF.  There  are  four  independent  scaling  constants  possible 
in  this  formulation.  The  model  is  assuming  as  an  approxima- 
tion that  the  various  violations  of  utility  independence  can 
be  ignored  but  that  preferences  for  tradeoffs  between  avail- 
ability and  any  response  time  depend  on  the  level  of  the  other 
response  time.  This  seems  reasonable  since  tradeoffs  between 
response  times  are  of  concern  after  the  user  has  logged  in. 

On  the  other  hand,  the  value  o^  logging  in  (e.g.,  the  amount 
one  is  willing  to  trade  to  ga  i a faster  RC)  may  depend  on  how 
good  RT  is. 

To  test  out  this  MUFCAP  approach,  we  can  calibrate 
the  MUFCAP  model  using  the  graphical  data  in  Grochow  ( 4 1 . 

This  data  provides  enough  information  to  attempt  setting  of 


-51- 


the  scaling  constants  for  the  MUFCAP  model.  In  calibrating 
the  scalar  attribute  utility  functions#  an  "average"  constant 
risk  form  for  each  attribute  was  estimated  from  the  data. 

After  calibrating  the  model,  various  points  in  the 
attribute  space  (i.e.,  alternatives)  were  evaluated  and  ranked 
to  see  how  closely  they  compared  to  the  graphical  data  in 
j Grochow  [ 4 ] . The  results  illustrate  ' in  the  computer  print- 

outs following  this  section  were  reasonably  close  to  the  gra- 
phical data  and  seemed  to  justify  the  MUFCAP  approximation 
scheme.  The  agreement  seemed  reasonable  in  spite  of  the  fact 
that,  constant  risk  forms  were  used  for  the  scalar  attribute 
' utility  functions.  The  graphical  data  exhibited  "jumps" 

^*'^^3ch  could  be  modeled  by  piecewise  linear  forms  in  ci  more 
refined  approximation. 

If  one  is  satisfied  with  the  MUFCAP  approximation,  we 
I can  immediately  proceed  to  perform  gradient  calculations  show- 

ing which  direction  one  should  take  for  maximum  improvement 
of  the  current  state  (in  the  attribute  space)  as  Grochow 
suggests.  Also,  expanding  the  model  to  include  more  attri- 
butes (e.g.,  cost)  seems  easier  with  the  MUFCAP  scher a than 
with  further  conditional  utility  functions  and  "corner  point" 
(i.e.,  extreme  value)  assessments  for  scaling  constants. 

To  summarize,  MUFCAP,  with  nesting,  may  be  used  to 
capture  the  essential  features  of  situations  which  may  not 
satisfy  some  of  the  independence  assumptions.  When  the 
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apprcximation  can  be  used,  gradient  calculations,  sensitivity 
analysis  and  expansion  of  the  model  to  include  more  attributes 
become  feasible  using  MUFCAP. 


5.2.4  Computer  Time-Sharing  System  Illustrations 


INPUT  grochow 

HOW  lUNY  ATTRIBUTES  ARE  IN  THIS  MUF?  : 2 
INPUT  IIAl'iE  AND  RANGE  FOR  ATTR  1 OF  UTIL  FUNC  Rrochow 
a .1  1 

INPUT  NAME  AND  RANGE  FOR  ATTR  2 OF  UTIL  FUNC  grochow 
response  0 0 

HOW  MANY  ATTR.  ARE  IN  THIS  MUF?:  2 
INPUT  NAME  AND  RANGE  FOR  ATTR  1 OF  UTIL  FUNC  response 
rt  9 2 

INPUT  NAME  AMD  RANGE  coR  ATTR  2 OF  UTIL  FUNC  response 
rc  120  2 

COMMAMP?  : HE  BUG 


STRUCTURE  FOR  zrochovi 


a 

0.500 

9.99999E-02 

l.OOOOOE+OO 

response 

0.500 

rt 

0.500 

9,00000F.+  00 

2.00n00E+00 

rc 

0.500 

1.20000E^fl2 

2.00003E+00 

COMMAND?  : 


Illustration  14 
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UtnSET  a CR 

tNPUT  ANY  50-50  LOTTERY  Itl  THE  FORM  OF  C.E.,01  A 02. 

.1  1 

COMMANO?  : urn  SET  rt  CR 

INPUT  ANY  50-50  LOTTERY  IN  THE  CORf<  OF  C.E.,Q1  A Q2 . 
5 9 2 

COMMAND?  : UN  I SET  rc  CR 

INPUT  ANY  50-50  LOTTERY  IN  THE  FORM  OF  C.E.,01  A 02. 
20  120  2 

COMMAND?  : 

Illustration  15 


INDIFl  rt  rc 

INPUT  AN  INDIFFERENCE  PAIR  PLEASE 
5 120  9 2 

INDIF  PAIR  YIELDS  INFO  ADOUT  REL  K'S 
REL  K CHECK.  CURRENT  RATIO  rc 
IMPLIED  RATIO  = O.SOO 
COMMAND?  ; KSET  rosponsf* 

•■t  » :.G67 

rc  - ;.553 


TO  rt  = 1.000 


BIOK=  0.000 
COMMAND?  ; 


Illustration  16 


INDIFl  rt  rc 

INPUT  AN  INDIFFERENCE  PAIR  PLEASE 
5 2 2 120 

IMPLIED  NEW  K'S  FACTOR(S)  1.000  (1254  O05) 
IMPLIED  NEW  DIGK=  0.004  nsj 

Illustration  17 


PLEASE 


PLEASE 


PLEASE 
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ADPALT  al 


IS  ALT.  PROB?  (YES  OR 

NO):  NO 

ALTERNATIVE  al 
a » : • 5 

SPECI  F. 

rt  «:5 

rc  ■ : 4 0 

nOMMAflD?  : ^nnALT  a 2 

IS  ALT.  PROR?  (YES  OR 

MO):  NO 

ALTERNATIVE  a 2 
a « : .4 

SPECIF. 

rt  ®:4 

rc  ■ : 4 0 

COMMANP?  ; A P HALT  a 5 

IS  ALT.  PROB?  (YES  OR 

NO):  NO 

ALTERNATIVE  a 3 
a » : . 7 

SPECI  F. 

rt  *:6 

rc  * : 4 0 

COMMAND?  : APPALT  a 4 

IS  ALT.  PROB?  (YES  OR 

NO):  NO 

ALTERNATIVE  a 4 

SPECIF. 

a . 8 

rt  =:7 

rc  » : 4 0 

rOMflANP?  : 


Illustration  18 
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EVU  a 


al 

0.279 

a2 

0.191 

a3 

0.501 

a4 

0.F41 

COMMAND?  : 

E’/AL  response 

al 

0.409 

a2 

0.511 

a3 

0.315 

a4 

0.223 

COMMAND?  : 

Illustration  19 


INDIF2  a response 

INPUT  UTILITY  VALUES 

INPUT  2 INDIFFERENCE  PAIRS  PLEASE 

• 

• 28  .III  .iq  .51 
.315  .fi4  .23 


BIGK=  l.R5n/K(?i 
IMHIF  PMR  yields  iriFO  ADOUT 
RPL  K CHECK.  CHRREMT  RAT  in 
IMPLIED  RATIO  = 1.345 

COMMAND?  : KSET  ;;rochow 
» ;.25 

response  = :.34 


) 

REI.  K'S 
respons^k 


BinK=  4,824 
COMMAND?  : 


TO  a = 1.000 


Illustration  20 


INDIFl  a respon.se 

INPUT  UTILITY  VALUES 

INPUT  AN  INDIFFERENCE  .^A I R PLEASE 

.501  .315  .64  .228 


IMPLIED 

IMPLIED 

COMMAND? 


NE'.I  K'S  FACTOR(S) 
NEW  BIGK=  5.239 


9.976  ( -2.301) 


Illustration  21 
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ADPALT  a5 

IS  ALT.  PROR?  (YES  OR  l\0) ; NO 

ALTERNATIVE  a 5 SPECIF, 

a ■ : . 4 


rt  ■ :3 

rc  » : 4 0 


COMMAND? 
a3  . 
a4 
aS 

COMflANP? 


EVAL  r.rocha/  a 3 a 4 a 5 
0.?97 
0.298 
0.308 


Illustration  22 


CHANGE  response  K .31 


COMMAND? 

a3 

a4 

a5 

COMMAND? 


EVAL  ^rochow  a3  a4  a5 
0.292 
0.296 
0.293 

CHANGE  response  K .34 


COMMAND?  : KSET  pjrochow  .75 


BIGK=  11.660 
COMMAND?  : EVAL 
a3 
a 4 
a5 

COMMAND?  : KSET 


,'^rochou  a 3 a 4 a 5 
0.262 
0.260 
0.261 

r.rochov/  1 .33533 


BIGK=  4.824 
COMMAND?  ; 


Illustration  23 


GRAD  ;»rochov/ 
al  0.255 

ATTRIO.UTIL.  GRAD  COMP, 
a 0.418 


response  0.454 

rt  0.303 

rc  0.151 


al 

AND  ATTR.  GRAD 
3.9E5E-01 

-4.461E-02 

-1.401E-03 


COMP. 


Illustration  24 
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IS  ALT. 

ADDALr  a 7 

PROB?  (yrs  or  NO) : NO 

ALTERNA7 

IVE  a7 

SPEC 

a 

»:.76 

rt 

«:9 

rc 

* :2 

COMMANH? 

: ADDALT  a 3 

IS  ALT. 

PROB?  (YES  OR  NO) : NO 

ALTERIMT 

IVE  as 

SPECI 

a 

= :.l 

rt 

»:2 

rc 

»:2 

COMMANH? 

: EVAL  srochow 

a7  a8 

a7 

0.338 

a 8 

0.340 

COMMAfJD? 

• CHANOEALT  rc 

a7 

rc 

“ :100 

roMMANn? 

: CHANflEAl.T  rc 

a 3 

rc 

»:100 

COMMAND? 

: EVAL  r'rochoi-/ 

a7  a8 

a7 

0.148 

aS 

0.2  2S 

COMMAND? 

Illustration  25 
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ADPALT  a 9 

IS  ALT.  PROB?  (YES  OR  NO):  YES 

alternative  nT  SPECIF. 

(ItlCL  0 VIP  lOO";)  FOR  a 

INPUT  THE  CUM  FUNC  F(X).  X'S  FIRST  THEN  F(X)'S 

.1  1 

0 1 

INPUT  THE  CUM  FUNC  F(X).  X'S  FIRST  THEN  P(X)'S 
2 9 
• 

0 1 

INPUT  THE  CUM  FUNC  F(X).  X*S  FIRST  THEN  F(X)*S 
2 120 

0 1 

COMMAND?  : EVAL  srochow  a9 
a9  0.281 

COMMAND?  : AOPALT  alO  .5 

ALTERNATIVE  alO  SPECIF. 

COMMAND?  : EVAL  grochow  a9  alO 
39  9.281 

alO  0.232 

COMMAND?  : 


Illustration  26 
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5.2.5  Comments  on  Computer  Time-Sharing  System 
Illustrations 

Illustration  14 

The  INPUT  command  is  used  to  structure  the  multiattri- 
bute utility  function.  "Response"  is  a nested  MUF.  The  DEBUG 
command  shows  the  defaults  present  after  INPUT. 

Illustration  15 

All  the  unit's  are  set  using  the  constant  risk  form. 
Illustration  16 

The  relative  k^'s  are  determined  between  "rt"  nd  "rc" 
using  INDIFl.  Notice  how  INDIFl  can  aid  in  calculation  when 
a Type  II  Question  (see  Section  3.2)  is  asked.  The  KSET  com- 
mand sets  the  relative  k^'s  based  on  the  output  from  INDIFl. 
The  absolute  k^'s  are  not  yet  known. 

Illustration  17 

INDIFl  is  used  to  determine  the  magnitude  of  the  k^'s. 
The  results  show  that  our  current  setting  is  close  to  the  one 
implied  by  these  indifference  points.  The  nested  MUF 
"response"  has  thus  been  assessed. 

Illustration  18 

Several  alternatives  are  set  up  using  ADDALT.  These 
will  be  used  in  assessing  the  scaling  constants  for  the  MUF 
"grochow. " 
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Illustration  19 

The  utility  values  for  "a"  and  "response"  are  evaluated 

for  the  alternatives.  These  will  be  used  in  the  subsequent 

commands;  e.g.,  u (.5)  = .279 

a 

u^(5,  40)  = .409 

alternative  al  is  the  consequence  (.5,  5.  40) 


Illustration  20 

INDIF2  is  used  to  assess  the  relative  between  "a" 

and  "response."  We  must  use  utility  values  in  specifying 
indifference  points  because  "response"  is  a vector  attribute; 
e.g.,  to  specify  that  (.5,  5,  40)'^(.4,  4,  40)  we  say 
(•279,  . 409)'^( . 191 , .511)  (See  Appendix  A,  Section  A. 2). 

The  KSET  command  is  used  to  set  up  the  relative  k^*s  implied 
by  the  output  frc  i INDIF2. 

Illustration  21 

INDIFl  is  used  to  assess  the  magnitude  of  the  k^^*s  for 
the  MUF  "grochow."  The  results  show  that  our  current  settings 
are  reasonable.  The  MUF  "grochow"  is  now  set. 

Illustration  22 

SVAL  is  used  to  rank  the  alternatives.  The  rankings 
here  are  essentially  the  same  as  in  Grochow. 
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Illustration  23 

Some  sensitivity  analysis  is  performed.  The  CHANGE 
command  alters  the  scaling  constant  for  response.  The  alter- 
natives are  evaluated  and  the  rankings  have  changed.  The 
original  model  is  restored  and  the  magnitude  of  the  k. 's  for 
"grochow"  are  changed  using  KSET.  Again,  the  rankings  change 
from  the  original  model.  The  original  model  is  restored. 

Illustration  24 

The  gradient  for  "grochow"  is  calculated  at  the  alter- 
native  a^^. 

Illustration  25 

Two  indifferent"  alternatives  under  the  current 
model  are  set  u?  using  ADDALT.  The  CHANGEALT  command  is  used 
to  alter  the  common  value  of  "rc"  for  the  two  alternatives. 
They  are  evaluated  again  and  are  no  longer  indifferent.  This 
shows  that  tradeoffs  between  "a"  and  "rt"  depend  on  the  level 
of  "rc."  Our  nesting  scheme  has  captured  this  facet  of  the 

problem.  The  tradeoff  value  of  logging  in  is  degraded  by  the 
poorer  "rc." 

Illustration  26 

A probabilistic  alternative  is  input  and  evaluated. 

In  this  case,  uniform  distributions  are  implied  by  the  cumu- 
latives  which  are  input. 
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A1 though  not  shown  on  the  computer  printout,  the  fol- 
lowing table  is  a comparison  between  the  MUFCAP  approximation 
and  the  graphs  in  Grochow  (4).  (The  scales  in  Grochow  [4] 
are  not  easy  to  interpret  and  the  following  uses  my  interpre- 
tation.) 


Consequence 
(a,  rt,  rc) 

*^MUFCAP 

^GROCHOW 

(1,9,2) 

500 

500 

(1,9,120) 

250 

290  (?) 

(1,2,120) 

750 

750 

(.5,9,2) 

221 

250 

(.5,9,120) 

70 

60 

(.5,2,120) 

373 

383 

(.5,2,2) 

524 

494 

(1,5,120) 

500 

490 

(1,5,2) 

750 

740 

(1,2,40) 

807 

915 

(1,9,40) 

306 

282 

Table  5.1 


A Comparison  of  MUFCAP  and  Grochow  Utility  Functions 


1 
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5.3  The  Comparison  of  Dial-A-Ride  Alaorithms 


This  section  presents  elements  of  a MUFCAP  application 
to  decide  between  two  algorithms  used  by  a computer  to  sche- 
dule Dial-A-Ride  service  which  is  a mode  of  transportation 
being  tried  in  certain  cities  today.  The  presentation  is  con- 
fined to  aspects  of  the  application  which  illustrate  further 
features  of  MUFCAP. 


5.3.1  Attributes  for  the  Problem 

The  attributes  of  interest  in  this  section  are  those 
for  which  preferences  are  not  monotonic.  These  include: 

pickup  time  deviation  = the  difference  in  minutes 

between  the  promised  pickup 
time  and  the  actual  pickup 
time 

travel  time  deviation  i the  difference  in  minutes 

between  the  promised  delivery 
time  and  the  actual  delivery 
time 

The  utility  functions  for  these  attributes  were 
assessed  and  input  into  MUFCAP  making  use  of  the  piecewise 
linear  form.  Two  other  attributes  along  with  these  were  used 
in  making  up  the  overall  utility  function  (see  Turnquist  [171). 

The  utility  function  parameters  were  assessed  and 
several  certainty  alternatives  were  evaluated  to  check  that 
the  utility  function  reasonably  represented  the  preferences 
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of  the  person  being  assessed.  For  this  application,  however, 
the  actual  alternatives  to  be  evaluated  were  outputs  from  a 
stochastic  simulation  program.  One  hundred  outputs  for  each 
algorithm  were  evaluated  using  the  utility  function  assessed 
via  MUFCAP.  That  is,  once  the  utility  function  was  assessed, 
it  was  coded  up  in  a separate  program  to  process  the  output 
from  the  simulation  runs.  An  estimate  of  the  expected  utility 
which  was  the  criteria  for  choosing  between  the  algorithm  was 
obtained  by  taking  the  average  of  the  one  hundred  output  eval- 
uations. This  represents  a way  for  evaluating  the  expected 
utility  in  a case  where  the  attributes  ar  not  probabilistic- 
ally independent  of  each  other.  Although  the  whole  evaluation 
was  not  done  through  MUFCAP,  this  method  for  handling  a case 
in  which  probabilistic  independence  did  not  hold  was  not  too 
difficult.  This  was  because  sanritivity  analysis  could  still 
be  fairly  easily  performed  since  the  utility  function  had  been 
conveniently  parameterized  into  the  multiplicative  form  via 
MUFCAP.  It  is  conceivable  that  MUFCAP  could  be  given  an 
option  for  reading  an  output  file  from  a simulation  model  in 
a future  version  of  the  program.  Then  evaluations  could  be 
performed  within  MUFCAP. 

The  results  of  the  evaluation  showed  that  one  algo- 
rithm was  slightly  superior  to  the  other  over  a wide  range  of 
parameter  variations  a*;d  different  simulation  runs.  Currently, 
a more  ambitious  effort  is  being  undertaken  to  assess  public 
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preferences  for  attributes  germane  to  this  problem  as  opposed 
to  one  particular  individual's  preferences. 


5*3.2  Dial~A~Ride  Illustrations 
UNISET  pickdev  PL 
HOW  MANY  PTS.  I?l  DTIL  FWIIC?  : 5 

INPUT  THE  FUNC.,  X'S  FIRST  THEN  U(X)'S 

-30  0 10  15  30 


.’75  1 .75  .5  0 

COMMAND?  ; UillCAL  pickdev 


U( 

30.000)a 

ci.OOO 

U( 

18.000)= 

0.400 

U( 

6.000)= 

0.850 

U( 

-B.000)= 

0.950 

U( 

-18.000)= 

9.250 

U( 

-30.000)= 

0.750 

COMMAND?  : 


Illustration  27 


INVERSE  pickdev 
30.000=IN'/(  0.000) 

27.000=INV(  0,100) 

21.n00=lN'/(  0.300) 

15.000  = IN'/(  o,500) 

11.000=IN'/(  0.700) 

4.000=INV(  0.900} 

0.000=INV(  1.000) 

COMMAND?  : 


Illustration  28 


-68- 


5*3«3  Coininents  on  Dial-A-Ride  Illustrations 
Illustration  27 

A non-monotonic  utility  function  for  pickup  deviation 
IS  input  using  a piecewise  linear  utility  function.  Some 
sample  ut..lity  function  values  are  tabulated  using  UNICAL. 

The  range  of  the  function  was  input  as  30,  -30. 

Illustration  28 

The  INVERSE  function  shows  only  positive  deviations  as 
attribute  levels  having  certain  utility  values.  This  is 
because  MUFCAP,  for  piecewise  linear  forms,  searches  the  range 
from  the  1st  range  value  to  the  2nd  range  value  until  it  finds 
a level  with  the  appropriate  utility.  This  same  feature  holds 
true  when  an  indifference  curve  is  generated.  This  has  no 
effect  on  the  proper  evaluation  of  alternatives. 
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5.4 A Sampling  of  Problems  to  which  MUFCAP  Has  Been  Applied 

This  section  surveys  some  of  the  areas  where  MUFCAP 
has  been  used  in  a preliminary  manner  to  develop  multiattri- 
bute  utility  functions.  In  all  these  applications,  the  var- 
ious commands  and  procedures  already  illustrated  in  previous 
sections  were  employed.  Chapter  6 further  discusses  some  of 
the  things  which  were  learned  from  these  experiences. 

5.4.1  Evaluating  Health  Plans 

Four  attributes  were  formulated  for  evaluating  health 
plans.  These  were  convenience,  quality,  cost  and  personalness 
of  the  service.  Psychometric  measures  were  developed  for  each 
of  the  attributes  and  questionnaire  assessments  were  used  to 
estimate  the  utility  function  parameters.  MUFCAP  was  then 
used  to  calculate  k in  the  multiplicative  form  and  generate 
indifference  curves  between  certain,  attributes  (see  Hauser  and 
Urban  [ 6 ] ) . 

5.4.2  Evaluating  Policies  for  Dealing  with  Prostitu- 
tion in  the  Boston  Area 

A class  project  in  a decision  analysis  course  at  MIT 
involved  evaluating  five  options  for  dealing  with  the  question 
of  legalizing  prostitution  in  the  Boston  area.  These  options 
were  strict  prohibition,  toleration  or  benign  neglect,  regula- 
tion of  prostitution,  licensing  of  individual  prostitutes  and 
decriminalization.  The  attributes  wore  chosen  to  reflect  the 
prostitute's  position,  the  public  attitude,  the  economics  of 


70- 


the  options,  the  criminal  justice  system’s  opinion  of  the 
options  and  the  political  implications  of  the  choices. 

The  class  divided  into  groups  which  concentrated  on 
the  specific  attribute  areas  defined  above.  The  groups 
assessed  expected  utility  values  for  their  individual  attri- 
butes for  each  option.  Pseudo-attributes  consisting  of  the 
five  attribute  areas  each  measured  by  a utility  value  on  a 
linear  scale  from  0 to  1 were  then  input  into  MUFCAP.*  A sen 
sitivity  analysis  concerning  ranking  of  the  options  was  then 
performed  on  the  magnitude  of  the  k.’s.  It  showed  that  regu- 
lation was  the  preferred  policy  for  the  particular  relative 

k.'s  used  in  this  problem  over  a large  range  of  their  magni- 
tudes. 


This  application  illustrates  how  a complex  problem 
can  be  subdivided  into  smaller  problems  and  the  outputs  from 
these  combined  in  an  overall  utility  function.  In  some  cases, 
the  overall  decision  maker  may  not  be  familiar  with  the  speci- 
fic attributes  used  to  represent  the  objectives  of  a particu- 
lar area  or  group,  if  he  has  a "feel"  for  associating  utility 
with  that  group's  preferences,  however,  he  may  be  able  to 
estimate  the  scaling  constants  and  conduct  reasonable  sensiti- 
vity analyses  in  a manner  analogous  to  what  was  done  in  the 
class  project  on  prostitution. 
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5,4.3  Evaluating  Police  Dispatching  and  Assignment 
Policies 

Attributes  for  evaluating  police  assignment  and  dis- 
patching strategies  include  cost  per  person  per  year,  response 
time  to  various  priority  calls  and  distribution  of  the  work- 
load among  the  different  police  units.  While  models  have 
been  formulated  to  predict  what  workloads  and  response  times 
will  result  from  implement!  g certain  strategies,  work  is  just 
beginning  on  evaluating  the  tradeoffs  between  the  various 
attributes  in  the  problem  which  go  into  deciding  upon  a strat- 
egy. MUFCAP  is  now  being  used  in  preliminary  attempts  to 
structure  a utility  function  for  such  strategy  evaluations. 

5.5 Other  Problem  Setting  Amenable  to  MUFCAP 

Many  problems  which  can  be  cast  as  multiobjective 
decision  making  problems  involving  risk  might  be  amenable  to 
analysis  using  MUFCAP.  This  section  presents  some  examples 
of  current  problems  and  how  they  might  be  structured  for 
MUFCAP  analysis. 

5.5.1  Nuclear  Power  Plant  Siting  and  Setting 
Standards  for  Air  Pollution  Control 

This  subsection  mentions  two  areas  which  have  been 
formulated  as  multiattribute  decision-making  problems  in  the 
literature.  In  Keeney  and  Nair  (10!»  general  objectives  are 
described  for  a nuclear  power  plant  siting  decision.  These 
include  minimizing  environmental  damage,  maximizing  human 
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health  and  safety,  providing  quality  service  for  the  customer 
and  maximizing  the  economics  of  the  company.  Explicit  attri- 
butes might  be  level  of  radiation  per  person  for  human  safety 
and  service  interruption  in  days  for  quality  of  service  to  the 
customer. 

Keeney  and  Ellis  ( 1 ] descrJbe  the  decision  problem 
faced  by  New  York  City  in  legislating  acceptable  levels  for 
sulfur  content  in  fuel  to  be  consumed  by  industry.  The  prob- 
lem is  organized  in  detail  into  a multiattribute  utility  func- 
tion structure  including  attributes  which  reflect  such  objec- 
tives as  the  cost  to  the  city  of  any  plan,  and  effects  on  the 
health  of  the  residents. 

In  both  these  cases,  good  descriptions  of  how  to  for- 
mulate the  problem  are  available.  The  actual  assessment  in 
detail  or  implementation  of  the  formulations  appear  to  be 
possible  through  the  use  of  MUFCAP. 

5.5.2  Anti-Stagflation  and  Energy  Policy  Decisions 

Two  of  the  most  important  multiobjective  problem  areas 
facing  the  United  States  are  how  to  deal  with  the  economic  and 
energy  crises  currently  plaguing  the  country.  A crucial 
aspect  in  these  problems  has  been  deciding  what  tradeoffs  to 
make  between  apparently  competing  objectives. 

In  the  economic  area,  some  of  the  measures  for  objec- 
tives include  the  unemployment  rate,  the  consumer  price  index 
and  growth  in  the  GNP.  The  energy  area  includes  cost  of  fuel 
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and  degree  of  dependency  upon  other  nations.  In  addition,  the 
problem  of  sharing  the  burden  equitably  among  the  different 
groups  in  the  United  States  such  as  labor,  management,  minori- 
ties,  lower,  middle  and  upper  classes,  residents  of  certain 
geographical  regions,  social  security  recipients,  etc.,  lead 
to  explicit  consideration  of  the  tradeoffs  between  these  dif- 
ferent groups  in  trying  to  decide  upon  a policy. 

These  problems  appear  to  be  very  difficult  and  a for- 
mal analysis  such  as  could  be  attempted  with  MUFCAP  might  shed 
some  light  on  comparing  alternative  solutions.  Perhaps  as 
important,  differences  of  opinion  concerning  tradeoffs  among 
the  objectives  might  also  be  clarified. 

5.5.3  Mu3 tiobiective  No-Risk  Contexts 

In  situations  where  no  uncertainty  is  present,  multi- 
attribute utility  theory,  of  course,  is  still  valid.  In  these 
situations,  however,  the  theory  of  value  functions  (ordinal) 
rather  than  utility  functions  (cardinal)  are  applicable  as 
well.  With  three  or  more  attributes,  preferential  indepen- 
dence implies  that  an  overall  value  function  exists  which  is 
a weighted  sum  of  the  individual  value  functions  assessed  over 
the  attributes.  How  one  assesses  value  functions  as  opposed 
to  utility  fun.ctions  will  not  be  discussed  here. 

MUFCAP,  while  designed  to  imple.ient  utility  theory, 
can  nevertheless  be  used  to  implement  a value  function 
approach  to  a problem.  The  value  functions  for  the  individual 
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attributes  are  input  as  if  they  were  utility  functions  using 
the  UNISET  conunand.  The  scaling  constants  are  input  using  the 

KSET  command  and  the  overall  "value"  function  is  deliberately 
made  additive  also  using  KSET. 

MUFCAP  can  then  be  used  to  evaluate  alternatives  or 
generate  indifference  curves.  Different  functions  based  on 
the  preferences  of  different  people  can  be  compared  using 
MUFCAP 's  filing  capability  and  sensitivity  analysis  varying 
the  scaling  constants  and  value  functions 


can  also  be  tried. 
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6.  AREAS  FOR  IMPROVEMENT  AND  FUTURE  RESEARCH 

This  chapter  discusses  various  improvements  which 
might  be  made  to  MUFCAP.  Many  of  these  were  anticipated  in 
the  sense  that  MUFCAP  should  be  considered  a first  edition  or 
a basis  on  which  to  improve.  In  addition,  through  the  use  of 
MUFCAP,  other  new  ideas  for  routines  and  commands  emerged. 

Besides  technical  improvements  which  can  be  made  to 
the  program,  several  theoretical  and  practical  issues  concern- 
ing types  of  assessment  questions  arose  during  the  course  of 
testing  and  using  MUFCAP.  These  issues  are  also  discussed  in 
this  chapter. 

6.1 Ideas  for  Improving  MUFCAP  as  a Computer  Program 

MUFCAP,  being  a computer  program,  can  be  improved  in 
the  ways  that  computer  programs  are  generally  improved.  These 
encompass  four  general  areas. 

The  first  would  be  more  testing  and  debugging  of  the 
existing  routines.  Currently,  a bug  exists  in  the  LOTTERY 
command  which  was  intended  to  perform  a particular  calculation 
when  there  is  a 2-consequence  lottery..  This  bug  can  be 
easily  corrected  when  a later  version  is  compiled,  hopefully 
including  more  than  just  the  fix  for  this  bug. 

The  second  area  concerns  better  program  documentation. 
In  programming  MUFC/^P,  less  attention  was  paid  to  documenting 
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routines  as  opposed  to  getting  them  to  work  properly.  Hope- 
fully, this  thesis  and  the  program  listing  are  sufficient  for 
a knowledgeable  programmer  to  successfully  modify  MUFCAP.  In 
addition,  the  documentation  for  program  usage  could  be  expand- 
ed into  a more  extensive  user's  manual  should  MUFCAP  ever 
attain  the  status  of  a standard  package  for  decision  analysis. 

A third  improvement  involves  making  the  program  more 
"fail-safe”  for  the  user.  Many  precautions  have  already  been 
taken  to  "protect"  the  user  against  leaving  out  necessary 
input  or  making  input  mistakes.  There  remains  room  for  im- 
provement, however.  One  special  area  concerns  generating  an 
indifference  map  involving  an  attribute  with  a risk  averse 
exponential  form.  With  this  form,  there  is  a limit  to  the 
utility  one  could  obtain  even  if  one  had  an  infinite  amount  of 
a desirable  attribute.  If  an  indifference  point  is  given  and 
another  is  desired  having  less  of  one  desirable  attribute  but 
more  of  the  risk  averse  one,  it  is  possible  that  no  amount  of 
that  attribute  will  make  the  new  point  indifferent  to  the  old 
one.  In  this  case,  MUFCAP  tries  to  extrapolate  by  taking  the 
log  of  a negative  number  causing  one  to  exit  from  the  program. 
Thus,  one  should  always  save  the  status  periodically  so  in 
case  one  is  forced  to  exit  from  MUFCAP,  the  program  can  be 
invoked  again  and  the  status  restored. 

Finally,  the  output  could  be  made  more  aesthetic  and 


easy  to  understand.  This  improvement  is  a necessary 
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complement  to  having  better  documentation. 


► 


6. 2 Expanding  Old  and  Adding  New  Routines 

Several  ideas  for  better  routines  concern  the  areas  of 
generating  indifference  curves  more  automatically,  expanding 
the  number  of  available  scalar  attribute  utility  function 
forms,  providing  an  easier  way  of  specifying  probabilistic 
distributions  and  providing  for  analysis  of  alternatives  where 
probabilistic  independence  need  not  be  assumed.  There  is  also 
the  area  of  more  automatic  sensitivity  analysis. 

One  should  be  able  to  generate  an  indifference  curve 
between  two  attributes  which  are  preferentially  independent  of 
all  the  other  attributes  after  obtaining  two  sets  of  indif- 
ference pairs.  Currently,  this  can  be  done  in  MUFCAP  in  three 
stages.  First,  INDIF2  is  used  to  obtain  the  relative  scaling 
constants  and  k in  terms  of  one  of  the  scaling  constants. 

Then  KSET  is  used  with  the  OVERIDE  option  to  set  one  scaling 
constant  arbitrarily,  the  second  in  terms  of  the  first,  and 
k in  terms  of  the  first.  Then,  IMAP  is  used  to  generate 
indifference  curves.  This  procedure  is  one  which  is  often 
requested  because  indifference  curves  are  a valuable  source 
of  feedback.  A needed  improvement  would  be  to  have  INDIF2 
stay  in  an  indifference  curve  generating  mode  and  automati- 
cally generate  indifference  curves  for  the  user  rigVit  after 
input  of  the  indifference  pairs.  This  should  be  fairly  easy 


\ 
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to  implement.  (Alas,  a computer  program  must  be  limited  to 
some  extent  so  a version  can  finally  be  produced.) 

MUFCAP  has  three  scalar  attribute  utility  forms  and 
more  could  be  added.  These  might  include  decreasingly  risk 
averse  forms  based  on  a single  parameter  which  are  very  con- 
venient to  assess  or  multi-parameter  forms. 

Currently,  specifying  probabilistic  alternatives, 
especially  for  a many-attribute  problem  is  laborious.  More 
automatic  setups  of  these  alternatives  are  possible.  Sugges- 
tions include  setting  all  attributes  with  uniform  density 
functions  over  their  ranges  automatically  or  setting  them  all 
with  normal  distributions  about  their  centers  and  having  the 
range  limit?'  be  several  standard  deviations  away.  Also, 
having  set  up  a probabilistic  alternative,  one  should  be  able 
to  copy  it  into  another  alternative  and  then  have  the  ability 
to  change  a particular  component.  A method  of  handling  pro- 
babilistically dependent  alternatives  has  already  been  dis- 
cussed in  Section  5.3.  Another  improvement  would  be  provision 
for  discrete  probability  functions  for  the  scalar  attributes. 

Presently,  in  doing  sensitivity  analysis,  a user  must 
input  the  parameter  changes  and  then  evaluate  alternatives. 

The  program  could  be  made  to  vary  a parameter  over  a range 
and  automatically  evaluate  alternatives,  or  generate  other 
feedback.  This  would  enable  the  user  to  perform  sensitivity 


analysis  more  rapidly. 
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6.3 Making  MUFCAP  Easier  to  Use 

MUFCAP  requires  an  intermediate  "decision  analysis 
person"  to  operate  the  program,  ask  assessment  questions,  and 
discuss  the  feedback  implied  by  the  output.  The  program  might 
be  upgraded  to  (a)  "prompt"  what  assessments  should  be  made  at 
various  stages  of  the  MUF  development,  and  (b)  print  more 
interpretation  about  what  certain  output  numbers  mean.  More 
will  be  mentionc:i  in  this  vein  in  later  sections  of  this 
chapter. 

To  develop  an  interface  dialogue  that  the  program 
would  *')e  completely  self-explanatory  to  decision  m.akers  in  any 
field  would  take  a lot  of  testing  and  v/or;c.  This  might  not  be 
desirable  either  since  discussion  with  a decision  analyst 
should  not  necessarily  be  avoided.  I have  found  that  users 
not  "immersed"  in  multiattribute  utility  theory  were  neverthe- 
less able  to  "order  me"  in  rapid-fire  succession  about  what  to 
do  next.  Setting  up  the  initial  model  is  the  hardest  part. 

But  sensitivity  analysis  should  be  fairly  :*asy  for  a "layman" 
once  he  is  reasonably  satisfied  with  the  initial  model. 

Another  suggestion  has  been  to  put  a graphics  capabi- 
lity into  MUFCAP.  This  would  enable  the  program  to  draw 
utilicy  functions  -^nd  indifference  curves  displaying  their 
shape  to  the  user.  Using  a MUFCAP  with  graphics  would  be 
more  stimulating  in  th,?t  information  would  bo  presented  to  the 
u'*er  in  a more  concise  manner.  Graaient  vectors  might  even 
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be  presented  on  a representation  of  a utility  surface.  Also, 
changes  to  utility  functions,  indifference  curves  or  parame*- 
ters  could  be  input  via  a light  pen  or  a joystick  cursor 
enabling  the  user  to  conduct  sensitivity  analysis  with  his 
hand.  An  advantage  of  the  non-graphics  current  package  is 
that  it  can  be  run  on  a portable  terminal. 

6.4  Assessment  Question  Issues 

Although  MUFCAP  is  a definite  aid  in  MUF  assessment, 
a great  deal  of  discussion  and  patience  is  still  necessary  to 
solicit  accurate  information  from  the  decision  maker.  The 
results  output  by  MUFCAP  are  completely  based  upon  the  input 
information.  In  the  early  use  of  the  package,  it  was  tempting 
to  input  numbers  which  were  not  reasonably  arrived  at  just  to 
see  some  output  from  the  package.  The  output  was  often  non- 
sensical from  the  viewpoint  of  certain  assumptions  about  the 
multiplicative  form.  For  example,  if  two  pairs  of  indiffer- 
ence points  are  input  to  INDIF2,  MUFCAP  essentially  solves 
simultaneous  equations  of  the  form  Ax^^+By^^+CXj^y^  = 
Ax2+By2+Cx2y2  where,  for  the  multiplicative  form,  A corres- 
ponds to  k^,  B to  kj,  and  C to  kk^k ^ . In  solving  these 
equations,  however,  arbitrary  input  can  lead  to  arbitrary 
values  for  k and  k^  in  terms  of  k^,.  For  example,  sometimes 
the  implied  k is  equal  to  -2/k^^  which  is  not  allowed  for  the 
assumptions  of  the  multiplicative  form  as  defined  in  Kecney(9l 
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since  it  is  less  than  -1.  When  this  happens,  new  pairs  of 
indifference  points  should  be  input. 

Besides  leading  to  nonsensical  output,  certain  forms 
of  indifference  pair  inputs  can  given  very  inaccurate  results. 
Indifference  questions  involving  extreme  attribute  levels  are 
very  difficult  to  consider.  However,  indifference  questions 
involving  consequences  which  are  not  very  different  from  each 
other  in  terms  of  attribute  levels  can  give  very  inaccurate 
results.  This  is  because  it  is  hard  to  discriminate  between 
what  is  preferred  and  what  is  indifferent.  The  best  questions 
seem  to  be  those  in  which  the  indifference  points  are  spread 
about  the  middle  of  the  attribute  ranges  and  in  which  attri- 
bute amounts  vary  halfway  between  the  middle  and  extreme  end  ' 
of  the  range.  Also,  specifying  two  indifference  pairs  which 
share  a consequence  point  in  common  [e.g.,  (aj^,b^)  ^ 
and  (aj^,b^)  ~(a2,b^)l  seem  less  prone  to  giving  nonsensical 
results. 

In  using  MUFCAP,  certain  indifference  pairs  appear  to 
be  more  "robust"  than  others  in  terms  of  the  implied  relative 
scaling  constants.  For  example,  the  type  II  question  men- 
tioned in  Section  3.2  is  very  robust  in  the  sense  that  if 
((x.*,  x^®)^(x.®,  X.*)]  implies  certain  relative  scaling 

^ j ^ J 

constants,  [(x^^'  + 6x^  ,x^® ) (x^®  , x ^ ' ) ] implies  almost  the 

same  relative  scaling  constants  provided  6x^  is  small  com- 
pared to  the  range.  This,  however,  is  not  always  the  case 
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when  INDIF2  is  used  with  two  sets  of  indifference  pairs. 
In  oases  where  the  difference  in  the  consequences  is  relative- 
ly small  and  it  appears  as  if  one  of  the  scaling  constants  is 

more  than  twice  the  other,  a 6x.  which  is  small  can  lead  to 

1 

large  changes  in  the  implied  relative  scaling  constants. 

Fortunately,  one  can  test  the  robustness  of  the  rela- 
tive scaling  constants  implied  by  two  sets  of  indifference 
pairs  using  MUFCAP.  One  merely  varies  one  of  the  attribute 
amounts  by  a small  percentage  and  observes  if  the  implied 
relative  scaling  constants  are  vastly  different  from  those 
implied  by  the  original  sets  of  indifference  pairs.  A nice 
improvement  to  MUFCAP  would  be  for  the  program  to  automati- 
cally test  the  robustness  of  certain  inputs  by  performing  the 
appropriate  variations  and  displaying  the  results  for  the 
user.  More  about  this  will  be  discussed  in  the  next  section. 

6.5  Areas  for  Future  Research 

One  area  for  future  research  concerns  the  specification, 
from  a theoretical  point  of  view,  of  assessment  questions 
involving  indifference  pairs  which  are  "robust"  as  discussed 
in  Section  6.4  A starting  point  might  be  to  examine  the 
indifference  curves  which  are  hyperbolas  in  the  utility  plane 
X . (See  Appendix  E.)  We  could  imagine  having  three 
points  on  an  indifference  curve  and  then  displacing  one  of 
the  points  and  plotting  a new  indifference  curve.  H<^w  much 
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the  new  curve  differs  from  the  old  might  depend  on  the  spread 
of  the  initial  three  points. 

A secon*^  rea  for  examination  is  how  to  interpret 
varying  output  during  sensitivity  analysis.  When  several 
pairs  of  indifference  points  are  input,  the  implied  k is 
often  different.  Interpreting  what  constitutes  a significant 
difference  is  not  very  precisely  defined.  For  example,  is  a 
k = -.50  significantly  different  from  a k = «.80.  Where  the 
relative  scaling  constants  are  concerned,  variations  here  are 
directly  related  to  the  size  of  the  differences  in  attribute 
amounts  necessary  to  maintain  certain  indifference  relation- 
ships. But  where  k is  concerned,  it  is  difficult  to  tell 
where  the  differences  will  be  because  k = -.50  as  opposed  to 
k = -.80. 

MUFCAP  can  be  used  to  empirically  examine  what  differ- 
ences result  when  certain  variations  are  perceived  in  the 
value  of  k.  In  addition  to  aiding  in  such  sensitivity  analy- 
sis, MUFCAP  might  also  aid  in  researching  the  area  of  robust 
assessment  questions  and  interpreting  what  constitutes  signi- 
ficant variations  in  parameters  implied  by  the  ansr^ers  to 
assessment  questions. 

A third  topic  for  future  research  would  be  methods  of 
verifying  preferential  and  utility  independence  assumptions. 

In  order  to  use  the  multiplicative  form,  we  must  test  that  the 
appropriate  independence  assumptions  are  satisfied.  This  can 
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be  done  by  asking  a lot  of  tradeoff  questions  and  lottery-type 
questions  (see  Keeney  [8]).  It  can  often  be  laborious  to 
rigorously  verify  the  requisite  assumptions,  however. 

MUFCAP  provides  another  means  for  testing  preferential 
independence.  If  tradeoffs  between  attributes  i and  j imply 
a negative  k,  but  tradeoffs  between  j and  1 imply  a positive 
k,  then  obviously  the  set  of  attributes  i,  j and  1 cannot  be 
combined  into  a single  multiplicative  form  and  are  not  pre- 
ferentially independent.  Earlier  in  this  section,  we  dis- 
cussed the  problem  of  what  constituted  a significant  differ- 
ence in  the  value  of  k implied  by  indif i.erence  pair  inputs. 

If  this  were  known,  preferential  independence  could  be  tested 
by  seeing  if  several  indifference  pair  inputs  implied  the  same 
k within  a certain  "confidence  interval."  If  so,  we  could 
assume  more  confidently  that  preferential  independence  was 
indeed  present. 

6.6  Summary  of  the  Chapter 

This  chapter  discussed  a variety  of  areas  for  improving 
MUFCAP  and  for  future  research.  These  included  improving  and 
further  documenting  the  computer  code  and  expanding  and  adding 
new  routines  to  improve  feedback  and  make  specifications 
easier.  The  issues  in  asking  the  "best"  kind  of  assessment 
questions  were  discussed.  These  included  asking  questions 
which  would  have  "robust"  answers  and  not  yield  results  too 


sensitive  to  small  deviations  in  the  answers.  Areas  for 
future  research  concerned  these  issues  of  robust  yet  reason- 
able assessment  questions,  how  to  interpret,  in  a statistical- 
like  fashion,  variations  in  parameters  implied  by  certain 
indifference  pair  inputs  and  further  ways  of  verifying  certain 
independence  assumptions. 
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1.  SUMMARY  AND  CONCLUSIONS 

This  chapter  sunmarizes  the  main  aspects  of  the  compu- 
ter package  MUFCAP,  The  current  version  provides  the  basic 
features  necessary  to  assess  and  use  multiattribute  utility 
functions  on  complf.x  decision  problems.  In  particular,  it  per- 
mits one  to  use  realistic  and  simple  questions  in  assessing 
the  decision  maker's  preferences,  in  addition  to  the  "extreme 
value"  types  of  questions  previously  used  for  computational 
reasons.  MUFCAP  provides  for  (a)  a variety  of  immediate  feed- 
back of  implications  of  the  decision  maker's  responses, 

(b)  evaluation  of  alternatives  and  sensitivity  analysis,  and 

(c)  analyzing  differences  of  preferences  and  judgements  which 
constitute  differing  models  of  the  same  problem  such  as  might 
arise  among  various  individuals  in  a decision-maKing  group. 

The  present  MUFCAP  should  be  considered  a first 
edition,  a basis  on  which  to  improve.  In  this  regard,  many 
possible  improvements  have  been  suggested  in  this  thesis  such 
as  new  routines  for  (a)  providing  more  readable  output,  per- 
haps even  graphical  displays,  (b)  promoting  easier  feedback 
such  as  more  automatic  computation  of  the  implications  of 
certain  input,  and  (c)  providing  more  aid  to  the  user  as  to 
what  to  do  next.  In  addition,  areas  of  research  were  sug- 
gested concerning  what  kind  of  assessment  questions  are  the 
best  to  pursue  with  respect  to  the  properties  of  being 
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reasonable  to  answer,  and  having  parameter  implications  not 
overly  sensitive  (i.e.,  robust)  to  the 


precision  of  the  answer. 
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APPENDIX  A 

LIST  OF  MUFCAP  COMMANDS  WITH  BRIEF  DESCRIPTIONS 


A.l  Notation  and  Command  Descriptions 


CE 

MUF 

UNIF 


Certainty  equivalent 
Multiattribute  Utility  Function 
Uniattribute  (scalar  attribute)  utility 
function 

Brackets  indicate  the  options  which  may  be 
chosen.  No  option  needs  to  be  selected. 
Parentheses  indicate  that  a choice  must  be 
made  among  the  options  given. 


INPUT  name  - Inputs  the  structure  of  the  multiattribute 

utility  function  to  be  referred  to  by  *name.*  The  dialogue 
requests  names  for  the  attributes  and  their  ranges.  Ranges 
for  attributes  over  which  preferences  are  monotonic  should 
be  input  with  the  least  desirable  end  of  the  range  first. 

A vector  attribute,  (and  hence  a nested  MUF)  is  signalled 
by  specifying  a range  whose  lower  and  upper  limits  are  the 
same.  After  INPUT,  the  default  for  all  MUF's  is  the  addi- 
tive form  with  = k^  for  all  i,  j.  The  default  for  all 
UNIF's  is  the  linear  utility  function.  The  user  is  set  to 
calculate  immediately  after  INPUT. 


SAVE  filename  - Saves  the  current  preference  and  alter- 
native specifications  in  file  named  'filename.* 


r 
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READ  filename  - Restores  the  information  which  was 

saved  in  'filename.* 

DEBUG  - Lists  all  the  attributes  in  the 

utility  function  structure  including  their  names,  scaling 
factors,  ranges,  and  UNIF  typer.  (0,  1,  and  2 indicate  res- 
pectively linear,  constant  risk  aversion,  and  piecewise 
linear) . A vector  attribute  has  its  name  and  scaling  factor 
listed  and  is  followed  by  its  component  attributes. 

ADDALT  altname  [factor]  - Initiates  dialogue  to  specify  an 
alternative  to  be  referred  to  by  'altname.*  Either  a pro- 
babilistic or  certainty  alternative  may  be  specified.  If 
the  former  is  the  case,  a piecewise  linear  cumulative  proba- 
bility distribution  is  requested  for  each  scalar  attribute. 
(Abscissa  values  for  the  cumulative  are  input  in  ascending 
order.)  The  option  'factor*  is  a number  which  sets  all  of 
the  scalar  attributes  at  the  factor  level  of  their  ranges, 
e.g.,  if  factor  = .1,  all  the  scalar  attributes  are  set  at 
one-tenth  of  the  way  from  the  1st  range  value  to  the  2nd 
range  value. 

DROPALT  altname  - ..emoves  the  alternative  'altname' 

from  the  status. 

EVAL  uname  (A,  B,...j  - Evaluates  the  alternatives  A,B,..., 


using  the  ut-’ lity  function  associated  vath  'uname.' 


If  no 


alternatives  are  specified,  all  alternatives  in  the  status 
are  evaluated  and  the  results  listed. 

UNISET  uname  (LIN,CR,PL)  - Sets  the  scalar  attribute 

utility  function  associated  with  * uname*  to  linear,  constant 
risk  averse,  or  piecewise  linear  form.  For  the  piecewise 
linear  form,  the  abscissa  values  are  input  in  ascending 
order . 

KSET  mname  (factor , ADD, OVERIDE]  - Sets  the  scaling  factors 
for  ♦he  MUF  associated  with  'mname.*  The  number  'factor* 
causes  the  current  scaling  factors  to  be  multiplied  by  that 
number.  The  program  automatically  calculates  the  k asso- 
ciated with  the  new  scaling  factors.  If  ADD  is  specified, 
the  current  factors  are  normalized  to  add  to  1.  The  user 
may  input  k directly  in  response  to  the  final  prompt  by  the 
computer  if  OVERIDE  has  been  specified. 

GRAD  uname  [A,B,...]  Calculates  the  gradient 

components  of  the  utility  function  associated  with  'uname* 
for  all  or  some  of  the  alternatives  A,B,...  . 

INDIFl  unamel  unamel’.  - In  the  unamel-uname2  attri- 

bute plane,  given  relative  k^'s,  (i.e.,  scaling  factors 
with  the  appropriate  ratio  relationship  to  each  other  but 
not  necessarily  the  appropriate  absolute  value)  the  k is 
specified  by  a single  pair  of  indifference  conseguences . 
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INDIFl  requests  a pair  of  indifference  consequences  and  uses 
the  current  k^*s  as  the  given  relative  k^*s.  On  output,  the 
k is  given  along  with  the  factor  by  which  the  current  k^'s 
must  be  multiplied  to  yield  the  k (see  KSET  command  with 
•factor*  option). 

INDIF2  unamel  uname2  - In  the  unamel-uname2  attribute 
plane,  with  scaling  factors  denoted  by  k^  and  k^#  inputting 
two  pairs  of  two  indifference  consequences  each  specifies 
the  ratio  k^/k^  and  k - constant/k^^ . After  INDIF2,  the 
KSET  command  may  be  used  to  fix  k^^,  and  then  k^  and  k in 
terms  of  k^.  The  command  IMAP  can  then  be  used  to  generate 
indifference  curves  in  the  unamel-uname2  plane.  (For  these 
indifference  curves,  the  values  of  k^,  i ^ 1,2,  are 
irrelevant) . 

UNICAL  uname  [n]  - Prints  a list  of  utilities  using 

the  UNIF  associated  with  * uname.*  Once  the  number  n is 
specified,  the  user  supplies  n attribute  amounts  and  the 
program  returns  the  n associated  utilities. 

INVERSE  uname  (n)  - Prints  a list  of  attribute  amount/; 

associated  with  utilities  using  the  UMIF  *unarne.*  Onco  the 
number  n is  specified,  the  user  supplies  n utility  amounts 
of  'uname*  and  the  program  returns  the  n associated  attri- 
bute levels.  If  n is  not  specified,  the  program  has  a 
default  printout. 
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CHANGEALT  uname  altname  - Routine  to  change  the 

•uname*  attribute  component  of  the  alternative  * altname* 
without  cha!iging  the  other  components. 

CHANGE  uname  (NAME,K,  RANGE;  param  - Routine  to  change  the 
n€une  or  scaling  factor  or  range  of  the  attribute  * uname*  to 
param.  When  the  range  is  changed,  param  is  not  required. 

The  program  requests  respecification  of  the  UNIF  type  when 
the  range  is  changed.  When  the  name  is  changed,  param  must 
not  be  left  blank. 

ALTLIST  - Lists  the  current 

alternatives.  The  probabilistic  alternatives  are  listed 
with  their  CE  equivalent  components. 

DISPLAY  uname  - Displays  the  charac- 

teristics of  the  ucility  function  associated  with  *uname.* 
The  scaling  factors  for  the  attribute  arguments  and  their 
sum  is  listed  for  a MUF  while  the  range  and  type  is  listed 
for  a UNIF. 

FRACTILE  uname  altname  - Displays  the  cumulative 

distribution  for  * uname*  in  the  alternative  * altname.* 

LOTTERY  uname  n - Calculates  the  CE  for 

a lottery  involving  the  scalar  attribute  * uname.*  The  num- 
ber n specifies  the  number  of  possible  lottery  consequences. 
These  are  solicited  with  their  corresponding  porbabilities 
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and  the  CE  is  calculated. 


IMAP  unamel  uname2  - Initiates  a dialogue  to  generate 

an  indifference  ‘curve*  in  the  unamel-uname2  plane.  A point 
through  which  the  curve  will  pass  is  solicited.  Then  values 
of  unamel  are  input  and  the  uname2  values  required  to  main- 
tain indifference  are  output. 

STOP  - Thanks  the  user  for  using  MUFCAP 

and  exits  from  the  program. 


ADDU  unamel  uname2  - Initiates  a dialogue  which  adds 

an  attribute  * unamel*  to  the  argument  list  of  the  MUF  asso- 
ciated with  *uname2.' 

DELU  uname  - Deletes  the  attribute  *unamo* 

from  the  structure. 


SWITCH  uname  uname2  - Adds  current  attribute  * uname* 

to  the  argument  list  of  the  MUF  associated  with  *uname2*and 
deletes  * uname*  as  an  argument  of  the  MUF  to  which  it 
originally  belonged. 


INTERBK  uname  - If  any  attribute  arguments  of  the 

MUF  associated  with  * uname*  is  a vector,  its  utility  func- 
tion is  a nested  MUF  with  its  own  internal  constant  k. 
INTERBK  calculates  the  theoretical  k for  the  nested  MUF 
which  would  make  the  nesting  of  the  inner  attributes 
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unnecessary  and  prints  it  along  with  the  current  internal  k. 

A. 2 Further  Notes  on  INDIFl,  INDIF2  and  IMAP 

The  INDIFl  comma,  i may  b'3  used  with  input  to  a Type  II 
Question  (see  Section  3.2).  It  will  then  give  the  relative 
k^*s  as  output.  An  example  of  this  is  shown  in  Illustration 
16  of  Section  5.2.4. 

For  INDIFl,  INDIF2  and  IMAP,  if  either  unamel  or  uname2 
is  a vector  attribute,  consequences  must  be  input  as  utility 
pairs  rather  than  attribute  value  pairs.  The  utility  for  an 
attribute  value  is  the  result  obtained  when  that  attribute 
amount  (vector  or  scalar)  is  evaluated  using  the  utility  func- 
tion associated  with  the  attribute  name.  An  example  of  this 
is  shown  in  Illustraions  19  through  21  of  Section  5.2.4. 

Unamel  and  uname2  must  be  explicit  arguments  of  the 
same  KUF  when  using  INDIFl,  INDIF2  or  ir4AP.  That  is  to  say, 
(unamel,  uname2)  must  be  preferentially  independent  of  the 
other  attributes. 

Finally,  on  output,  INDIFl  prints  a number  in  paren- 
thesis as  a second  factor  by  which  to  multiply  the  current 
kj^'s.  If  multiplied  by  this  factor,  uhe  new  k^*s  will  not  be 
consistent  with  the  indifference  pair  input.  Hov’ever,  these 
new  kj^'s  will  yield  a k identical  to  that  of  the  new  k^^'s 
derived  by  using  the  non-parenthesized  factor.  In  practice, 
although  not  consistent  with  the  indifference  pair  input,  the 
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"alternative"  k^*s  come  close  to  being  consistent.  Sometimes, 
the  nor*-parethesized  factor  will  yield  k^*s  which  are  not 
allowed  in  the  multiplicative  form;  e.g.,  k^^  > 1 for  some  i. 
When  this  happens,  the  parentesized  factor  can  be  tried 
instead.  Using  IMAP,  with  these  alternate  k^'s,  we  can  see 
if  the  indifference  pair  consistent  with  these  alternate  k.*s 
is  close  enough  to  the  original  pair  used  in  INDIFl  to  justify 
use  of  the  parenthesized  factor. 
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APPENDIX  B 


MUFCAP  PROGRAM  LISTINGS 


«n»CAF:  PPOC  OPTIONS  CNATN) ; 

DCL  DONflT  *?!!*  (1280'?)  BUSED  (DOT®)  i 
OCL  (E«PTY,NOLL,OMSOn?Cf)  PDTt.TTM; 

DCl  (Bir,K,TNIPXP)  'NTFY  FTTK^NM.; 

nCL  WORD(IO)  CHXP  (12)  ,PFt.Hr,  BTT(1),nF  BTT(1); 

DCL  coin \NC  (2S)  CH^P(12)  I NIT  (•  I NPOT*  , ' SX  VE*  , * ®^»D  • , 

•DEBOr.*,  • ADDMT*  , *P"DPALT*,  'FYAL*  ,Mf  U3ET*  , * FS  ET  ' , • G?  A D*  , 
•IMCIfl'  ,•  INDIP2*  ,*  DNICA  » , • INV , 

•CHANGEALT*  , '‘CHANGE*,  'AL  1ST*  ,*  DISPLAY*  F*>ACTTl?»  , 

•LOTTERY*,*  INAP*  ,*STt  t-*  , * DEt.'J  * , * ADDU  * , * SWT  Ten  • , fTHT^pniT  •)  . 
OCl  BOFILE  FILE  P2COFD  3FCUENTIAL; 

OCL  PROC(26)  LABEL; 

DCL  OLIST  (3^'*  OFFSFT(rUfm)  STATIC; 

OCL  OTHABCt,',  STATIC  CHAF(12),  N*JTN  STATIC; 

OCL  HGSAD  STATIC  ,GPAP  (RC)  ST>T£r^j^p  bIT(1); 

OCL  NSflB  (33)  STATIC  , XIN(15),GP  PIT(1); 

DCL  (MAT,  NC,TCALT,IPALT)  STATIC,  .1APT(2)  FIXED; 

OCL  CARD  CHAP(80)  , ANAS2  CHAP  ( 1 2) ,FN A" F CHAP  (12); 

OCL  (J1 ,J2,X1 ,K?)  FIXED, YIN  (15)  ; 

/♦  OBFTNE  AREA  ’OH  BAS’n  ALLOCATTDNS  •/ 

OCL  1 IIST  BASED  (LISTPTP) , 

2 FIRST  CP’SET  (OHNNY) , 

2 MASPALT  (5)  CHAP  (12)  , 

2 MABCALT  (10)  CHA^»  (’2)  , 

2 BODY  AREA  (12R00)  ; 

/*  BOLTIATTRIRnT?  OTILTTY  FONCTIOM  STROCTORF  */ 

OCL  HOFF  PTR  STATIC; 

OCL  1 BOP  BASrn  (BO^P), 

2 CAPK, 

2 BMABE  CPAH  (12) , 

2 M01AT, 

2 S0BAT(12)  , 

3 CHAIMF  CF’SET  (DOBBY) , 

3 SBALIK, 

3 OMIPTR  OFFSET  (DOBBY)  , 

3 OMABf  CHAR  (12)  ; 

/♦  OMIATTHigurP  OTII.ITY  FUKCTTON  STROCTOFE  •/ 

OCl  DMIFP  FTP  STATIC, 

DCL  1 ONIF  BAS’D (MNIpn) , 

2 OLO,  2 MHI,  2 aiTPE  FIX’D, 

.2  CALT(K),  /*  CEPT.  ALTFPNATTYFS  •/ 

3 CATX,  3 FMC, 

2 OXP  (IS) , 2 UYP (IS) , 2 wnp, 

2 PALT(5),  /•  PPOR.  ATT’RNATIY’S  •/ 

3 MF,  3 XF(9)  , 3 CF(9)  , 3 EOP; 

/*  IMITIAII7P  •/ 

ALLOCATE  LIST;  NnTN=0;  Dt»TP»  A DOR  (BODY)  ; NC0n*2S;  GF»*''*D; 

DO  1*1  TO  S;  NA«PAL?|I)=*  *;  FND; 

DO  1*1  TO  10;  NABCALT(i)»»  •;  ^no; 

POT  SKIP  lIST(*rCflnANO  HOPD  AND  FILE  BAB»S  .1057  DE  TV  CA’S*); 

CM  CONfERSTDN  P’GTN; 

OCL  P FLOAT;  IF  ONSOOFC  E*  • orilT*  THEN  DO; 

FOT  SNIP  T.TST  (*  POSSIBLE  STATOS  CHANGE.  ONOO  PARTI.M  OP.*); 


CDC03C1D 

0300002^ 

oooorci? 

OOOCOOUH 

OICO'^C‘^0 

OCO'^D’A? 

00C03C7C 

OOOOOOPO 

coroim 

oimoiG-? 

ODOOGII'^ 
00C0312C 
0000013^ 
000  0314  :• 
CC0C01S3 
C3CCOVO 
0''000  VIC 
00030'<RC 
0000310"' 
0300370'' 

ooc:*»2io 
00000220 
00CC0230 
000^024'' 
0300''75'' 
OJOO''?6'' 
0000027C 
O‘'C002P0 
0300029  ■> 
0000030: 
00000313 
0300332: 
000  70330 
OOCC3340 
00C702S7 
0300O3S? 
0*30003  70 
OOOOOIP'^ 
03000390 
C30CO407 
000074  10 
C3000420 

000OO4  7P 
03003440 
C0C0"'4S0 
030004S.7 

00C07U*»C 

07000480 

03000407 

03COOSOO 

030C0S13 

00003S27 

03C07S30 

03000S40 

C30C35S"' 

00307Sfo 

07O00073 
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0.820225. 11 940. laPCI^P.PLI 
GO  TO  GBTC05;  FMO; 

POT  S^IP  £CIT  (nSSOURCE.  • IS  <OT  k fUlID  W09B?P.», 

•IHPHT  TH?  CORPF.CT  HUHRf!??  :«) 

ICOLC2)  (t.ENr.TH(ONSOtIRCE))  ,\,k)  ; 

GET  LIST(P);  ONSIUPCE^P;  E»in; 

OM  OHCEPINFOPTLE  C*U»TLE) 

POT  SKIP  :iT(»PILES  10ST  PE  ALLOCATED  AND  PILE  NAH'S  IN  CAPS'); 

GO  TO  GIlCJS;  PND; 

/•  COH8AND  PSOCPSSOP  SFCTTON  •/ 

GPTC05: 

PUT  SKIP  LIST('CC1!1ANO? 

CALL  GETLIS»; 

DO  NC*1  TO  NCOR; 

IF  WORD  (1)=CONSAMn  (NC)  TH’M  GO  TC  t>prC(NC);  END; 

POT  SKIP  EDIT  (»CPD(1),'  IS  NOT  A VMTD  COH1AND.', 

• (C0H8AND  WORD  “GST  BP  TV  CAPS)') 

(COM2)  « A,  A,  A)  : GO  TO  GFTCOB; 

PBOC(1):  /•  'INPUT'  */ 

IP  WCFD(2)  ='  ' TVEV  DC; 

POT  SKIP  LIST  ('NAVE  PCS  «UP  PLEASE?  ;'); 

CALL  GETLIVE:  WORD  ( 2)  =WO»D  (1)  ; GO  TC  P10C(1):  ENT>: 

POT  SKIP  LI3TCH0W  VANY  ATTPIDUTRS  AP’’  IN  THIS  PUP?  ;'); 

GET  IIST  (NAT):  IF  NAT  <*0  THEN  GO  TO  GETCON; 

DO  1*1  TO  5;  NAVPALT(T)*'  ';  END; 

DO  1*1  TO  IT;  NAVCALT(I)*'  •;  END; 

BODT*EKPTY;  ALLOCATE  MUH  in  (BODY);  "NA VE=«OPn  (2)  ; 

PIPST*iaPP;  CAPK*0;  NUSAT*NAT;  CALL  GflMOLT;  CALL  SETOPP; 

, 60  TO  GETCOK; 

PROC(2)  ! /*  'SAPP'  */ 

TP  »ORD(2)«'  ' THEN  DO; 

POT  SKIP  LIST('»ILF.  NAH^  FOR  SAP*  PLEASE?  :'); 

CALL  GRTLTNE;  WORD  (2)  *MCPn  (1 ) ; SO  TC  P»OC(2);  END; 

OPEN  PILE  (HD^'ILE)  TITLE  (WORD  (2)  ) OUTPUT; 

NRITE  PILE  (nUPIL'*)  PRCH  (LIST);  CLOSE  PILE  (BUPTL2)  ; 

GO  TO  GETCOH; 

PROC  (3)  : /•  'PEAC  •/ 

IP  WOBD(2)*'  ' TH’^N  DC; 

POT  SKIP  LISTCPILE  NAHE  FOR  P.BAC  PLEAS5’  ;'); 

CALL  GETLINE;  WOVC  (2)  =WO?D  (1 ) ; GO  TO  r'»OC(3);  END;  . 

OPEN  PILF(HDPILP)  TITLE  (WORD  (2) ) INPUT; 

READ  PILE(HUPILE)  INTO  (LIST);  CLOSE  PILE ( BUPIL?) ; 

CALL  SBTOFP;  GO  TO  GETCOfi; 

PROC  (4)  ; /•  'DEBUG'  •/ 

POT  SKIP  EDIT (' STRUCTURE  FOR  ' ,PNAHP)  (COL (5) , A, A) ; 

00  1*1  TO  NUTN; 

aOFP»OLIST(I) ; PUT  SKIP  EDIT (UTNA V!  (I) ,SOBAT  (NSUB  (I) ) . SVALLK) 

(COL  (2)  ,A,P  (H,3)  ) ; I?  SUBAT.CHAI  IP(VS0B  (I)  ) =NriLL  '^REN  DO; 
ONIFp«SOBAT. UNTPTR (NSns  (I) ) ; PUT  SKIP  LIST  (ULO,MHI ,UT Y PE)  ; SND;END; 
CO  TO  GETCOH; 

PROC  (5)  ; /•  'ADDAl  ‘'  •/ 

IF  N0PDC2)  *'  ' THEN  DO; 
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POT  SKIP  lISTf'NAflF  ALTF'MATI7F  PLFAS"?  :•); 

CALL  GFTLTNE;  WO? 0 (2 ) =W OF 0 (1 ) : f.O  TO  P^OCC*^);  F.k’P; 

AMAH2*WOPP  (2)  ; I'*  WOPD(3)-.=  * • THEN  nO;  A?=M»D; 

PPLAG=*1*P;  AP^C  = Uopt  (3)  ; OO  TO  CZ^T;  EN^;  FI.SF 
ANAHEsWOPD  (2)  ; PUT  SKId  LTST(*TS  Al.T.  p-^OR?  JYFS  OF  NO);*): 
REPLY5;  CAtl  07TLINF:  T»  W OPO  (1 ) = ' T E 5 • THEN  PPLAO-*1*B; 

ELSE  IP  WOPO  (1)  N"^*  TMFN  -•  0*  F ; 

USE  DO;  PUT  SKIP  Li!Sl  (•''EPLY  •’tlST  BF  TZS  nr  NO  TN  CAPS  •): 

GO  TO  PSPLYS;  FNO; 

IP  PPIAG  THEN  rn  I=1  -Q  S; 

IP  NAHPALT  (I)  =•  • THEN  DO; 

NAHPALT  (I)  =ANAHF.;  ir\)[T  = I;  GO  TO  PP0C5B;  END; 

TP  1=5  THEN  DC; 

POT  SKIP  EDIT  ('ONLY  5 PPOP.  ATTEPNATIVES  ALLOWED.  •, 

•ONE  SROHLO  P?  DELETED.  (STATUS  CAN  PE  SAVED,  TOC)*) 

(COL  (2)  , A,  A)  ; 

GO  TO  GETCCP;  rflD;  FND; 

ELSE  CEP'T*;  DO  T=  1 TO  1C; 

IF  NANCALT(I)  =•  • THEN  no; 

NA*1CAIT(I)  =AN.A«!E;  ICALT=I;  GOTO  poccSB;  END; 

IP  1=10  TM’N  CO: 

POT  SKI?  LTSTCONLY  10  CFPT.  ALT,  AUnWED*); 

GO  TO  GETCCf;  ?NC: 

PPOC5B:  PUT  SKIP  ''niT  (•  ALTEPNATTV**  '^ANAPP,*  SPECIF.*)  (A.  A, A): 

00  1=1  TC  NUTN; 
nnPP=OLIST  (T)  ; .1  = NSUB(T); 

IP  SUBAT  (J)  .CHAIN’  = NHLI.  '^HFN  OO; 

OHTPP=SnBAT(J)  . riNIPTR; 

• CALL  ALTCONP;  '‘ND;  ENC;GO  TO  GETCOPI: 

ALTCCHP:  PFOC; 

/♦  NEEDS  I,PPLAG,TC.ALT  OR  IPALT,AE  AND  UKTPP  S^T  •/ 

/♦  SETS  THF  COMPONENT  POP  AN  Al TEPNATTVE  */ 

IP  PFLAG  THEN  DO; 

POT  SKIo  FDIT(MIOW  HANT  PFACTILE  (INCL  0 AND  10C«)  »*0P 

OTNAHE(I),*  (2<=N<=H)  : •)  (A,  A, A); 

GETLYST(N):  PALT ( IP  A f T)  . N?  = V; 

POT  SKIP  IISTCINPHT  the  CUN  PUNC  » ( X)  . X**S  PIPST  THEN  P(X)**S*); 
GET  IISI((XIN(J)  OC  1=1  TO  N)  ) ; 

If  OHKOLO  THEN  DO  7=1  TO  V;  OAL T (IF ALT)  . XP  (.7 ) = 
CnN(N-J»1)-niC)/(URI-TTL0)  ; ; 

ELSE  DO  J = 1 TC  N;  PALT  (I  FALT)  . X P (,7)  = ( XI  N (J) -OLO)  / 

(OHI-OLO) ; FND: 

GET  LTST((XTN(J)  DO  J=1  -^0  N) ) ; 

IF  OHKULO  THEN  DO  J=  1 TO  N; 

- PALT(IFALT)  .C?  (J)  =1 -X 1 N (N-.7^1 ) ; ^ND; 

EISE  DO  .1=1  TC  N;  PAL  I ( T P AT,  T)  . C?  (.1)  * TIN  (.1)  : '’ND; 

/•  ABOVE  INSOPTS  THAT  INTTPNAL  '>EP  C*  CUP  PU  NC  ^ s OKAY  */ 

CALL  ONIFU  (IPALT,ANSJ  ; PALT(I^ALT)  . EUP=  ANS; 

END: 

ELSE  IP  -AP  THFN  DO:  'UT  SKI®  ZDI T( CTNANF (I) , * = : * ) (A, A): 

CETLIST(X);  X=  (X-ULO)/ (UHI-ULO)  ; CALL  UNT  CAL  (X  , A NS)  . 

CALT  (ICALT)  .CALT  = T;  cAI.T  (ICAL^)  , 2UC  = ANS; 

END;  ELSE  no;  CALL  UNTC AL  ( A» AC , ROC ( TCA LT) ) ; 

CALX  (ICALT) =A '^C;  FND;  END  ALTCOHP; 

PFOC  (6):  /*  *OP'^'PALT  •/ 

IP  N0PD(2)  =*  * THEN  DO; 
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POT  SKIP  LIST(*N^nF  FC9  ULTEPMATIV*  PI.EASP?  :•); 

CKll  CETLTM!;  BOP P (2 ) =«0?0  (1 ) : no  tq  PPOr<#i);  END; 

IMA1»«WOPP  (2)  ; DO  T*1  TT  10;  TP  H AHC  ALT  (I)  -A  NAfI?  THPM 
KABCALT  (I)  =•  •;  INP;  no  1*1  *no  S;  IP  NAIPALT(l)* 

AMAHE  NABPALT(T)«  • *;  PMD;  GO  TO  GFTCCP; 

PrOC(7)2  /*EVAL  •/ 

CALL  PF0C7A;  GO  G??CC*1; 

PPOC7A:  PROC;  TP  WORD  (2)  *PNAf1E  THEM  DO;  nnpp*PI?ST; 

DP*«0«P; 

CO  TO  PRCC7C;  FND; 

DO  1*1  TO  NOTN;  T»  OTNA^in  (I)  =Bopp  (2 ) THFN  GO  TO  PR0C7B:  EMT); 
PUT  SKIP  LIST  ('ATTFin  M 3*^  POUHD*):  =ETM®N: 

PK0C7B:  BtlPpsOLIST  (T)  ; I?  SUBAT  (SSPB  (I)  ) .THATV^s-^faL  THEN  DO; 
OKIPP*SOBAT  (NSUB  (I)  ) . GNICto;  END;' 

ILTE  CO;  0?*»0*B;  pnppssnB-^r  (NSHB  (T)  ) .CHATNF;  »»ND; 

PH0C7C:  IP  BOPD(3)**  • TH!N  GO  TO  PFOC7F; 

J*2; 

PPOC7D:  J*J^1;  IP  WOPn(J1  *•  ’ THFH  PE'^U»»J; 

PPlAGa*0*B;  DO  1=1  to  10;  ICALT  = T;  MA  PCM.T  f I)  =WOPD  ( J) 

TBEH  DO;  I»  fl?  THEN  A »IS*CALT  (T)  . EUC;  ELSE  CAL.  K-ULTF.7  (AMS)  ; 
POT  EDTT(MAPCAL?  (I)  ,AN5)  (COL  (2)  .A  , X (1)  ,F  (B  , , ) ; 

IP  (GPC(-OP))  THEM  CALL  GFTG^aD; 

GO  TO  PPOC7D;  END;  SMn; 

PPIAG*«1»P;  DC  1=1  TO  5;  IOAT.T=I;  IP  MABPALT  (I)  ^ BOFD  (J)  THEM 
DO;  TP  THEN  AMS  = PALT  (T)  . EGP;  PT.SP  CALL  riH.TFV  (AMS)  ; 

POT  EDIT(NAflPALT(I)  ,AMS)  (COL  (2)  , A , X (1)  ,F  (B,  3) ) ; 

IF  (CPC  (-HP))  THEN  CALL  (.ETC-PAR;  GO  TO  PRCC7D; 

EHD;  END; 

PF0C7E:  PPLAG=*n»B;  Do  1=1  TO  10;  ICAIT=I;  IP  NAWrALT (I) -=•  • THEM 
DO;  IP  OF  THEN  ANS--CALT(T)  .PrjC;  ELSE  CALL  nOLTRV(ANS); 

POT  PniT(NAMCALT  (I)  , ANS)  (COL  (?' ,A  - X ( 1)  ,F  (R , 3) ) ; 

IP  (GPC(-.UP))  THEN  CALL  GETGPAD; 

IHD;PND; 

PPLAG*»1»B;  DC  T = 1 TO  5;  IPALT=I;  TP  NAHP  ALT  (I)  ^.a*  • THEN 
DO;  IF  OF  THEM  AMS=PA  L*^  (I)  , -’flp;  FLSF  CALL  nOLTEV(AMS); 

POT  FDIT(NAnFALT(I)  ,aNS)  (COL  (2)  , A,  X (1)  ,F  (fl,3))  ; 

IF  (GPC(-»OP))  THEM  CALL  G»»Tr,PAD; 

END;  EfJD;  RETO»N;  END  PPOC7A; 

PKOC(R):  /•  •ONISET*  •/ 

CALL  ONIGPT; 

PROCHC:  IP  BOPD(3)=*  * THEM  DO; 

POT  SKI^»  LISTCTTPr?  ;•);  CALL  GETIIMF;  WORD  ( 3)  aVORI)  (1)  ; 

GO  TC  PROCOC;  ^WD; 

IF  HOPO  (3)  a»CP*  THEM  DO; 

POT  SKIP  EDir(*IMPOT  AMT  50-*ia  LO’^TERT  IM  THE  FOPH  OP  . 

• C.I.,01  6 02.  PLEASE*)  (A  ,A)  ; 

GFT  11ST(CF,X  1,  T2)  ; IP  UHI>MLO  THEM  TO; 

ILOaHIN  (fl  ,X2)  ; X‘<I  = NA  X (XI , X2)  ; EMD; 

ELSE  no;  T!C  = BAX  (X1,T2)  ; XHl  = 11  ?MX  1 , X2)  ; END; 
pa(OHT-nLO);  Cn=  (CF"rTTO)/r  ; XLO=  (XLC-ni,0) /P  ; 

IHI*  (fHI-ULC)/P  ; CALL  IINIEXP  (CS,TLO,  XHI , UXP  (1)  ,UTP  (1) ) ; 
OTTPBal;  END; 

ELSE  IP  BOPD  (3)  =*LI  N»  THEM  f?TTPE*0; 

ELSE  IP  BOPO(3)a*PL*  THEN  DO; 
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POT  SUP  LTST(»HO«  HANT  PTS.  TM  OTTl  fOMC?  00CO22Pv0 

G^T  IIST(N);  PUT  STIP  EDIT(*IMPHT  THE  PUNC,,  I*'S  PIPST  000C23C0 

•THEN  U(X)**SM  (A, A);  000n23n 

GIT  IIST((XIN(I)  ?C  1=1  TO  N)):  01002320 

GET  LIST((OYP(T)  DO  1 = 1 TO  N));  0000231'' 

llOr*M;  00002.3JJO 

IF  OHI>OLO  THEN  DO  1*1  TO  N;  OXP  (T)  = |XTM  (T) -010)  / {UHI -ULO)  ; 0^CC2?S0 

END;  ELSE  DO;  DO  i*1  TO  N;  HT**  (I)  = (TIN  (N-t  41) -ULO)  / OnoC2360 

(ORI-OLO)  ; XTN  (N-:^  1)sfIYP  (I)  ; ^0  1*1  TO  N;  0000237^^ 

OTP  (1)*XTN  (I)  ; END;  ESP;  010023P'' 

0TTPB*2;  END;  PLSE  DO;  PUT  SKIP  irST(»UN:»  TYPE  NOT  VALID*);  00002  3^'^ 
GO  TO  OFTCd;  END:  0)00240'' 

/T  GPDATE  EXPFCTPD  UTILITY  POP  ALTSPNATIVFS  •/  0000241^ 

00  1=1  TO  5;  IP  NAHPAIT  (I)-.**  • THEN  CALL  nNIEH  (I  , PA  LT  ( I)  . EOP)  ; 0 30P  2 42  0 
END;  00  1*1  TC  10;  IP  NA "CAL (I i -** • • THEN  CALI  0 )002430 

ONICAL(CALT  (I)  .CALX.CALT  (I)  .FIJC)  ; END;  GO  TO  GFTCON;  00002440 

000024*^0 

PPOC  (9)  : /♦  *KSFT  •/  0?0024fiC 

IF  HOPO  (2) *PNAnP  THEN  DO;  NUPP*P^RST;  GO  TO  PPOr9C:  FND;  0100341C 

00  1*1  TO  NUTN:  IF  OTN  AH  E (I)  = NOR  D (2)  th£K  GO  TO  PROCPH;  PNO;  0G0''24H'' 
PUT  SKIP  LIST  (•  AT'^PIB  HOT  «^OUHO*);  GC  TO  G«’TCO":  O0CO24G3 

PPOC9B:  HUPP*OI.TST  (II  ; IP  SUB  AT  (NSIJB  ( I)  ) . CKA  INP  = N PLL  THEN  DO;  0')002‘^00 

POT  SKIP  LIS?(*ATTPTH  IS  NOT  A NUP*);  GO  TO  GETCOH;  FND;  0''0G2sn 

HOPP*SORAT(NSUB (I)  ) .CHAIN?;  00002S20 

p1B0C9C!  OOC0253  0 

IF  H0RD(3)  **f.DD*  TH’^N  DO;  SUHK*C;  10  T*1  TO  MHHAT;  OOOOPFtiP 

S0HK*S0?1K^SUPAT  (I)  .SPALLK;  END;  *ACTCR=  1. /SUHK  ; END;  0)002*)S') 

ELSE  IP  N0PC(3)  **OVEPIDE»  THFH  DO;  NOPD  (4)  «HOPD  (3)  ; 0)002S60 

tOPO(3)**  •;  END;  0O002570 

* ELSE  IF  H0P0(3)-***  • THEN  PACTOP*  HOP  D (3 ) ; 000025H0 

to  1*1  TO  HOHAT;  030C2590 

If  W0P0(.3)**  • TH«'N  DO;  POT  EDI  T (SOBAT  (I)  . ON  AHE,  • * ;•)  (COL  (2)  rOOOO  2600 
E,A);  00002610 

GET  LIST  (SUBAT  (I)  .SPIALLK)  ; END;  ELSE  SHBAT  (T)  . SH ALLK*PACTCP * 03Q''262G 

SOBAT(I)  . SB  ALLK;  END;  IE  WORD  (4) -’OVPPI  DE*  TH’^N  GET  LIST  (CAPK)  ; O0CC263O 
ELSE  CAPK*BTGK(SBALLX,  BUNA’*’)  ; PUT  SKIP  EDIT  (•  BIGK*»  ,CAPK)  0)002640 

(COL  (2)  ,A,y  (1)  (9^31)  ; 0000265'' 

GO  TO  000C2660 

0')002670 

PPOC(IO):  /♦  GR?D  •/  00002680 


GF»*1*B;  CALL  PP0C7A;  GP**0*3;  GO  10  GZTCOB; 

/♦  PROCECOPE  TO  RESET  OPPSpT  LIST  •/ 

SFTOFF:  PRCC;  NOTrur;  nOFP*PTPST;  PNAHE*PNAHE;  CALL  RPFETOF; 
END  SETOPP; 

RESETOP;  PROC  RFCUPSt^i*. 

DCL  TFNP  PTF,T  PTXED; 

DO  1*1  TO  NDBAT; 

HUTN*NOTN»1;  ri'»’NABS  (NOTH)  =S^BAT  (I)  . ONAp  E; 

OITST(NOTN)  *BnrP:  NSDfl  (NOTN}=T; 

IF  SHBAT  (1)  .CHAlNF-^sSULL  “HEN  DO; 
tFPP*HUPP;  BUPP^SOBAT  (I) .CHAIPP; 

CALL  PPSFTCP; 
f fOPPING  n?  */ 

HOFPaTBBP;  IMF; 

END; 
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FETtlRM:  EMC  RFSFTCF: 

/•  PBOCEDUP*  TO  HAKE  S^’ST  •/ 

GFTBOIT:  PPOC  PFCU®SIVE; 

DCL  TF?1P  PTP,I  Fnro; 

11:  DO  1=1  TC  PUPAT; 

POT  SKIP  FDIT  (*  IMptr  SA'^E  AND  ®APr;^  FOR  ATTR  •,I,»CF  OTTL  FONC  *, 
BNABE)  (CCl(2).A,PO),T|1)  ,A.A)  ; 

CALL  GETLINB;  SOEAT ( I) . UHAFE^wCP^  (1)  ; R1=W0P9  (2) ; 

•2*iORD  (?) ; 

IF  P1=P2  THEN  ^0: 

POT  SKIP  LIST(*HOW  NANY  ATT®.  A® R TN  THIS  HOP?:*): 

GET  LIST  (NAT);  •^EHP=H»IFP; 

/•  CPEATE  A HOP  •/ 

ALLOCATE  PUP  IN  (PCOY)  : NnHAT^NAT; 

TEf1P->SH3AT(I)  .CMATNP  = NU?P: 

HRAHE=TP?1P->SUEAT  (I)  .»JNA“E:  CAPK=0; 

/*  PF.CUPSTVZ  CALL  •/ 

CALL  GETHOLT; 

/*  POPPING  rj?  AGAIN  A®T®R  HPC^IPSTON  •/ 
f!nPP*TEflP;  END; 

ELSE  CO;  SOP  AT  (I)  .CHA  IN®  = NTJ  LL;  ALLOCATE  ONT?  TN  (BOOT)  ; 

SOBAT  (I) .UNirT?  = HNlPP;  UTToE-0;  nL?=®1;  091=92; 

END; 

SPBAT  (I) .SHAllKs 1./N0®AT; 

END  LI; 

BETOPN; 

BHD  GETNOLT; 

/•  FRFE  F0R1AT  PFAD  CASD  RCOTTNE  */ 

GETIINE;  FPCC; 

DCL  (T,J,K)  FTXFD; 

DO  1*1  TC  10;  WOPC(I)=*  •;  END; 

GET  EDIT  (CARP)  (A  (PO)  ) ; 

1*1;  K*1; 

DO  HHILE(T<=0C) ; 

IF  SOBSTR  (CAPn.I  ,1)  =•  » THFN  GO  TO  CONT; 

J«1;  IF  I=PO  THEN  GO  TO  GOT; 

HORE:  I*  SOBSTR  (CARD,  !♦, 1,1)  =•  » THEN  GO  r,CT; 

J*J*1;  TP  (T^J=P1)  TRFfI  GO  TO  GOT;  GO  TO  HOPE; 

GOT;  HORD(K)  =SnBST?  (CARD,T,J)  ; K = K4-1  ; I=IfJ-1; 

COPT:  I=I^1;  END; 

BITORN;  ENC  GETT IN®; 

OHICAL;  PROC(T,ANS): 

/•  OHIFP  IS  ASSO®?0  POINTING  AT  TH®  PROP®®  D’^ITTT  ®nNCTION  */ 

/•  PROCEDUFE  TC  CALCOLATF  THE  OTTLTTY  OP  A VALD®  •/ 

DCL  J PIT®C; 

IF  OTTPE=C  ;H®N  PC;  ANS=X;  ?®TU®N;  END;  /•  LIHEA?  H(X)=X  */ 

/•  CCNSTANT  RISK  »/ 

TF  UTYPE=1  THEN  CO;  A NS =OTF  (1 ) • (1 -BXP  (-OYP  (1) ®T) ) ; ®PTO?N;  END; 

/•  PTFCEHISB  LTN®AF  •/ 

If  0TTPE=2  TR®N  CO; 

IF  I>«1  THEN  ANS=1.^  (HYP  (NUP) -UTP(NOP-I)  )/(OXP(NnP) -aXP(NnP-l))* 
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00003130 

0000.3140 

00003150 

OOC03160 

C0003170 

00003190 

0000319'' 

0000320'' 

(?9nn32i'> 

0000322*' 

00003730 

0090324'' 

00073250 

00003250 

0009377D 

000032P0 

0000329? 

00073.3CO 

00003319 

0000332'' 

090‘'333'' 

0*'0C3‘>40 

000973s: 

00003359 
O02C337" 
000 '13  3 09 
0 3:0.139  7 
OOCC3409 
00003419 
7 3C0  lu^** 
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(T-1.); 

ELSE  IF  X<*0  THFN  ANS=  !UTP  (2) -HYP  f1))/ |HXP  (2)-OXP  ( 1)  ) *X  ; 

ELSE  IF  (X<1  C X>C)  THE.’  TO; 

CO  J»1  TO  NTIP:  TF  fJXPCJ)>X  TH*H  GO  TC  CAL;  EMO; 

^Al; 

A9S>nTP(J)  ♦ (UYP  (.141) -IITP  (J)  )/ (OXP  (J^I)-OXP(J))  •{t-nXP(J))  ; 
EETOSN;  END; 

EID  ONI CAL; 

OHIFOtPPOC  (HALT, AMS); 

/•  PFOC  TO  CALCULATE  '‘XP.  UTIL  FOR  OMIAT.  UMIFP  ASSUMED  SET  •/ 
DCl  (J,NX,JU,JP,  KL)  FTXED; 

OCL  (SO(*)  ,S*(*)  ,TX  (•)  (*))  CONTPOLLFD; 

DCL  1 ALT,  2 MPA,  2 XPA(9),  2 CPA  (9),  2 FOA; 

ALT*PALT(NALT) ; ANS=0; 

IP  nTTPE*0  THFN  DC; 

DO  J**2  TO  MPA;  ANS=  AMS*  (CF A (.1)-CFA  (J- 1)  ) / |XPA  (J) -XFA  (J- 1)  ) * 
(IPA  (J)*XFA  i.7)-TFA  |J-1)  *XFA  (J-1))/2.; 

EMO;  RCTOFN:  END; 

ELS?  IP  UTT«F=1  THEN  DO; 

DO  J*2  TC  NFA;  AMS*  AN  S*  (CPA  ( J) -CPA  |J- 1)  ) / |XFA  (J)  - 7FA  *.1- 1)  ) • 

UXP  (1)  ♦ (XFA  (J)-XFA  (.1-1)  ♦ (EXP(-1¥0(1)  •XFA  (J)  ) -FXP  (->?YP  (1)  * 
XPA(J-1)))/nYP(1))  ; FND;  RFTORN;  ESP; 


00003430 

0''C0344n 

00003450 

00CC34SC 

0000347C 

00CO34P0 

00003490 

OOrC.3500 

00003510 

00''03S29 

0C003530 

00003540 

OOC03550 

00003560 

0*'00357O 

000035flC 

000035Q'' 

O3C036C0 

0000361'' 

00003620 

OOOC3639 

n0003640 

0OCO3650 

0000366'' 

00003670 


ELSE  IP  nTYPE=2  THEN  00; 

JP»1;  DO  JU-1  TO  MOP;  IF  (UXP  (JO)  >X PA  (.TF) ) THFM  GO  TO  ALOC?  SMD; 
/♦  INTEGRATE  CNLY  WHF’E  SCHB  P?OB.  TS  •/ 

ja»Mnp;  /•  THIS  LAST  IN  CASE  HE  »ALL  OUT  OP  LOOP  •/ 

ALOC:  Nl*NPA^N0P-.10;  ALLOCATE  SU ( ML)  , SF  (NL)  , XX  (NL)  ,B  (NL)  ; 

/♦  ML  IS  HAX  NOPEEP  0*  INTERVALS  PFOOIRFD  •/ 

; 

DO  J*1  TO  ML;  /•  PUT  INTERVALS  IN  ORDER  ♦/ 


TO  XNC.Tf;  IP  OXPI.Tn)<-rXFA  (JF) 
I^  UXP  (JU- 1)*XPA  (JF)  TH'N 


END; 


NI«NX4l;  IF  T?A(.1P)>1  THEN  GO 
IX(NX)»OXP(JO)  ; JU*.THf1; 

JP*.Tf4l;  END; 

ELSE  IMCJN;  no. 

XX  (NX)*XFA  (JF)  ; .TF=JP»1  ; 

IP  (JP>MFA)  THEN  GC  TO  LOEND; 

IP(n(MX)<15XX  (NX)>C)  THEN  DO; 

SO  (NX)*  (UYP(.TU)-UYr(.lU-l))  /(OXO(.TD)  -UXP  (JO-1))  ; 
H(MX)*OYP(JU)-Sn(NX)  *OXP  (JU)  ; »^ND; 

ELSE  If  XX(NX)<^0  THEN  DO; 

SO(NX)«  (OYP  (2)-UYP  (1) ) / (UXP  (2) -UXP  (1))  ; 

B(MX)*OYP  (2)-SU  (NX)  *UXP  (2)  ; 7MD; 

ELSE  IF  XX  (NX)  >-«1  THEN  DO; 

SO  (NX)*  (OTP  (NUF)-UYP  (NO  P-1 ) ) / (UX  P ( NO  P)  nXP(NOP-l))  ; 
B(MX)*OTP(NUF)  -SU  (NX)  •OXF(NUP)  ; END; 

SP  (MX)«  (CFA(JP)  -CPA(JF-1))/(XfACJP)-TNA(JF-1))  ; 

END; 

LPENO;  DO  J*2  TO  NX- 

AiS«ANS^SO  (J-  1)  *SF(.1-  1)  * (XX(.T)  •XX(J)  -XX  (J-1)  *XX  (.1-1)  )/2.; 
ANS*ANS^SF  (J-1)  *H  (J-1)*  (XX(J>-XX  (J-1))  ; 

IMT  (MflX*3) ) *K  (NX**2/2  ♦BX)  •/ 

END; 

FREE  SU,SF,XX,P:  RFTORN;  END; 

E*D  ONIEO; 


03C0368C 

03C''369O 
DCC037rn 

07C03710 
00003720 
03CC3730 
00003740 
0O003750 
THEN  DO;C0003760 
C300377C 
07003780 
00r03790 
0000 3H00 
0C003H17 
C0003H20 
0C003R37 
00003B40 
00003B50 
00073660 
00003R70 

OOOC3HHn 

00003H90 

0700390'' 

C700371C 

00OC3O27 

O0CO3Q3O 

O0C03«40 

0C003950 

00003760 

00C73977 

0770398'' 

00003970 
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BO^TEf:  PPOC(ANS);  NGH^Ds<';  CML  KTflTCAMANS)  ; PFTOPH  ;FNO  flULTFV; 
ROiTCAi:  «FOC(ANS)  rErn^sivK: 

OCL  TIH»  PTP»  I PITFD,  P PLOA'^*, 

OCL  TPHPn  (12)  , NGO(  12)  ; 

IP  CAPK=0  THEN  DO;  /*  IDDITIV?  POPH  •/ 

INS*0;  DC  1 = 1 TO  Nlin.AT;  NG  D=Nr,  1 ; NGP  (I)  =MG?  AD ; 

IF  SnBAT(I)  .CRATMF --NMI.L  THEN  TP  PPLAG  THFM 
P*S0BAT  (I)  .UNIP':^->PALT  (IPALT)  .EUO; 

ELSE  F*SUDA’f  ri)  . CN  TPT»J->CALT  (ICALT)  . EUC; 
ftSf  CO;  /•  NEED  TO  FVAL  A P.U?  •/ 

TEBP'  — np;  POFP=SOBAT(T)  .CHAIN’’:  CALL  NfJLTCAI.  (R)  ; 

/•  t'OP  lip  •/  nUFP=TENP;  FND; 

' .<SsANS4S0BAT(I)  .SflALLP*?;  END; 

1 1=1  TO  NIHAT;  GRAD  (NGP  (I)  ) =SnBAT  (I)  . SPAT.T.K;  END; 

TORN;  ®ND; 

I I DO;  /•  nOLT.  ?OR*1  */ 

AX  ’1.;  DC  T=1  TO  NUNAT;  NGOAOsNGPAC^I;  NG» ( I? =NGPA D; 

If  'FA*’’ (I)  .CHAINF=N»ITL  '”HFN  IF  PFLAG  TH’’N 
R*50w'  ''(I)  .nNirTP->PAT  T(IFALT)  . ED^;  ELSE 
R*SOBAT(I)  . DNTPTP-'CALT  (ICAL*”)  .*nC; 

ELSE  DC; 

TE«P*BD»P;  HUPF=SnpAT  (T).CHATNF;  CAIl  finLTCAL(R); 

H0PP*TENP;  END;  A MS -ANS • (1 +CAPK* SUBAT  (1) . SH ALLK*P) ; 

TEHPn(T)rR;  FND; 

DO  1*1  TO  NONAT;  GPA C (NG P (I) ) =AN S/ ( 1 ♦CAPE* SUB AT (I) . SR ALLK^ 
TEBFT)  (I) ) *Sn?AT(T)  o SRALLK;  END; 

. ANS«  (ANS-1 )/CAPK: 

BET04N;  END; 

END  flOLTCAL; 

GBTGRAD:  PPOC; 

OCL  TIHP  PTP,PAC’*OP  *10AT,I  FIXED, .1  PIIED; 

*IGRAD»0;  PACTCF='!;  CALL  SE7GPAD(PACTOP)  ; 

CO  1=1  TO  NUTN; 

J*I-1;  TEHP»OLIST(I)  ; TP  RUppsTF'!?  THEN  GO  TO  GC.PAD2;  END; 
GGRAD2: 

IP  -PPLAG  TH*N  POT  SFIP  LIST 

(•  AITRIB,OTIL.  GPAD  COUP.  AND  ATTR.  GPAO  CCfP.'); 

ELSE 

POT  SKIP  LISTCATTPin  T.  DTIL.  GPADIENT  CONPOHFNTS*)  ; . 

00  I«1  TO  NGPAD; 

IP  -«PPLAG  C OLIST  (.7*T)->CHAIV»  (NS?I9  (I)  ) *NULL 
THEN  00;  L NIPF=OLI ST (Jfl) -^ONt PTP  (N SOB  (I) ) ; 

CALL  ONICAL  (CALX  (TCAT.T)  ,P1)  ; CALL  T FTCAL  (CALX  (ICALT)  ♦.  0 1,  R2)  ; 
OERIV=  (P2-P  1)  /(CALX  (ICALT)  ♦.01-CALX  (TCALT)  ) ; 

. OERIV=DEri  V/ (UHI-OLO)  : ?flT  ’PIT  (UTNA-®’ (.T^I)  ,GFAO  (I)  , 
CaAO(I)*DE?IV)  (COL  (2)  ,A  ,X  (1 ) , P (Ar  3)  , X (2)  ,F  (1  0,  J)  ) ; 

END; 

ELSE  POT  EDIT  (OTNAH?  (.UI)  ,GPAD(I) ) 

(C0L(2),A,X(1)  ,P(B,3));  1^0;  pn?  SMV(2);  PFTTIRN; 

PNC  GFTGRAD; 

SElGRiiD:  PPOC(PACTOH)  PECUPSIVF: 

OCL  TEHP  FTP,  I 'IXFD,  FACT  FLOAT; 

00  1*1  to  NDNAT;  NGR AD=MGPAD^ 1 ; 


0000400^ 
moouol  3 
00004020 
00004010 
00004040 
00004050 
00004050 
0OC04070 
Oi0^40H'' 
00004090 
OOC041CO 
00C041 10 
00C04120 

00004  no 

00004140 
0O0O41S0 
00004160 
00CO4170 
00004  no 
0irC41Q'' 
00004200 
00004210 
00004270 
010^4231 
O0CO474D 
O0CO425C 
0'>0C4?6C 
0000427^' 
00004280 
00004240 
000^^4  **04 
00004310 
00004320 
00004334 

04004340 
00C043S4 
00004 lAG 
00004370 
00004380 
03004390 
00004400 
00004410 
00004420 
00004433 
00004440 
C30CU450 
03004460 
00004470 
03004480 
03004490 
0O004500 
00C04S10 
000C4520 
03004S34 
03004540 
00004550 
030^4560 
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If  SOBAT(I)  .CHAINP-^aHOLL  THEN  DO; 

TEHP«HUFP;  HOFP*SUPAT  in  .CHAINF;  •AC2*FACTOR*r.FAD  (NGR AP)  ; 

Cill  SFTGRAD  (FAC2)  ; J1JPP-TEHP;  FNO: 

CRAD(HGPAD)  *FAC»'OR*GRAn<NGPAD)  ; END; 

ICTORN; 
fID  SETGRAD; 

OilGET:  PPOC; 

no  I«1  TO  NHTN;  IF  UTNAN**  (I)  sWOPD  (2)  TH»N  GC  TO  HG'TP; 

END;  UGETC:  TUT  SKI?  lIST(*nsir  NOT  »OUND»):  GO  TO  GFTCCH; 
OCETB;  BUFP*OLTST  IT) ; TP  SUBAT  (NSHB  | T) ) .CH A T NF-= NULL 
THEN  GO  TO  OGETC;  0NIFP=  SURAT  {NSOB  (T) ).  UN TPTP ; FE70«»N; 

FBO  ONTG7T; 

PBOC(11):  Pt^OC?12);  /*  INOTF1  AND  INDTF2  •/ 

/•  EASED  ON  'OUATION  U (X1  ,T1 ) =U  |X2,Y2)  WHEN 
CE1.T1)  IS  INDIF»»EPEMT  TO  CX2.Y2)  •/ 

CALL  GE‘»’2;  GO  TO  PPCC‘»iC; 

GET 2!  PPCC; 

/•  ROOTING  GPTS  2 n*«S  AND  S»^TS  0 1 , J2,K  1,  K 2,  UF,  HUFP 
JARTC  ).  AND  T */ 

JA8T  (1) ,JAPY (2) DO  T*1  TO  HUTN;  DC  J»1  TO  2; 

IF  »OPD(n.l)*nTNANE  (I)  THEN  DO;  JART(J)*I; 

60  TC  GET2B;  END;  END; 

6ET2Q:  IP  (JA‘»Y|1)>r  r.  JAPY(2)>C)  then  go  to  GET2C;  END; 

POT  SKIP  LIST  (•  ATT?.  NOT  BOTH  POUND'*);  GO  TO  GFTCCN; 
GFT2C:J1«JAPY  (1)  ;.12*JARY  12)  ; TF  (OT.TSTI.TI) -***CT  IST|J2)  ) THEN  DO; 

POT  SKIP  LISTI'ATT?.  NOT  IN  SABF  NUP») : GO  TO  GBTCOB;  END; 

. HnfP*OLTST (01) ; 

R1«BSnB(.T1)  ;K7=NSUP(J2)  ; TF  (CRAINF  |K1 ) -.*Nni,l  1 CHATNF  (K2)-*NnLL) 
THEN  CO;  0?*'C»R;  PUT  SKIP  LIST(»INFUT  UTILITY  VALUES*) : 

END;  ELSP  UF**1*B;  FE'^’UPN;  END  GfT2; 

fR0C11C:RrFP=0LIST(T)  ; R-SUPAT  (NSUB  (J2)  ) . SBALLK/SBALLK  (K  V)  ; 

IF  iOPOl  1)  **INDTP2*  THEN  GO  TO  poOClIF; 

POT  SKIP  LTSTI'INPUT  AN  IN0IPFEP2NCE  PAIP  PLEAS**); 

GET  IIST  (IXIN(I)  no  1*1  TO  4));  IF  -UF  THEN  GO  TO  ncALI; 

DC  J*1  TO  2;  0NIPP*SU3AT  (NSUB  (JAPY  (.1)  ) ) .UNTPTR;  DO  T*JrJ^2; 
IIB(I)*  (XIK  (T)-ULC)/(nHI-ULO)  ; CALL  UNICAL  (XIN  (I)  , ANS)  ; 
IIB(1L)«ANS;  END;  END; 

OCAll;  OESisXIH  (1)  ‘XIN  (2)  -XIM  (1)  pyTH  (4)  ; 

TF  0ES1«0  THEN  GO  TO  PROClID; 

1ATK«  (R*  {XIN  (4)  -TIN  (2) ) ♦XIN(3)  -XIN(1)  ) / (PPPPSl)  ; 

IF  RATK«0  THEN  DO;  EK=C;  GO  70  PR0C11**;  END; 

BK*1;  DO  T«1  TO  NUBAT;  BK=BK*  (U 0ATK*SU3AT|I)  , SHAILK/ 

SOBAT  (NSOB  (.11)  ) , SBALLK)  ; END;  BK  = BK-1; 

FAC1«P ATK/ (SNALLK (K1) •PF) ; 

FAC2«=- (NBA  UK  (K1)  ♦SBALLK  (K2)  ) /(BK»SBALI  K (K1)*SNALLK  (K2))  -FACI  ; 
POT  SKIP  EDIT  (*  INPLIED  NFN  K**S  FACTOP(S)  * , *AC  1 , * ( * , *AC2r  *)  *) 
(C0L(2),A,F(R,J),X(1).A,F(«,1)  ,A)  ; 

FBC-11F: 

POT  SKIP  EDIT(*IN?LIFD  NFN  3IGK*  *,BK)  (COL  |2) r A , F (B ,1) ) ; 

60  TO  GFTCCB; 

FB0C11D:  POT  SKI®  LIST(*INDIF  PAIR  YIELCS  INFO  ABOUT  PEL  K**S*); 

POT  SKIP  EDIT(*REL  K CR*C*r.  ClIFRENT  RATIO  •,  N0PD(3),*  TO  \ 


tORD(2) 


r®)  (COL(2),A,A,A,A,A.F(H,3))  ; 


f(IIN(3)  -XIN  (1))/(X  .N(2)-XIN(4)); 
POT  SKIP  FD1T(*IBPLIE0  PATIO  * 


(COL(2),A,F(S,3)); 


03C04S70 

0*^004580 

03004D9C 

00C04f0C 

0*'004610 

o:'C04ft?D 

0000463'^ 

0300464'^ 

CD0C4fS0 

00n04#^f»3 

00C046‘*O 

C00046RO 

030046Dn 

00004700 
00C04710 
00004720 
0O0C473C 
00C0474D 
0D0047«i0 
00004760 
00004770 
00004780 
0OC04790 
0^00480'^ 
C0004810 
O0OC4820 
0'>004h'»D 
000C48U0 
003048*50 
OOC^'JRfO 
00004870 
0000488'^ 
0^004890 
00004900 
00004910 
00004920 
0000493v0 
O0CO494C 
0000495*' 
00004960 
00004970 
00004980 
00004990 
00005000 
O0OOSO1O 
07005C20 
0 )005030 
010*'5C40 
00005050 
00005 060 
00005070 
00005080 
0O00509*' 
000051C0 
000C5110 
000C5120 
00005130 
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GO  TO  GETCOP: 

PROC11E:  POT  SKIP  LTST(*IVpnT  2 IHnTFPERFNCF.  PAIRS  PiFAS?*)  ; 

GET  IIST((XTN(T)  DC  1=1  TO  fl)); 

IP  -»0P  THFN  GO  TO  DCAL2: 

DO  J=1  TO  2;  UNIPF  = SUBAT  (MSnn  (JAPT  (.1)  ) ) .riNIPTF; 

DO 

IIM  (i)=  (iTK  (T)-nLO)/ inni -fjLO) : call  onical  (xin  (I)  , a'js)  : 

XTH(I)=ANS;  END:  FNO; 

/♦  CHICK  BOTH  DESCPIMINANTS  F^'P  P''D0MDANCY  ■/ 

DCAL2:OES1=XTN  ( 1)  *XIN  (2)  -XT»M  1)  •XTM  (U)  *. 

DES2*XIN  (5)*XIN  (6  ) -XIN  (7)  •Xi  . 

IF  (PES1=0| OES2=0)  THFM  DO; 

POT  SKIP  EDIT(*0'IC  TNDI®  PATH  vDpVTOUSLY  YI*‘LDS  ?"L  K»*S.  OS^*, 

• THE  CCftPAND  IN0IP1  WITH  IT  TO  CCHFA*»F  WITH  cnPRFNT  »ATTO*  ) 
(COL  (2)  , A,  A)  : GO  TO  GETCO**:  «?Nn; 

Q1*XIN(7)-XIN(S)  -(XIN  (3)  -XIV  (in  -CES2/DES1; 

02«XIN  (6)-XTN  (a)^  (XIN  |U)-XTV  (2))  *DES2/nES1: 

IF  (01  = 0|02  = C)  *^HFN  DC; 

POT  SKIP  LIST  ( 'CANNOT  D''TF?HINF  PFL  K'»S  PRCH  TH  ES PTS.')* 

GO  TO  GETCC-;  END; 

PAC1=  (Q1/02)*  (TIN  (U)-x:n  (?)  ) ♦XTV  (31-XTN  (1)  ; 

FAC1*PAC1/  ( (01/02) 'CES1)  : 

POT  SKIP  FDIT('9IGK-',?AC1,*/K  (• ^WOPO(2) »•) •) 

(COl  (2)  ,A  ,F  (8,3)  ,A,A,A)  ; 

XIN(3)*01;  XIN(2)=02;  XIN  (4)  ,x*IM  (1)  =0:  GOTO  PR0C11D; 

PROC  (13) : /•  UNTCAl  •/ 

CALL  OMIGET;  IF  W0?D(1)='  • THFN  GC  TO  PR0C13C; 

. N*H0PD(3);  GET  LIDT((XIN(I)  DO  1=1  TO  N)); 

PR0C13B:  rO  1=1  TC  N;  P=  (XIN  (I)-MLO)/(nHT-ULO)  ; 

CALL  UNICAHP,ANS)  ; FUT  EDITC  U(' , XIN  (I)  , • )=  • ,ANS) 

(COL  (2), A, F (10,3)  ,A,F  (9,3))  ; 
mO;  GO  TC  GETCC; 

PR0C13C:  XIN(1)=0L0;  DO  1 = 2 TO  1C  RT  2;  XIN  (1^I/2) =010X 1 . F-1*I* 

(OHI-OIO)  ; FNF;  N=6 ; GOTO  ^P0C1.3H; 

PROC  (14):  /•  • INVERSE'  •/ 

CALI  ONTGBT; 

IP  WORD  (3)='  ' TH'-’H  GC  DR0C14r; 

K*WORD(  j)  ; 

GET  LIST  ((XTN(I)  DO  T=  1 TO  N) ) ; 

PR0C14B:  CO  1=1  TO  N;  CALL  (JNIN7  (XIN  (I)  , ANS)  ; ANS  = ULO  ♦ (HHI-ULO) 

♦ ANS;  PUT  ECIT(.ANS,'  = INV  (' ,TIN  (I)  ,' ) •)  (COL  ( 2)  , P ( 10,  ^)  , A, 
F(«»3)fA):  END;  GC  1C  G*”^CDN; 

PR0C14C:  XTN(1)=0;  DO  1=1  TO  9 PT  2; 

Xiri(2^I/2)  =1.F-1*I;  END;  XIN(7)  = 1;  N=7;  GO  TO  PR0C14B; 

OHIMV;  PROC  (T,ANS)  ; 

/•  PROC  TO  GP*^  INVERSE  OR  ANS=X|  H(X)=T  •/ 

DCL  J FIXED; 

IF  OTTPE*0  THEN  DO;  ANS=Y;  ^FTORN;  END; 

IF  0TYPE=1  THEN  DC; 

AKS*L0G(nXP(1) /(GXP(1) -T) ) /0YP(1)  ; FETORN;  END; 

IF  0TTPP»2  THEN  DO; 

IF  (T>1(Y<0)  IHEM  HESS:  DO; 

POT  SKIP  FDIT(Y,'  IS  CUT  OP  RANGE')  (COL  (2)  , P (9,  3)  , A)  ; 


00005140 

oocnsiso 

00005160 

00005170 

000051R0 

0300519r> 

00005230 

O0CO521O 

0*'00522'^ 

03005230 

03005243 

030rs25o 

030DS?f 

0''005270 

0''CC5293 

03005290 

03005303 

00CD5310 

0DC05‘»2O 

0300533'' 

00005343 

0300535C 

03035360 

03005370 

0'rC053R0 

00C053OC 

330054C3 

03005413 

0)005423 

000C543O 

0O0C544C 

030C5453 

0300546O 
C"-:  054  70 
030^5483 
0300549O 
0300*^5C'' 
00005510 
00005523 
030^553C 
OOOOSSf'i 
0-»0:5553 
G000556C 
03005*^73 
H3CC55B0 

03035590 

0300560'' 

03CC561'' 

03005620 

03C05630 

03005643 

0'*005650 

00035663 

00005670 

03005600 

03005693 

030057CO 
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RIIS*-1;  RPT0?»N; 

00  J«2  TO  HO**;  I»  (UYP(.T-1)-r)*(nY*»(.T)-V)<«0  TriEtl 
GC  TO  COT;  FHO;  00  TO  H?SS; 

GOT:  INS-  (Y-OTP  (J-1))/(UYP(J)-nYP(.T-1))*(MXF  (J)-UYP(,l-1)) 

♦ OXF(J-I);  RTTOPN;  ENT;  END  nHiHV; 

PHOC(IS)!  /•  'CHINGOU.T*  •/ 

PFUG*«0»D:  AP«*0*n.'  nr.  1*1  TO  10; 

IP  HinCILT  (I)  sWO’OC  3)  then  nc;  ICILT^I;  GO  TO  CALLILT; 

END;  FHD;  PFLAG-M*3;  DO  T*1  TO  S;  IP  VAH?nT  |X)  =WOPD  ( 1) 

THf.n  DO;  TPALT=I:  GO  TO  CALULT;  ENt;  F?ID; 

POT  SKIP  LIST  (*  ALT^PM.  NCT  POUND*)*.  GO  TO  GETCO-*; 

CULALT:  CALL  UKIGET;  /♦  SPTS  UNIP?  ASC  I •/ 

CAil  ALTCOHP;  GO  TO  GETCCH; 

PROC(lfi)  : /•  change  •/ 

IF  HORO  (2)*PHAn^  THEN  GO  TO  PR0C16C; 

DO  I«1  TO  NOTH;  IP  OTNAnF  (I) -NORD  (2)  THEN  GO  TO  poOClKB;  END; 
POT  SKIP  LTST(*ATTPID  M0‘»*  •OUMD*);  GC  TO  GPTCOM; 

PR0C16B:  BnPP*OlIST(T) : IF  cnATNE(NSUR(T) )=NOLL  THEN  PO; 
OHI»P*nNIPTS  (NSUO (I) ) ; UP-*1*B;  END;  ELSF  UP**0*B; 

IF  HCRD  (3)  **  NAME*  THEN  DO; 

OTHAHI  (I)  , nNAM'’ (NSns  (T)  ) =NORD  (4)  ; IF  -.U** 

CHAINF  (HSUB  (I)) ->MNAf1E=B0*>D(4)  : END; 

ELSE  IF  H0PD(.1)**K*  THEM  DO; 

SRA11K(NSUB  (D) -W0?L(4)  ; CAPN*BIGK  (SMALLK, NOHAT)  ; END; 

ELSE  IF  (HORD  (3)=»FANGE*  F,  HP)  THEN  DO; 

POT  SKIP  LIST(*ALTP*NATIVE  C0M*>0N?M1S  NEED  CHANGING*); 
rOT  SKIP  LIST(*RANGE  PLEASE:*);  GET  LIST (UL0,UHI)  ; 

• ■ORO(3)«*  *;  GO  TO  FPOCHC;  END;  GO  10  GPTCON; 

PR0C16C:  HO»psPIPST:  IP  NORD (1) ** NAME*  THPN  PNAHE, HN AHF»iOPD ( 4) ; 

GO  TO  GETCCH; 

PS0C(17):  /*  ALTIIST  */ 

AF«*0*B; 

J1«0;  00  I»1  TO  1C;  IF  NAHCALT (I) -»•  ' THEN  DO; 

1F«*1*B; 

PCT  »DIT  (NAMCAL'*' Cl) ) (COL  ( 12*J  1)  , A)  ; END;  END; 

IF  -»AF  THPN  DO;  PUT  SKIP  IIST(*N0  CERT.  ALTEPN.*); 

GO  TO  ITSTPPB;  END; 

POT  SKIP(2)  ; DO  I»1  TC  NOTH;  NUF*>-OLT ST  (I)  ; IF  CHAINF  (NSUH  (I)  ) 

■ NOLL  THEN  00;  UN  IFP-UN  IPTR  (NSUR  (I) ) ; PUT  FDIT  (UTMA'*E  (I)  ) 
(COL(?)  ,A  (1C))  ; J1-0;  00  J*1  TO  10;  IP  NAHCALT  (J) -.*«  « '»*HPN  DO 

J1-JU1;  X*nLOtCALX  (.1)  • (UHI-ULO)  ; PUT  ^OIT  (X) 

(C0L(12*.11)  (10,3))  ; END;  END;  END;  FNO; 

LISTPBB:  A*»*d*B; 

POT  SKIP  (2)  lIST(*CE?i'  EQUIV.  TABLE  FOP  P»OB  AT.TEPN*); 

J1*0;  00  I«1  TO  5;  IP  HA HPALT (I) -»*  * THEN  DO; 

AF»«1«B; 

J1-JU1;  PHT  PDTTCNAHPAIT  (I)  ) (COL  ( 12M  1)  , A)  ; FNO;  -NO; 

IF  -»AF  THEN  DC;  OUT  SKIP  LIST(*MO  PFDB.  ALTFRN.*); 

GO  TO  GETCCH;  '’ND; 

POT  SRIP(2;  ; PO  1*1  TC  NUTM;  HUPP-OT.IST  (I)  ; IK  CHAINF  (NSUH  (T)) 

■ Hnil  THPN  DO;  tINI*p*nNIPTR  (NSnp  (T)  ) ; PUT  EDIT  (nTNAHF(I)) 
(C0L(2)  ,A  (1  C)  ) ; .11*C;  DO  J»1  TO  5;  IF  NAPPALT  (J)-^**  * THEN  DO; 
J1»JU1;  CALL  UNTNV  (EUP  (J)  ,X)  ; X*ULO»X*  (OHT-ULO)  ; 

POT  EDIT(X)  (COL  (12*J1)  ,P(10,  D)  ; knD;  ?ND;  END*.  •ND; 


OOCOSTn 

0K0S720 

OOOCS73D 

0700S74U 

0')CC57S0 

OOOOS7f3 

00OCS770 

000US7HO 

0000S7U0 

C??C5RC1 

o*»oosRn 
0)0C5R20 
OUOOSH3n 
OOOOSH40 
OUOOSRSO 
07CC5BR'' 
oocc')H7r 
COCFSflH'' 
OlODSRur. 
01COSUO? 
0000501'* 
00005923 
0000593'* 
00005940 
OOC05950 
000059^;c 
00005970 
030059RO 
0000599C 
00035C30 
0*035010 
0O0U6U20 
00006030 
0DCC6*'4n 
00006950 
000C606'* 
0000607? 
900060R0 
0000609? 
0000610^ 
03C961 10 
0300612? 
;0)C?6  130 
00006140 
0300615O 

00O06160 

C000617? 

0C0061R0 

03036190 

00006200 

00005210 

0O0C6270 

C3006310 

3*ro6243 

03C0625? 

00006760 

O0C05270 
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CO  TO  CfTTCP; 

PR0C(18):  /•  niSPLAT  •/ 

IF  fOPO(2)  =FN^riF.  THRN  DO;  ’»fIPP-7T^ST;  CO  TO  npociflC;  FMO; 

00  1=1  TO  NUTN;  T®  U"'NA*1^  d)  =tf  OPO  (2)  THPN  CO  TO  PP0C1PP;  *'VD: 
PUT  SFI»’  TIST(*ATT"I9  NOT  ^OUMD*):  GO  TO  RFTC0»1  ; 

FPCC18B:  8nPP=0lIST  (T)  ; IF  CHAIN*  (NCUD  |T)  ) =NI?LL  TH*M  DO; 

OMI?P=aNT?T'' (NSliri  m ) ; PUT  SKIP  EDIT  (•  T ANC*:  • ,ML?,U”r) 

(COM2)  ,A,  (2)  P (1'' ; I?  UTYP*=0  TH*N  Pn?  SKIP  LIST 
(•UTTPS  IS  ITN?AP»)  ; 

USE  IP  nTYP-  = 1 TH*N  DO; 

PUT  SKIP  LTST(MITYP?  iS  COMSTAN'"  U (X)  =P  f 1 -*XP  (-CX)  ) • ) : 

POT  TDIT(* P=* ,OXP (1) ,*C  = * .HYP  (1)  ,*VAPIA  ILF  NORNALIZFO*) 

(COM2)  .A.*(fl,3)  ,X  (1)  ,A.P(fl,3)  ,X(2)  , A)  ; 

IF  TIXP(1)>0  THEN  AN\*-  = * f VERSE*  ; *ISF  A NA  HE  = * PPCNE* ; 

POT  EDIT  I *PISK*,  ANANF)  (C0L(2)  , A,  X (1 ) , A)  ; 

IMP; 

flSE  I*  UTYF*  = 2 THFN  P'^; 

POT  SKIP  lIST(*nTYPE  IS  PI*f!EWlSE  LTNPAP*); 

00  J=1  TO  NftP;  IF  UHKUL^  TH  *N  CO;  X = UXP  (Srj?-J4l ) ; 

T«OTP  (NT1P-J41)  ; *Nn;  ELSE  DO;  X = IIXr|J);  T = UY?(.T);  END; 

X»nL04X*  (UHI-ULO)  ; FMT  "Cl T ( • U { • , X, • ) = • , Y) 

(COM2)  , A,P  (10, 3)  ,A  ,F  (P  ,1) ) ; *»ND;  END;  CO  TO  GETCOH;  **11; 
KOFP=CHATNF  (fJSlJP  (I)  ) ; 

PROCIflC:  POT  SKIP  LIST  (*ITSTING  OF  T *ACTORS»);  SUflK^C; 

00  .1=1  TO  NCMAT;  SflRKsSUFK^S-ALL*  (J)  ; 

POT  E0IT(UNA*1E  (.1)  ,S**ALLK  (.!)  ) (COL  (2 ) ,A  , * (8 , 1)  ) ; END; 

POT  SKIP  FDIT  (**'IGKs«  ,CA«'K,*  SUN  K*»S  =*,SU«K) 
(COL(2),A,F(8.1),xn),A.P(8,l));  GO  TO  GFTCTH; 

PPOC  (19) ; /•  PPACTTL*  4/ 

CALL  UMIGET;  DC  J=1  TC  5;  IF  NA  NPA LT  (J)  =WOPT>  ( 3)  THEN 
GC  TO  eP0C19N;  FND; 

POT  SKIP  LIST  (• alTEPN.  voT  PrUND*);  GO  TO  G*TCOH; 

PF0C158:  POT  SKI?  l:st(*cun  distp!9  pop  the  ALTEPN.*); 

DO  1*1  TO  NP(J);  IF  IJHKOLO  TH*N  DC;  X=  PALT  (.1)  . X F (Nf  (J) -!♦  1)  ; 
T*1-PALT(J)  .CP(N*(.1) -I4l)  ; END;  FLSI  DO; 

X*FALT(.1)  .XP  (I)  ; Y=  OA IT  (.1)  .CF  (I)  ; END; 

I«0LO4X*  (UHI-OLO)  ; FMT  EDI  T { • F ( • , X,  • ) =•  , Y) 

(COM2)  ,A,P  (n,3)  ,A,F  (8.3)  ) ; END;  GC  TO  GETCCH; 

PFOC  (20) ; /•  LOTTERY  •/ 

CALL  UNTG*T; 

IF  HOPD(l) -•=*CE*  IWOFD  (3) -i=*P*  THEN  10; 

N-NOPC  (3)  ; I*  N>1  THEN  DO; 

POT  SKIP  IIST(  *T,CTT?«Y  FNDFTS.  PLEASE?*); 

GET  LIST  ( (XIN  (.1)  DO  ,1  = 1 TO  N) ) ; 

POT  SKIP  LIST(*rO^PESP.  PPO 3A0TLTTI ES  PLFASE?*); 

GET  LIST  ( (YIN  (,1)  DO  .7  = 1 TO  N) ) ; X = 0 ; DO  J*1  TO  W; 

Xm  (J)=  (TIN  (.7)-UL0) /(U»I-nLO)  ; CALL  ONICAL  (X  IN  (.7)  , ANS)  ; 
X«X4AKS*YIN  (J)  ; FND;  GO  TO  PP0C2CP;  ;ND;  END; 

POT  SKIP  LIST(*IN?UT  L'.'TTEFY  FNDPTS.  (POTTON.TOP)  AND  THE  CF  OP  P*) 
GET  LIST  ( (XIN  (,1)  DC  .1  = 1 TO  *») ) ; 

DO  J«1  TO  2;  XIN(J)  = (XIN(J) -IILO) /(DHI-OLO)  ; 

CALL  ONICAL  (XIN  (.1)  , ANS)  ; TIN(J)=ANS;  END; 

IF  fOPD(3)  =*CE*  THEN  CO; 

X*XIN  (3)*XIN  (2)  ♦ (1-XlN  (3) ) *XIN  (1)  ; FPOC20B:  CALL  ONTNT  (X,  AN  S)  ; 


07008280 

00008280 

03''C830C 

03C0A3ir' 

03078328 

80008318 
f»088f,  lu  n 
C8CC83*^C 
03C0636C 

0C008378 

08708?98 

33078  3Q-' 
08C76U  )0 
03CO8U 10 
O^COGUE^' 

00008t{  10 

0''006U'iC 

03C08US0 

08C06U83 

00C08U70 

ODOOfiunO 

770CEUQ8 

08008588 

07008510 

07076527 

OOOOf 510 

OrO('f  5U0 
700C85S'' 
00086580 

C700857O 

07C  385P'' 

07706590 
00006800 
07008810 
0708f  8?7 
000868^7 
070C88U0 
000C6657 

09008680 

080868*70 

00088668 

08086807 

07C067C7 

008C8718 

07078720 

OOOC8738 

0'‘0068a7 

070067SO 

00008780 
07CC8770 
08C867PO 
; 01088707 
00C06PC8 
07078610 
07086628 
C8C08P30 
03068698 
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X»OIO^ANS*  (OHI- niO)  ; PtIT  SKIP  EDIT  (*CE  FOP  LOTTEPY*  *,X) 
|COI(2),A,F(10,^))  : ?ND; 

Ei.^F  If  BOPPn)s*P*  TR'‘»I  DO; 

Xm(3)  »(XTN(3)-ULO) /(flHI-OLO)  ; CALL  HNICAL  (XIN  (3)  ,AMS)  ; 
X*(AKS-XIN  (1))/(XIN  (2)-XIS  (1))  ; PUT  SKIP  FDIT 
(•P  FCP  lOTTfPY-  *,X)  (COL  (2),A,P(fi,3))  ; END; 

GO  TO  GETCCH; 

PPOC  (21)  : /•  IFAD  •/ 

CALL  GET2;  snFP*Ol  1ST  (I)  ; PUT  SKIP  LIST 

(•IMOUT  TNOIF  FT.  THFCUGH  MHTCR  CDH3!  WHL  PASS:*)*. 

GET  LIST((XIN(.l)  00  J-1  TO  ?)); 

IE  -^UF  TH»^N  GC  *^0  PPCC21P;  DO  J=1  TO  2; 

OlfIFP*SnPAT  (NSUP  (JAR7  (J)  ) ) .UNIPTF; 

Xlir(J)  a(XIN(J)-nLO) /(nni-ULO)  ; call  ONICAL  (XIN  (.1)  ,AHS)  ; 
XTH(J)*^MS:  '•ND; 

PP0C21B:  X=SFA  LLK  (K1 ) •XIN  (1) ♦SFALLK (K2) *TTS (3) ♦ 

CAPK*STALIK  (K1)  *SS»I.LK  (K2)  *XTM  (1)*XTM  (2)  ; 

PUT  SKIP  LIST  ('INPUT  NU-DER  OP  PTS.  *0P  HAP:*); 

G''T  LIST(N)  ; PUT  SKT»  ^'OI''’ (•  INP’.T  •,UNAF'(K1), 

• fALOES  PCP  WAP*)  (COL  (2)  rA,A,A)  ; 

GET  LIST  ((YIM(J)  DO  J*  1 T>0  N) ) ; 

IP  OP  TH^'N  CC;  UNI  FF^nNIPT’ (K1)  ; DO  .1*1  TO  M; 

XIK(J)^  (YIM  (.D-tlT.n)/  (UHI-ULO)  ; 

CALL  OHICAL  (XIN  (.1)  , ANG)  ; XIN(J)*ANS;  END;  END; 

ELSE  DO  J=1  TO  N;  XIN  (0 ) * YiN  (J)  ; ENT; 

DC  .7*1  TO  M; 

XIU(J)*(X-SRUIK  (K1)*XIN  (J) ) /(CAPK*SFALLK  (K1)  *S  NALLK  (K2) 

, *X\H  (J)  ♦SFAII  K (K2) ) ; IF  UF  THEN  DO; 

0PIFP*0NIPTR(K2)  : CALL  UNINV (XIN  ( J) » AM S) ; 
xm  (J)*OL04ANS*  (UHT-fflO)  ; END;  »ND; 

VOT  SKIP  LIST  ( 'INDIFFEofNCF  PTS*); 

POT  SKIP  FOIT  ( (»  (*  ,YIN  (J)  ,»  ,Xl  i (J)  ,» ) • DO  .1*1  TO  H) ) 

(COL(2),A,?(1C,3),A,f(10,3),A); 

POT  SKIP  EDIT (* UTIL  FOR  CORVF  WITH  CTHER  ATTR.  AT  O', 

X)  (COL(2),A,X(3),F(fl,3)); 

GO  TO  GETCON; 

PEOC  (22)  ; /<*  STOP  */ 

POT  SKIP  LIST  ('THANKS  FOR  USING  NUFCAP') ; STOP; 

PEOC  (23):  /•  DELU  •/ 

CALI  DELUT(W0ED(2))  ; CALL  SETOFF; 

POT  SKIP  LIST('K"S  NEED  NOPNALIZTNG  AND  RIGK  NEEDS  SETTING*); 
GO  TO  GETCCN; 


DILOT:  PROC  (TNA"E)  ; 

OCL  TNARF  CHAF(12),I  FIXED,  IS  FTTEt; 

00  1*1  TO  Nirs;  IS=M5UB(I);  IF  UTNA P®*  (T)  *TNARE  THEN  DO; 
RUFP*OLIST(I)  : GO  TO  »OUND;  END;  ENE;  POT  SKIP  EDIT 
(THABE,'  NOT  IN  MS^')  (COL  (2)  , A , A ) ; GO  TO  GETCCH  ; 

POOMD;  IF  (NUFAT-1)=0  THEN  DO;  PUT  SKIP  LIST 

(•PLEASE  DFLi'TE  TRF  NU»  Tn  WHICH  THIS  ATTR.  BELONGS'); 
BBTUR'.;  •’ND; 

POHAT*MnBAT-l;  DO  T*  1 TO  NOBAT;  IF(I>*IS)  THEN  SUBAT(I)* 
SOBAT(I^I);  END;  PETOFN; 

FiO  DBLOT; 


112 


O.R20225,1194C'.f1UPCHP.PLI 


PROC(24):  /•  Horn  ♦/ 

CALL  ADDUT(«G?D(?)  ,«TF!)(3)  ) ; C\I.L  SETOFF; 

POT  SKIP  LTS*^  (MLTERN.  CO«P.  "AY  S^TTINO*); 

PUT  SKIP  LIST(»K**S  NOOHAM ZING  AND  BTGK  NEEDS  SETTING*); 

GO  TO  GETCCN; 

ADDOT:  rPOC(TN1,TN2)  ; 

D:T,  (TN1,TN2)  CHAP(12),T  '•IT-O.rs  FTTFD.’^ENP  P'*’E; 

TF  TH2~FNAHE  TH*'N  DC;  NOPp=FTPST;  GC  TO  FOUND;  ’^NP; 

DO  T=1  TO  NIITN;  IP  (IT NA  N'' (T)  =“^  v2r.CI  T ST  ( I)  - >CM \I  (N 3MB  (1 ) ) -= 
null  TH»N  DC;  FMFP=CLIST (I) ->CHAINP  (NSUB  (I) ) ; GO  TO  POUND;  END 

PND; 

POT  SKIP  FniT(TN?,*  NOT  A HI*'*)  (COL  (2)  , A , A)  ; GO  TO  GPTCON; 
FCOHO:  NUN  AT=NDNAT^1  ; IS-MnNA-i  ; UNA  HE  (NM  PAT)  =TM  1 ; 

SNAT.lK(KUnAT)  =1. '•O/NU-AT;  ^MT  SKI®  mT( 

•INPUT  PANGE  FO?  ATT^.  '^NUrAT,*  OF  MTIL  PmnCTION  *,TN2) 

(CCL(2)  ,A,F(2)  ,A,A)  ; GET  LIST(M,32); 

IP  R1  = ?2  THEN  DO;PUT  SKIP  LIST(»HOW  PANY  ATTE.  IN  THIS  HUP?:*) 
GET  IIST  (NAT);  T'‘NP  = '*m ?p;  ALLOCATE  NM»  IN  (ROPY); 

NUBAT=NAT;  TENP->CWArNP  (IS) =-M^P;  NNANF=TEN0->UNANE (IS)  ; 
CAPK»C;  CALL  GFT-MLT;  1UFF=TFir; 

fLSE  DO;  CHAINS  (NUPAT)=VUT.L;  ALLOCATE  UNIP  TN  (PODY)  ; 

ONTPTR  (NiJNAT)  =nNTFP;  MTYPE  = :;  ULO-<»1;  UHI  = P2;  7NO; 

2BTUBN;  END  ADDJIT; 

PBOC  (25) ; /•  SNITCH  */ 

* CALL  SHITCRO  (SORD  (2)  ,HCRO  (3)  ) ; ..’ALL  S^TOPP; 

POT  SKIP  LIST(»K**S  IN  B^TH  HUPS  NFft  NOfH AL I ZING* ) ; 

GO  TO  GETCOH; 


SVB: 


SVC: 


SND: 

SVE: 


2; 


G’TCON; 

(J2)  ; 


SfITCHO;  PP0C(TN1,TN2)  ; 

DCL  (TN1,TN2)  CHAO  (12), I PIT^O,  TP^p  PTP; 

.7ART  (1)  ,.7APY  (2)  =0;  DO  1=1  TO  NUTN;  tO  .1=1  TO 
IF  VORD(1  ♦.!)  =0TNA'*'' (I)  THEN  00;  JARY(J)=I; 

GO  TC  SVR;  FND;  END; 

If  (JARY  (1  )>0F.JAPY  (2)  >'')  THEN  GO  TO  S?C;  END; 

IP  (JART(1)>0  GTN2  = FNANE)  THEN  GO  TO  SWD; 

POT  SKIP  LIST(»ATTP.  NOT  BOTH  POUND*);  GO  TC 
J1=JAPT  (1)  ; J2=JA?Y  (2)  ;K1  = NSUB(J1)  ; Hnpr=OLIS- 
TF  CHAINF  (N’^UR  (J2)  ) =n»!LL  TmEN  DO; 

POT  SKIP  EDIT(TN2,*  IS  NOT  A NUP*)  (COL  (2) , A , A ) ; 

GO  TO  GETCCH;  END;  f*UPP=CHA  I (NSUB  |J2)  ) ; 

TFflP=CLIST  (.11)  ; GC  TC  EVF; 

J1*JAPY(1);  K1  = NSfJ?  (.11)  : EnFP  = ?IoST;  TENP  = CLT3T  (.11)  ; 
HUNAT=NT*AT^1  ; UNA.NE  (NUflAT)  =TN1 ; SNALLK  (NOSAT)  = 
KPO/NUHAT; 

IF  TENP->CHATN?  ;K1) =NULL  THEN  DO; 

OHTPTP  (NUHAT)  =TF-?->UN  I PTP  (K  1)  ; CHATNP  (NUNAT)  =NULL;  EN’*-; 
FLSE  CHAINF  (NUflAT)  =TErP-X:HAINF(Fl)  ; CALL  D-!LUT(TN1); 
RETURN;  END  SWI^CRU; 


PRCC  (26) ; /•  TNTfPPK  •/ 

IP  V0RD(2)  =PNAFr  tRPN  00;  «1U**P=PIPST ; GO  TO  PR0C26D; 
00  1*1  TO  NUTN;  I*  WOPD  (2)  =UTN,AH  E (I)  THEN  DO; 
H0PP=0LIST (I) ; GO  TO  PP0C26P;  SND;FND; 


PNG; 
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GO 

TO 


TO  r,»TCCn; 
PPOC2^C; 


PB0C26C:  POT  SKIP  lISTpjnP  NOT  FpriND*); 

PP0C26B:  IF  CHAIN?(NS09(I))  =NIIIL  THEM  GO 
HnfP»CHATNf  (MSOP  (T) ) ; 

PP0C260:  DO  1=1  TO  NHIAT;  I?  CHAIN' (I) ^=N0LL  THEM 

PUT  SKIP  EDITfONIFFd)  ,*  9TGK=  • ,CHA  I Np  (i)  - >CA  PK. 
•mTEPBK«  SCAPK'SHALLK  (I)  ) 

CCOl  (2)  , A,  A,»  (P,  3)  ,x  (2)  ,A,P  (8,3)  ) ; FND;  GO  TO  G'TCOH; 
BUT)  HOPCAP;  /•  ♦/ 
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O.B20225.119H0.3TGK.PIT 


/•  CAICOUTB  A IK  noL*^. 
BIGK;  PROC  (?K,Nn’.AT)  R 
DCL  RK  (•)  , IT3RATE 


»CPR  */ 
''THPVS  (F-^OKT) 
LABFL; 


/• 


CIICOLATB  SOM  OP  PK'S  •/ 
SOHK«0;  PO  1*1  TO  KUMAT;  :> 
IP  APS  (SOMK-1 ) < 1 , P-S  TH'*N 
IP  SOPK<1.  THEN  GO  TC  POSK 


nMK*PK(l)  ♦SriMR; 
PPTnpMfT.  FO)  : 


FND; 


ITEPATE=HOnEIM:  GO  TO  TEST; 


/•  -1  < K < 0.  TRY  BK*-.S  */ 

MBGK:  BE*-, 5;  ADJ*-.*5; 

/♦  0<E  , Rr=i.  •/ 

S2!I;  iTrp»ir  = POSM:  no  -o  t»st- 

POSKl;  I'  S|l>Sl  THEN  (50  TO  n0SK2* 

BP«PK40P;  go  to  T'’ST; 

P0SK2:  iTESATE=HCNElH;'AD,i*.2S.BK:  I»  8lt»1  -fH-H  mi-  e 

BK»EK-AD.t;  GO  -O  TEST;  ’ 


MCBEIII:  ADJ*,S*ADJ; 


S?  < ST.  TREV  OK=BE*ADJ; 


/•  EfAtOATF  SIDES  OF  1fP* 
TFST:  SL*1.  ♦ np;  s?*1.; 
SR*SF*(  1.  4RK*pv  (T) 

IP  ABS(SR-Sl)  <1.E-1 
GO  TO  ITFPATE* 

IBD  BIGK; 


PPOD  (UKP  (T)  ) •/ 

DC  1*1  TO  Nnw.'T; 
PVO; 

THE  < PFTGPN  . 


USE  BK*BE-ADJ; 
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0.N20225.  »19a'i.UHI.PlI 


/•  flTS  THF  *0»»!  U(X)=P*  (1-PXPC-CJf))  . 

POMCTIOM  IS  flONOTCNlC  TNCPPXSING  ON 
0»1.  C IS  THE  ?TSK  AVE3STC*!  C^JMSTAN^ 
OHIFIP:  PROC  fXniO,XLJ,THI.B,C)  ; 

OCl  ITR?A‘»'F  LABFL; 

Sl*X»!in; 


TH?  TNTPPVAL 

•/ 


/•  CHRCX  OH  RAHG*  TO  SFAECH  FOR  C •/ 

/•  TPT^|C||"*l"**/*'^  t^SIOMe-l.;  ELSr  CSlr,M=1.; 

C*1,*CSIGN;  ITER AT»=P^NGS?IKO;  GO  TO  TEST; 

PAHGEPIND;  I?  (SP- SL)  •CSIGN<0 , THEN  GO  TC  »ANGE»CffND* 

C«C4C;  GO  TC  TEST; 

HAHGFPOnHC:  IT»P  ATE*H01^  IN ; 

IP  ABS(C)«1  THFH  AD*1*,S*r;  *TSE  Anj*,25*C; 

C*C-AOJ;  XOJa Anj*CS IGN;  GO  To  TEST; 

HCHEIH;  ADJ*.5*Af>J;  I?  SS<SL  THEN  C*C-Anj;  ELSE  C-T^ADJ; 

/•  EPALOATE  T^’ST  POP  C •/ 

TEST:  SPa-LOG (. 5E (FX? (-C-XLO) ♦^XP (-C*XHI) ) ) /C; 

IP  AES  (SR-SL)  <1  . E-3  THEN  GO  To  CUT;  GC  TO  IT'RAT^* 

OOT:  B*  V(1-EXP  (-C)  I ; PETUPN; 

IHD  ONIEXP; 
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APPENDIX  C 

SOME  ALGORITHMS  USED  IN  M/JFCAP 


Apart  from  implementing  the  formula  definitions 
necessary  to  calculate  particular  quantities,  certain  MUFCA? 
routines  make  use  of  some  numerical  analysis  techniques  or 
algorithms.  These  are  discussed  in  the  appendix. 


C«  1 Calculation  of  the  Parameter  k in  the  Multiplicative 

Utility  Function 

A subroutine  called  BIGK  calculates  the  k in  the  mul- 
tiplicative utility  function  using  (3)  described  in  Section 

2.1.  The  algorithm  employed  is  an  iterative  one  suggested  in 

n 

Keeney  (9).  Essentially,  depending  on  the  value  of  E k.  , 

i=l  ^ 

an  interval  is  isolated  where  the  value  of  k must  lie.  Once 

a finite  interval  has  been  found  where  k lies,  the  bisection 

method  for  finding  a real  root  as  described  in  Hamming  [ 5 1 is 

used  to  calculate  k to  the  desired  accuracy, 
n 

When  E k.  > 0 , we  know  -1  < k < 0 and  we  have  our 
i=l  ^ n 

interval  immediately.  When  E k.  < 0 , BIGK  tries  succes- 

i=l  ^ 

sive  powers  of  2 until  a comparison  of  the  two  sides  of  (3) 
indicates  that  a real  root  lies  in  the  interval  (2*^  2^) 

where  n is  as  large  as  necessary  for  the  particular  case.  The 
bisection  method  is  then  used  on  this  interval  to  calculate  k 
to  the  desired  accuracy. 


-U7- 


Hairaning  I 5 ] explains  why  the  bisection  method  is  a 
good  one  to  use  as  opposed  to  other  methods.  Aside  from  being 
easy  to  implement,  it  is  less  vulnerable  to  ill-behavior  and 
round-off  error  than  other  algorithms. 


Calculation  of  the  Constant  Risk  Scalar  Utilit 


Function 


A subroutine  called  UNIEXP  calculates  the  parameter  c 


in  the  constant  risk  form  u(x)  = a + b(l  --  e" 


where  the 


conditions  that  u(0)  =0  and  u(l)  = 1 impose  the  values  a = 0 
and  b = 1/(1  - e ^) . Internally,  MUFCAP  "normalizes"  all 
scalar  attributes  to  run  between  0 and  1.  For  constant  risk 
attributes,  MUFCAP  internally  has  the  attribute  increasing  on 
the  interval  [0,  1).  On  input  and  output,  the  appropriate 
scale  conversions  are  always  made  so  the  internal  normaliza- 
tion is  transparent  to  the  user  except  in  displaying  the  para- 
meters b and  c. 

One  reason  for  normalization  is  that  calculating 
utility  values  using  the  computer's  exponential  algorithm  is 
made  more  accurate  when  the  argument  for  the  exponential  func- 
tion is  not  excessively  large.  This  consideration  is  dis- 
cussed in  Schlaifer  (16). 

UNIEXP  is  very  similar  to  BIGK  in  its  algorithmic 
method.  The  equation  used  is  similar  to  that  in  Schlaifer  [16] 
where  he  discusses  fitting  constant  risk  forms.  Again,  the 
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bisection  method  is  used  because  of  its  nice  "idiot-proof" 
properties. 


C.3 


Calculation  of  Gradient  Comoonents 


The  formula  for  the  quantity  is  derived  in  a 

straightforward  manner  from  either  (1)  or  (2)  in  Section  2.1. 


3u 


The  quantity 


, where  x.  designates  a scalar  attribute 

X X 


amount  and  x is  a "certainty"  alternative,  is  calculated  via 

3u  3u 


the  chain  rule 


3x.  3u . * dx. 

Ill 


Because  of  the  various 
du^ 


forms  possible  for  u.,  the  quantity  is  calculated  by 

X • 

dui  1 

using  the  approximation  - = [u.  (x.  + .01)  - u.(x.)]/ 

QX^  11  11 

[(x^  + .01)  - (x^)].  Remember  (as  explained  in  C.2)  that  inter- 
nally, MUFCAP  scales  all  variables  to  run  between  0 and  1. 

This  approximation  was  felt  to  be  adequate  for  the  purpose  of 
the  program.  When  u^  is  a piecewise  linear  form,  the  expres- 
sion for  the  derivative  when  x^  is  a breakpoint  represents 
the  change  in  the  function  when  moving  in  the  dJ  rection  from 
the  first  range  value  to  the  second  range  value. 
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APPENDIX  D 

MUFCAP*S  OVERALL  PROGRAM  DESIGN 

This  appendix  gives  an  overview  of  the  of  crating 
characteristics  and  programming  design  of  MUFCAP. 

D. 1 Language  and  Operating  System  Considerations 

The  package  is  composed  of  three  procedures  which  are 
compiled  separately  and  then  linked  together.  The  main  proce- 
dures is  called  MUFCAP  and  contains  the  bulk  of  the  package 
making  use  of  internal  procedures  sharing  common  data  bases. 
The  two  external  subroutines  are  BIGK  and  UNIEXP  which  are 
described  in  Appendix  C. 

The  entire  .ackage  is  written  in  PL/1  using  IBM*s  PL/1 
optimizer  compiler.  Features  of  PL/1  which  are  used  heavily 
are  its  based  storage  capabilities  for  managing  linked  lists 
and  its  recursive  function  capabilities  for  dealing  with 
nested  multiattribute  utility  functions.  It  is  conceivable 
that  a MUFCAP  without  nesting  or  a single  level  of  nesting 
could  be  written  in  a language  like  FORTRAN,  but  a more  power- 
ful language  such  as  PL/1  seems  much  more  suitable  for  the 
general  nature  of  this  programming  task.  A helpful  reference 
for  PL/1  is  rike  (21. 

MUFCAP  currently  runs  on  an  IBM  370/165  using  IBM's 
Timesharing  Option,  TSO.  It  runs  in  a partition  of  300K  when 
using  files  for  input  and  output  although  I believe  it  could 
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get  by  with  less  memory.  MUFCAP  stores  information  on  files 
with  a fixed  record  format  of  blocksize  13000  bytes  using  IBM 
3330  disk  drj.vis.  These  file  characteristics  correspond  to  a 
structure  in  the  program  designed  to  have  room  for  roughly 
twenty  scalar  attributes.  These  can  be  adjusted  if  certain 
data  structures  in  the  program  are  made  larger  or  smaller  and 
if  a track  overflow  option  is  used  on  the  IBM  system  for 
blocksizes  larger  than  13000  bytes.  To  create  a dataset  for 
MUFCAP  use,  the  following  TSO  commands  work  for  the  current 
version: 

attrib  trib  recfm(f)  blksize  (13000)  Irecl (13000) 
allocate  file(namel)  dataset (name2)  using (trib) 

* space (5  2)  block (13000) 

The  parameter  'namel*  is  the  name  MUFCAP  uses  in  the 
READ  and  SAVE  commands.  After  a dataset  has  been  created,  new 
datasets  may  be  more  easily  created  by  copying  an  old  one  into 
a new  one  using  the  TSO  COPY  command.  Before  using  MUFCAP, 
all  datasets  which  are  to  be  read  or  saved  should  be  allocated 
using  the  TSO  ALLOCATE  command.  This  is  illustrated  in  Sec- 
tion 5.1.4. 

MUFCAP  ifj  861  cards  long.  Some  estimates  of  relevant 
costs  are: 

compilation  of  program  package  $12  - $15 

linking  the  programs  into  a load  module  $2  - $3 


a one-hour  assessment  and  use  session 


$5 
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D.2  Data  Structures  in  MUFCAP 

There  are  two  central  data  structures  in  MUFCAP;  one  is 
for  MUP*s  and  the  other  for  UNIF's.  For  any  MUF  required 
during  the  program,  a data  structure  is  allocated  with  provi- 
sion for  the  following  information:  the  parameter  k for  the 

function,  an  associated  function  name  and  the  number  of  attri- 
bute arguments  of  the  function.  Each  MUF  has  room  for  12 
attribute  arguments.  For  each  of  these  arguments,  the  MUF 
structure  contains  the  following  information:  a pointer  to 

another  MUF  structure  if  an  attribute  argument  is  a vector,  a 
pointer  to  a UNIF  structure  if  the  attribute  arguments  is  a 
scalar,  the  for  that  attribute  and  the  name  of  that  attri- 
bute. 


When  a scalar  utility  function  or  UNIF  is  required 
during  the  program,  a data  structure  is  allocated  with  provi- 
sion for  the  following  information:  two  range  boundary  values 

for  the  scalar  attribute,  the  utility  function  type,  room  for 
10  attribute  amounts  and  the  utilities  of  those  amounts  for 
"certain”  alternatives,  location  for  up  to  30  parameters  to 
specify  the  utility  function  (e.g.,  15  abscissa  and  ordinate 
values)  and  room  for  5 probabilistic  alternatives  each  denoted 
by  a cumulative  piecewise-linear  distribution  which  may  be 
specified  by  as  many  as  9 points. 

Along  with  these  data  structures  are  three  arrays 
which  contain  the  names  of  ail  the  attributes,  a pointer  to 
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the  MUF  where  the  attribute  is  "located”  and  the  argument 
number  of  the  attribute  in  that  MUF.  By  scanning  these 
arrays,  the  program  finds  the  desired  attribute  name  and  then 
has  pointers  to  all  the  information  necessary  to  perform  cal- 
culations involving  that  attribute  name. 

Data  structures  are  allocated  when  needed  in  a desig- 
nated area  which  can  be  written  out  on  a file  using  the  SAVE 
command.  The  relevant  pointers  are  expressed  as  offsets  to 
the  beginning  of  this  area. 


D.3  Recursive  Functions  and  Nestin 


The  data  structures  and  PL/l*s  recursive  procedure 
capability  enable  the  same  algorithms  to  handle  any  level  of 
nesting.  An  example  will  illustrate  the  point.  Suppose  the 
program  needs  to  evaluate  a MUF.  A routine  is  called  for  this 
purpose  using  (1)  or  (2)  of  Section  2.1  after  a pointer  has 
been  set  pointing  to  the  appropriate  MUF.  Now,  suppose  during 
the  course  of  evaluating  (1)  or  (2),  a vector  attribute  is 
encountered  having  an  associated  MUF  of  its  own.  At  this 
point,  the  routine  merely  saves  the  pointer  to  the  current 
MUF, sets  up  a pointer  to  the  nested  MUF,  calls  itself  to  eval- 
uate the  nested  MUF  and  takes  that  value  and  uses  it  as  it 
resumes  its  previous  calculation.  PL/l*s  recursive  procedure 
capability  handles  all  the  appropriate  bookkeeping.  MUFCAP 
uses  recursive  routines  to  perform  MUr  evaluations,  to 
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calculate  gradients,  to  chain  through  the  multiattribute 
utility  function  structure  in  setting  up  the  three  arrays  men- 
tioned in  Section  D,2  and  in  setting  up  a nested  MUF. 


D.4  Evaluating  Alternatives 


As  explained  in  Section  D.2,  each  UNIF  structure  con- 
tains room  for  specifying  the  scalar  component  for  each  of  the 
various  alternatives.  Whenever  an  alternative  is  specified  or 
a scalar  utility  function  is  set  or  change'^,  MUFCAP  automati- 
cally calculates  the  expected  utility  of  that  scalar  attribute 
for  the  alternative  affected.  By  saving  the  value  of 
E(u^{x^)]  as  well  as  x^,  MUFCAP  saves  a lot  of  redundant  cal- 
culations when  sensitivity  analysis  is  performed  involving 
only  changes  in  the  k^*s.  There  are  separate  routines  for 
calculating  expected  utilities  for  scalar  utility  functions 
depending  on  the  scalar  utility  function  type. 

Various  flags  in  the  program  enable  MUFCAP  to  keep 
track  of  when  it  is  dealing  with  a certain  alternative  or  a 
probabilistic  alternative.  The  names  for  alternatives  are 
contained  in  appropriate  arrays  and  are  saved  when  the  SAVE 
command  is  used. 


Program  Flow 


Program  flow  in  MUFCAP  revolves  around  the  command 


processor  section.  This  section  determines  what  kind  of 
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command  is  requested  and  then  transfers  to  the  appropriate 
command  execution  section.  Aftti  it  is  finished  executing  the 
command,  the  execution  section  transfers  back  to  the  command 
processor  section  for  another  command. 

The  execution  sections  are  not  internal  procedures  but 
invoke  procedures  as  is  necessary.  Operations  which  are 
invoked  by  more  than  one  execution  section  or  are  repeated 
fairly  often  are  incorporated  into  internal  procedures. 
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APPENDIX  E 

TRADEOFF  PROPERTIES  OF  THE  ADDITIVE  AND  MULTIPLICATIVE  FORMS 

Tradeoffs  between  attributes  and  X2  with  the  other 
^^tJ^ibutes  (X^^.-./X^)  held  fixed  can  be  represented  by  an 
indifference  map.  An  indifference  map  is  a set  of  indiffer- 
ence curves  each  having  the  property  that  no  point  on  a par- 
ticular curve  is  preferred  to  any  other  point  on  that  same 
curve.  That  is  to  say,  all  the  points  on  a particular  curve 
are  indifferent  to  each  other.  The  "points”  here  are  conse- 
quences X with  (X2#...#x^)  held  fixed  but  x^^  and  X2  allowed  to 
vary.  An  indifference  curve  is  generated  when  we  choose  a 
pair  (x^,  X2)  and  display  all  the  allowable  X2)  pairs 

which  are  indifferent  to  it. 

When  the  requisite  assumptions  to  imply  either  (1)  or 
(2)  are  satisfied  (Section  2.1),  an  indifference  curve  is 
represented  analytically  by  (x^,  X2)  pairs  satisfying 

(Xi)  + kk^k2Uj^  (x^)u2  (X2)  = constant  (E-1) 

This  equation  results  from  the  fact  that  when  two  consequences 
X*  and  x”  are  indifferent,  u(x*)  = u(x”).  When  k = 0 in  (E-1), 
this  corresponds  to  the  additive  form.  When  k ^ 0,  this  cor- 
responds to  the  multiplicative  form. 

From  (E-1)  we  can  see  that  (x^,  x^)  pairs  which  are 
indifferent  to  each  other  remain  indifferent  regardless  of 
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the  level  at  which  (x^/.-./X^)  happen  to  be  fixed.  Suppose  we 
wished  to  generate  an  indifference  curve  using  only  tradeoff 
information  between  and  X2.  Since  k in  general  depends  on 
the  other  k^  via  (3)  (Section  2.1),  we  can  generate  two  inde- 
pendent equations  using  two  sets  of  indifference  pairs  varying 
x^  and  X2«  Using  these,  wc  can  express  k and  k2  in  terms  of 
k^.  Setting  k^^  to  an  arbitrary  number  corresponds  to  setting 
the  constant  on  the  right  hand  side  of  (E-1)  to  an  arbitrary 
constant.  This  does  not  affect  which  points  are  indifferent 
to  each  other.  Thus,  two  sets  of  indifference  pairs  which  are 
independent  enables  us  to  calculate  the  parameters  of  an  equa- 
tion for  indifference  curves.  Then,  if  we  are  given  any  point 
(x^>  ^2^  • generate  all  the  {x^^,  X2)  pairs  which  are 

indifferent  to  it.  To  summarize,  indifference  curves  repre- 
senting tradeoffs  between  X^  and  X2  can  be  generated  using 
only  information  concerning  preferences  over  (x^,  X2)  pairs 
and  need  not  require  any  specific  tradeoff  information  con- 
cerning the  other  attributes. 

If  we  let  y^  = u^^(x^)  and  y2  = U2(x2),  equation  (E-1) 

becomes 

^1^1  ^2^2  ^^1^2^1^2  “ constant 

An  indifference  curve  in  (y^^,  y2)  space  as  opposed  to 
X2)  space  is  always  a hyperbola. 
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Now  let  us  examine  the  effect  of  nesting  on  indiffer- 
ence curves.  We  will  examine  a three-attribute  case  of  the 
form  u = u(u^,  u^^)  where  u^  = u^(a)  and  Uj^  = (u^,  u^)  . Thus, 
the  three  single  attributes  involved  are  A,  S and  T. 

In  the  multiplicative  form,  we  have,  symbolically 
(where  the  arguments  of  the  utility  functions  have  been  left 
out  for  more  concise  notation) , 

1 + ku  = (1  + kk^u^) (1  + kk.  u,  ) (E-2) 

a a DO 

3 + k*Uj^  = (1  + k*kgUg)  (1  + k'k^u^)  (E-3) 

Substituting  (E-3)  into  (E-2)  yields 

1 + ku 

= (1  + kk^u^)(l  + kk,/k*[(l  + k»k  u^)  (1  + k»k.u.)-l])  (E-4) 

a a D s s t:  t: 

Now,  note  what  happens  if  k*  = kkj^ 

We  then  obtain 

1 + ku 

= (1  + kk^u^)  iX  + ((1  + kk.k^u^)  (1  + kk.k.u.  )-ZJ) 
aa  OSS  Dtt 

= (1  + kk  u Ml  +kk'u,M  1 + kk'  u^)  (E-5) 

a a s s t t 

where  k*  = k.  k 

S D S 

k*  = k k 

Equation  (E-5)  is  nothing  but  the  multiplicative  form 
for  three  attributes.  Thus,  if  k*  - kkj^,  any  pair  of 
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attributes  has  the  preferential  independence  property  and  the 

indif f crenc-'  cur*’-  rrcpr.rties  cf  (E-1)  apply.  However,  if 

K ' kk^,  this  is  no  longer  true.  We  can  no  longer  factor 

thr  expression  for  1 ku  into  three  factors  each  dealing  wit 

a single  attribute.  Because  of  this,  if  u(a*,s*,t)  = 

r.  ^“",3”,t),  it  IS  not  necessarily  the  case  that  u(a',s',t')  “ 

li  (a"  , s”  . t * ) where  t*  t.  That  is  to  say,  indifference  curv^ 

^ -tvr.cn  a and  s depend  on  t when  there  is  nesting  and 

• ^ kk.  . 

b 

MUFCAP  has  a corjnand  IN’TERBK  which  calculates  the 

quantity  kkj^  and  corrpares  it  to  k*  where  b is  any  vector 

attribute  in  a particular  MUF  and  k,  , k and  k*  are  the  analo- 

b 

st/s  parar.etcrs  to  those  in  our  cxar.plc.  If  kkj^wk*..  then 
the  nesting  of  attributes  into  their  own  internal  MUF  may  be 
unnecessary.  Section  5.1.4  has  an  illustration  of  the  use  of 


i::terbk 


